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Abstract 

A  new  approach  to  calculating  nucleon-nucleon  scattering  matrix  elements  using 
a  proven  atomic  time-dependent  wave  packet  technique  is  investigated.  Using  this 
technique,  reactant  and  product  wave  packets  containing  centripetal  barrier  information 
are  prepared  in  close  proximity  to  the  nuclear  potential  energy  well.  This  is 
accomplished  by  first  using  an  analytic  equation  to  determine  the  wave  packets  in  a 
suitable  intermediate  asymptotic  state  where  the  centripetal  barrier  is  negligible.  Then, 
the  split  operator  technique  is  used  to  propagate  the  wave  packets  back  to  their  original 
positions  under  the  full  Hamiltonian.  Here,  the  product  wave  packet  is  then  held 
stationary  while  the  reactant  wave  packet  is  allowed  to  evolve  and  explore  the  nuclear 
well.  Scattering  matrix  elements  are  computed  from  the  correlation  function  between  the 
stationary  wave  packet  and  the  evolving  wave-packet  after  it  has  interacted  with  the 
nuclear  potential.  Determination  of  nucleon-nucleon  phase  shifts  follows  directly  from 
computation  of  the  scattering  matrix  elements.  This  technique  is  ideally  suited  for 
determining  nuclear  scattering  matrix  elements  and  phase  shifts  as  it  provides  a  high 
degree  of  energy  resolution  with  lower  computational  effort  than  traditional  time 
independent  methods.  These  advantages  will  lead  to  a  greater  understanding  of  nuclear 
reaction  dynamics. 
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TIME  DEPENDENT  CHANNEL  PACKET  CALCULATION 


OF  TWO  NUCLEON  SCATTERING  MATRIX  ELEMENTS 

1.  Introduction 


1.1.  Motivation 

In  1953,  Hans  Bethe  published  an  article  in  Scientific  American  where  he 
discussed  the  force  that  holds  the  nucleus  together.  In  this  article,  he  makes  the 
following  comment. 

“/n  the  past  quarter  century  physicists  have  devoted  a  huge  amount  of 
experimentation  and  mental  labor  to  this  problem  -probably  more  man¬ 
hours  than  have  been  given  to  any  other  scientific  question  in  the  history 
of  mankind  [1].” 

Over  the  50  years  that  have  passed  since  Hans  Bethe  penned  this  article  for  Scientific 
American,  not  much  has  changed.  Many  of  the  world’s  greatest  minds  are  still  engaged 
performing  experiments,  examining  data,  and  developing  models  in  an  attempt  to 
characterize  the  nuclear  strong  force  that  binds  nuclei  together.  Currently,  a  first 
principles  description  of  even  the  two-body  nucleon  interaction  does  not  exist.  Quantum 
chromodynamics,  our  closest  first  principles  approach,  performs  well  at  high  Giga- 
Electron  Volt  (GeV)  energies  but  has  not  been  extended  to  the  lower  energies  for 
computational  reasons  [2,3]. 
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In  lieu  of  a  first  principles  approach,  the  nuclear  community  has  developed 
phenomenological  nucleon  -  nucleon  (NN)  models  to  describe  observed  features  of  the 
strong  force.  Early  attempts  to  describe  the  nuclear  interaction  phenomenologically,  such 
as  the  one  by  Gammel  and  Thaller  in  1957,  were  only  able  to  recreate  aspects  of  the 
experimental  data  set.  The  experimental  data  set  consists  primarily  of  scattering  matrix 
elements  derived  from  observables.  Gammel  and  Thaller’s  as  well  as  all  subsequent 
phenomenological  model’s  are  all  based  on  parameterized  fits  to  experimental  data  [4- 
31].  Of  the  more  recent  models,  three  in  particular,  have  been  able  to  reproduce  almost 
all  of  the  features  of  the  experimental  two-nucleon  data  set:  the  high  precision  models  of 
Nijmegen,  Bonn,  and  Argonne  National  Labs  (ANL)  [2,  4,  32,  33]. 

A  numerical  technique  is  required  to  calculate  scattering  matrix  elements  from 

any  modeled  potential  surface  for  comparison  to  experimental  data.  In  order  to  compute 

scattering  matrix  elements,  the  nuclear  community  has  traditionally  relied  on  techniques 

based  on  solutions  to  the  time-independent  form  of  Schrodinger’s  equation, 

y/"  =  k^y/  (1.1) 

where  k  is  the  wave  vector.  Indeed  with  the  exception  of  some  time-dependent 

momentum  space  calculations  performed  by  Holz  and  Glockle  [12,34]  and  a  couple  of 

deuteron  breakup  calculations  in  the  80s’  [12],  almost  all  nuclear  scattering  matrix 

elements  have  been  computed  via  time  independent  means.  In  the  time-independent 

approach,  the  coordinate  or  momentum  space  approximations  to  the  wave  function,  y/ , 

provide  s-matrix  elements  at  a  single  energy  for  a  variety  of  asymptotic  quantum 

numbers.  Although  many  techniques  exist  to  solve  a  second  order  linear  differential 

equation  such  as  the  Schrodinger  equation,  the  one  that  is  utilized  typically  by  the  nuclear 
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community  is  the  Numerov  method.  In  this  approach,  a  Taylor  series  expansion  of  the 
wave  funetion  is  used  to  develop  a  three-term  reeurrenee  relation.  From  this  reeurrenee 
relation,  seattering  matrix  elements  are  then  eomputed  by  matehing  boundary  eonditions 
where  solutions  are  assumed  to  be  linear  combinations  of  free  space  solutions. 

While  the  nuelear  community  has  focused  primarily  on  eomputing  seattering 
matrix  elements  via  time-independent  methods,  the  atomie  and  moleeular  eommunity  has 
researched,  developed,  and  implemented  time-dependent  teehniques  as  an  alternative 
approach  for  calculating  atomic  reaction  probability.  References  to  time-dependent 
seattering  theory  first  started  to  appear  in  the  atomie  and  moleeular  literature  as  early  as 
1956  [35].  However,  it  was  not  until  the  early  80 ’s  that  the  eomputing  power  beeame 
available  for  effeetive  employment  of  time-dependent  algorithms.  Sinee  then,  these 
teehniques  have  made  signifieant  eontributions  to  our  understanding  of  few-body  atomie 
interaetions  [36-50]. 

Time-dependent  approaches  rely  on  propagation  of  wave  paekets  via 
approximations  to  the  time  evolution  operator, 

=  (1.2) 

where  h  is  the  redueed  Planek’s  eonstant  (a.k.a.  Dirae’s  eonstant)  and  His  the 
Hamiltonian.  In  eontrast  to  the  previously  mentioned  time-independent  methods,  whieh 
provide  s-matrix  elements  at  a  single  energy  for  a  variety  of  asymptotie  quantum 
numbers,  time  dependent  methods  provide  s-matrix  elements  for  a  single  set  of  quantum 
numbers  as  a  funetion  of  energy  in  a  single  ealculation. 
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One  time-dependent  method  that  has  been  sueeessfully  applied  to  a  variety  of 
atomie  and  moleeular  problems  is  ealled  the  Channel  Paeket  Method  (CPM)  [36-40,  52], 
The  CPM  allows  for  ealculation  of  a  single  seattering  matrix  element  in  eontrast  to  time- 
independent  teehniques  that  eompute  the  entire  column  of  the  scattering  matrix.  Since 
typically  only  one  scattering  matrix  element  for  a  range  of  energies  is  desired  in  a  nuclear 
scattering  calculation,  the  CPM  is  an  ideal  choice  for  nuclear  scattering  calculations. 

CPM  calculations  begin  by  preparing  two  complex  Gaussian  wave  packets  over  a 
range  of  energies  for  which  we  wish  to  determine  s-matrix  elements.  One  of  these 
contains  negative  momentum  and  is  designated  as  the  reactant  state.  The  other  contains 
positive  momentum  and  is  designated  as  the  product  state.  These  two  wave-packets  are 
then  typically  propagated  using  the  split-operator  method  first  away  from  the  interaction 
region  under  the  asymptotic  Hamiltonian  then  back  to  their  original  positions  under  the 
full  Hamiltonian.  At  this  point,  the  two  wave-packets  are  handled  differently.  The 
product  wave-packet  is  held  stationary  whereas  its  reactant  state  counterpart  is  allowed  to 
continue  to  evolve  and  probe  the  characteristics  of  the  nuclear  potential.  As  the  reactant 
wave  packet  leaves  the  interaction  region,  its  correlation  with  the  product  wave-packet  is 
assessed  to  determine  the  scattering  matrix  elements  between  the  initial  and  final  states. 

In  this  dissertation,  the  CPM  technique  is  used  to  compute  nucleon-nucleon 

scattering  matrix  elements  for  proton-proton,  proton-neutron,  and  neutron-neutron 

scattering  events.  Of  course,  for  this  application  of  the  CPM  technique  to  be  successful,  a 

complete  understanding  of  the  potential  energy  basis  and  form  of  the  Hamiltonian  is 

necessary.  Since  ANL’s  model  has  enough  supporting  documentation  to  give  insight  into 

a  natural  choice  for  potential  energy  basis  and  Hamiltonian,  this  research  effort  uses  the 
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AVI 8  model  to  determine  two-body  nuelear  seattering  matrix  elements.  This  novel 
approaeh  to  performing  nuelear  seattering  ealeulations  not  only  provides  a  signilieant 
improvement  in  eomputational  effieieney  over  existing  time  independent  teehniques,  it 
also  provides  the  nuelear  eommunity  with  intuitive  tool  for  visualizing  nuelear  reaction 
dynamics.  These  characteristics  will  be  invaluable  to  future  scattering  research  efforts  in 
nucleon-deuteron  scattering  and  muon  catalyzed  reactions. 

1.2.  Overview 

This  research  effort  examined  the  suitability  of  using  the  time-dependent  channel 
packet  method  to  calculate  two  nucleon  scattering  matrix  elements.  The  effort  was 
broken  into  two  phases.  Phase  1  evaluated  the  outputs  and  failure  modes  of  the  AVI  8 
FORTRAN  subroutine  to  gain  insight  into  a  suitable  choice  for  the  coordinate  system, 
Hamiltonian  and  basis  set.  Phase  2  assessed  how  well  the  CPM  technique  could  be 
adapted  to  the  nuclear  scattering  problem.  A  FORTRAN  CPM  code  was  developed  and 
validated  against  an  analytic  solution  to  a  square  well  problem  of  similar  dimensions  as 
the  nuclear  well.  Scattering  matrix  elements  were  then  calculated  from  the  AVI 8 
potential  surfaces. 

The  significant  research  accomplishments  presented  in  this  document  are: 

1 .  Establishment  of  new  and  more  efficient,  time-dependent  nuclear  technique, 
which  is  ideally  suited  for  determining  nuclear  scattering  matrix  elements. 
(The  numerical  effort  associated  with  the  Time  Dependent  Channel  Packet 
Method  scales  as  N  or  better  whereas  the  time-independent  methods  currently 
employed  by  the  nuclear  community  scale  as  N^.) ^  [50] 


*  Here,  N  refers  to  the  number  of  grid  points  in  the  calculation. 
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2.  Provides  the  nuelear  eommunity  with  a  unique  intuitive  tool  for  visualizing 
seattering  dynamies  in  both  momentum  and  eoordinate  representations. 

3.  Provides  a  eomplete  set  of  nuelear  seattering  information  for  a  speeified  value 
of  angular  momentum  with  a  signifieantly  higher  degree  of  energy  resolution 
that  is  provided  by  the  eurrent  time  independent  methods  employed  by  the 
nuelear  eommunity. 

This  doeument  is  organized  as  follows.  Chapter  2  outlines  how  seattering  matrix 
elements  are  obtained  from  experimental  observables.  Chapter  3  foeuses  on  the  AVI 8 
potential  and  the  time-independent  ealeulation  of  seattering  matrix  elements.  A  brief 
exploration  into  the  physieal  basis  behind  the  One  Pion  Exehange  and  the  strong  foree 
parameterization  is  ineluded.  In  Chapter  4,  the  theoretieal  baekground  behind  using  the 
CPM  to  ealeulate  S-Matrix  elements  is  deseribed.  Later  in  Chapter  5,  eoneems 
pertaining  to  using  the  CPM  teehnique  to  ealeulate  seattering  matrix  elements  from  the 
AVI 8  potential  surfaees  are  addressed.  Topies  reviewed  here  inelude  referenee  frame, 
basis,  Hamiltonian,  1/r  cutoff  correction,  units,  coupled  basis  diagonalization  and  s- 
matrix  parameterizations.  Also  presented  in  this  section  is  a  comparison  between  a  time- 
independent  and  a  time  dependent  solution  to  a  1-D  Square  Well  problem.  Then  in 
Chapter  6,  the  phase-shifts  obtained  from  the  nuclear  CPM  calculations  are  presented  and 
compared  to  ANL’s  published  results.  Chapter  7  reviews  the  technical  achievements  of 
this  research  effort  and  presents  possible  avenues  for  future  research.  The  appendices 
contain  a  complete  set  of  CPM  nuclear  phase  shifts  through  J  =  5 
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2.  Experimental  Background 


This  chapter  provides  a  heuristie  diseussion  of  how  S-Matrix  elements  are 
ealculated  from  experiment.  It  opens  with  a  general  discussion  of  why  examining  nuclear 
structure  requires  development  of  high-energy  partiele  aecelerators.  The  aceeleration 
methods  discussed  here  are  limited  to  the  basic  theory  of  charged  particle  acceleration 
and  the  leading  teehniques  of  producing  high-energy  neutron  beams.  This  seetion  then 
eoneludes  with  a  discussion  of  how  scattering  matrix  elements  are  obtained  from 
scattering  high-energy  polarized  particle  beams  off  target  nuelei. 

2.1.  Electrostatic  Accelerators 

Based  on  Einstein’s  1905  explanation  of  light’s  particle  like  behavior,  De  Broglie 
theorized  that  particles  analogously  might  also  exhibit  a  wave  like  behavior.  His 
argument  was  simple.  If  photons  have  momentum  like  particles,  should  not  particles 

h 

have  a  wavelength,  =  —  similar  to  photons?  His  proof  eame  in  1927  when  Davisson 

P 

and  Germer  observed  an  angle  dependent  interference  pattern  from  the  scatter  of  an 
electron  beam  off  a  Ni  erystal  target.  It  was  learned  that  in  order  to  obtain  information 
about  atomie  structure,  the  De  Broglie  wavelength  ( )  must  be  roughly  on  the  same 

order  or  smaller  than  the  dimensions  of  the  target  under  study.  As  Table  2.1  indieates, 
obtaining  information  about  underlying  nuelear  structure  and  the  short-range  nuclear 
foree  (range  1-2  fm)  that  binds  nueleons  together  requires  development  of  particle 
accelerators  capable  of  aecelerating  particles  to  MeV  energies. 
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Table  2.1.  The  De  Broglie  Wavelengths  of  the  Photon,  the  Electron,  the  Proton,  and  the  Nentron. 


De  Broglie  Wavelength,  fm 

Energy 

Photon 

Electron 

Proton 

Neutron 

0.1  MeV 

mam 

\MeV 

2.9x10^ 

WMeV 

■HH 

9.0x10' 

100  MeF 

1.2x10' 

m 

2.8x10' 

2.8x10' 

Early  accelerators  sueh  as  the  one  used  by  Davisson  and  Germer  to  prove  De 
Broglie’s  theory  were  quite  simple  (Figure  2.1).  Eleetrons  were  aeeelerated  from  a 
grounded  heated  fdament  loeated  at  the  eathode  aeross  a  gap  toward  a  positive  anode. 
Instead  of  impinging  upon  the  faee  of  the  anode,  a  small  hole  in  the  anode  allows  the 
partiele  beam  to  pass  and  interact  with  an  intended  target.  Eleetrostatie  generators, 
however,  were  unreliable  and  prone  to  eleetrostatie  breakdown  in  the  KeV  energy  range, 
well  below  the  MeV  energy  range  required  to  perform  nuelear  seattering  experiments. 
The  inability  of  eleetrostatie  generators  to  provide  suitable  partiele  beams  for  nuelear 
scattering  experiments  spawned  development  of  two  new  types  of  aecelerators  still  used 
today,  the  linear  aeeelerator  and  the  eyelotron  [53-56]. 
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Figure  2.1.  Schematic  Diagram  of  Electrostatic  Accelerator.  A  grounded  heated  filament  accelerates 
electrons  across  a  gap  toward  a  positively  charged  anode.  A  gap  in  the  anode  allows  the  electron 
beam  to  pass  on  to  an  intended  target 
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2.2.  The  Linear  A ccelera  tor 


The  basic  concept  of  linear  particle  acceleration  was  developed  by  Ising  and 
improved  upon  by  Wideroe,  Sloan,  Lawrence,  and  Alvarez  [57].  As  previously 
mentioned,  electrostatic  accelerators  have  problems  exceeding  breakdown  voltage.  The 
linear  accelerator,  or  linac,  overcomes  this  by  accelerating  a  particle  in  a  series  of  stages 
(Figure  2.2). 


Drift  Tubes 


Figure  2.2.  Alvarez’s  Basic  Linac  Concept.  A  pulsed  waveform  is  applied  to  each  cylinder  (drift 
tube)  to  accelerate  particles  across  each  gap.  The  drift  tubes  also  shield  the  particle  beam  from 
decelerating  effects.  [57] 


Here,  a  charged  particle  emitted  from  an  ion  source  passes  through  a  series  of  hollow 
cylindrical  electrodes  called  drift  tubes.  As  the  particle  passes  across  each  gap,  it  is 
accelerated  by  a  pulsed  waveform  applied  to  each  tube  from  a  high  power  radio 
frequency  oscillator.  Each  of  the  drift  tubes  also  functions  to  shield  the  particle  from  a 
variety  of  decelerating  effects.  The  beauty  of  this  approach  is  that  higher  particle  beam 
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energies  ean  be  aehieved  by  simply  adding  more  stages.  Eleetron  linacs  sueh  as  the  two- 
mile  Stanford  Linear  Aeeelerator  Center  (SLAC)  have  been  able  to  aehieve  eleetron 
energies  in  the  50  GeV  range.  Proton  linacs  have  been  able  to  achieve  a  maximum 
energy  of  about  a  GeV  [57]. 

2.3.  The  Cyclotron 

Extending  the  concept  of  the  drift  tube,  Ernest  Lawrence  used  a  magnetic  field  to 
reduce  the  size  of  the  accelerator.  His  idea  was  to  manufacture  a  pair  of  drift  tubes  in  the 
form  of  a  bisected  compressed  hollow  sphere  (Eigure  2.3). 


Figure  2.3.  Basic  Cyclotron  Schematic.  Particles  injected  between  two  strong  magnets  are  forced  to 
follow  a  circular  path  and  are  accelerated  by  an  rf  oscillator  across  the  gap  between  the  magnets. 
[57] 
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These  “dee’s”,  as  he  ealled  them,  were  plaeed  between  two  strong  magnets  whieh  forces 
the  particles  to  follow  a  circular  path.  An  RF  oscillator  matched  to  the  particles  exact 
circulation  frequency  excites  the  cyclotron  drift  tubes  and  ensures  the  sign  of  acceleration 
is  always  in  the  same  direction  when  the  particles  hit  the  gap.  In  this  manner,  cyclotrons 
have  been  able  to  produce  proton  beams  up  to  25  MeV.  Above  25  MeF however, 
relativistic  effects  become  significant.  The  relativistic  mass  increase  requires  more 
energy  to  maintain  the  specified  velocity,  which  results  in  a  drop  in  revolution  frequency 
and  a  loss  of  synchronization  with  the  accelerating  potential.  The  synchrotron  maintains 
synchronization  by  systematically  varying  the  cyclotron  RF  frequency  and  magnetic 
field.  Synchrotrons,  such  as  CERN’s  Super  Proton  Synchrotron  have  obtained  energies 
in  excess  of  400  GeV  [57]. 

2.4.  Production  of  Polarized  Neutron  Beams 

Having  introduced  basic  accelerator  theory  for  charged  particles,  the  question 
remains  how  to  produce  a  polarized  neutron  beam  given  that  neutrons  do  not  have  charge 
and  cannot  be  accelerated  by  any  of  the  traditional  means  identified  in  the  last  section.  A 
reasonable  first  choice  as  a  source  of  neutrons  is  a  nuclear  reactor  since  copious  free 
neutron  production  is  a  byproduct  of  the  fission  process.  Reactors  however  can  only 
produce  neutrons  with  a  maximum  energy  of  about  1  MeV,  far  too  low  for  nucleon- 
nucleon  scattering  experiments  which  typically  require  neutrons  ranging  from  \MeV up 
to  a  GeV  [57].  This  section  summarizes  two  of  the  leading  techniques  of  producing  a 
polarized  neutron  beam  with  a  higher  kinetic  energy  than  a  reactor  can  generate. 
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The  first  method  of  generating  a  high-energy  neutron  beam  is  through  the  breakup 
of  polarized  deuterons^.  In  this  technique,  scatter  of  a  polarized  deuteron  beam  off  a 
beryllium  target  induces  dissociation  and  generates  a  polarized  proton-neutron  beam. 
Passing  this  beam  through  magnetic  fields  generated  by  spin  rotating  solenoids  and 
processing  magnets  not  only  strips  away  the  residual  charged  particles  but  also  provides  a 

-5 

method  of  orienting  the  neutron  beam  along  the  vertical,  horizontal,  or  longitudinal  axis  . 

Since  neutron  beams  generated  by  deuteron  breakup  have  a  high  flux  (25x10^  — 

s  cm 

polarization  (60  %),  and  small  energy  spread  (FWHM  ~  50  MeV),  this  is  considered  the 
preferred  method  of  neutron  beam  generation  [58]. 

Few  particle  accelerators,  however,  are  capable  of  accelerating  deuterons,  so 
proton  -  neutron  exchange  serves  as  an  alternative  method.  In  this  technique,  a  stream  of 
polarized  neutrons  is  produced  by  scattering  an  incident  polarized  proton  beam  off  a 
liquid  deuterium  target.  Here,  the  polarization,  flux,  and  energy  spread  of  the  neutrons 
depend  on  the  scattering  angle.  If  the  proton  beam  has  a  180-degree  angle  of  incidence, 
the  resultant  neutron  beam  has  optimal  energy,  intensity  and  minimal  spread,  but  beam 
polarization  is  only  -40%.  Neutron  beam  polarizations  upwards  of  60%  can  be  obtained 
if  the  incident  proton  beam  has  a  160-degree  angle  of  incidence  and  a  polarization  normal 
to  the  scattering  plane.  If  polarized  neutron  beam  with  energies  greater  than  800  MeV  are 

^  A  polarized  deuteron  has  both  the  proton  and  neutron  spin  components  aligned  in  the  same  direction. 

^  Although  the  neutron  is  electrically  neutral,  it  does  have  a  magnetic  moment  of  -1.913  magnetons.  (The 
proton  magnetic  moment  is  2.973  magnetons.)  This  is  one  of  the  leading  indicators  of  underlying  nuclear 
structure. 
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desired,  the  polarized  proton  beam  is  seattered  off  a  light  nuelei  target  sueh  as  Li,  Be,  B, 
or  C  instead  of  a  liquid  deuterium  target.  Seatter  of  protons  beams  off  these  targets  ean 
produee  a  50%-60%  polarized  GeVheam  with  a  peak  neutron  flux  of  approximately 

50x10^  — .  Table  2.2  summaries  these  methods  of  produeing  high-energy  neutron 
s  cm 

beams  [58]. 


Table  2.2.  Properties  of  the  Neutron  Beams  Produced  by  Different  Neutron  Beam  Productions 
Methods  [58] 


Production  Method 

Liquid  Deuterium 

Deuteron  Breakup 

Carbon 

Primary  Beam  Intensity 

0.25-5.0 

3x10^'  (J  T ! spill) 

2-\QpA{p'[) 

Neutron  Flux  [lO  n  /  (^5 

0.9-5.0 

25 

10-50 

Neutron  Energy  MeV 

180-788 

300-1150 

200-580 

Neutron  Polarization  |%] 

40-60 

59 

40-50 

FWHM  Resolution  MeV 

15-20 

40-60 

11-50 

2.5.  Experimental  S-Matrix  Elements 

As  we  shall  see  in  Section  5,  it  is  rather  easy  to  obtain  cross  sections  or  other 
observables  given  the  S-matrix  elements.  The  converse  is  not  true  as  it  is  quite  difficult 
to  construct  the  S-matrix  given  observable  information.  Here,  an  outline  of  the  process 
by  which  experimental  S-matrix  elements  are  determined  from  observables  is 
summarized  [58]. 
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To  determine  nucleon-nueleon  seattering  matrix  elements  from  observables 
single,  double,  and  triple  scattering  parameters  are  required.  The  measurement  of  the 
single  scattering  parameter  is  illustrated  in  Figure  2.4. 


Figure  2.4.  First  Scattering  Parameter  -  Differential  Cross  Section.  A  Circular  Detector  Array  is 
Positioned  Around  a  Target  to  Determine  the  Angular  Distribution  of  Scattered  Particles  from  an 
Unpolarized  source. 


Here,  an  incident  unpolarized  nucleon  beam  is  scattered  by  a  target  located  in  the  center 
of  a  nucleon  detection  ring.  Detectors'^  placed  at  discrete  angles  and  fixed  distance,  r, 
record  the  scattered  beam  intensity  5,.  as  a  function  of  angle,  0 ,  and  in  most  cases  as  a 


4 


Physical  limitations  on  detector  size  limit  the  amount  of  angular  resolution. 
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function  of  azimuth,  (j) .  The  number  of  partieles  seattered  into  a  solid  angle  dQ.  at 


angle  0  divided  by  the  ineident  flux  is  defined  as  the  differential  eross  seetion 


d<7 

dQ. 


[30], 


2 

dQ  ~  S.  ^ 


(2.1) 


In  most  oases,  the  first  soatter  also  polarizes  the  randomly  oriented  spins  of  the 
ineident  nuoleon  beam.  The  amount  of  polarization  of  the  nuoleon  beam  resulting  from 
the  first  soatter  is  the  seoond  scattering  parameter.  The  Polarization  observable,  P, 
seattered  at  an  angle  0  is  defined  as. 


(2.2) 


where  is  the  number  of  partieles  with  spin  up  and  is  the  number  of  particles  with 

spin  down.  The  amount  of  polarization  is  scattering  angle,  energy,  and  particle 
dependent. 

The  triple  scattering  parameter  oharacterizes  how  a  seoond  soatter  alters  the 
polarized  beam’s  direotion  and  magnitude  [28].  The  orientation  of  the  ineident  beam 
polarization  with  respeot  to  the  secondary  scattering  plane  determines  the  experimental 
quantity  to  be  measured.  For  example,  the  depolarization  observable,  D,  refieots  the 
amount  of  polarization  remaining  in  the  beam  perpendicular  to  the  soattering  plane  after 
the  seoond  soatter.  The  depolarization  triple  soattering  parameter  measurement  is 
represented  diagrammatioally  in  Figure  2.5  [28,  31]. 
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Incident 

Unpolarized 


Figure  2.5.  Diagrammatic  Representation  of  the  Triple  Scattering  Parameter  -  Depolarization,  D. 
Angles  1  and  2  denote  the  first  and  second  scatter.  The  arrows  indicate  the  direction  of  beam  travel 
and  the  circles  indicate  polarization  normal  to  the  page.  [28,  31] 


The  triple  scattering  parameters  {D,  A,  K,  M,  N),  double  scattering  parameter  (P), 
and  single  scattering  parameter  ( )  can  be  expressed  as  combinations  of  five  invariant 

do. 

amplitudes^  a,  b,  c,  d,  and  e  and  the  scattering  angles  from  each  stage  0^ .  Here,  A  refers 
to  the  spin  correlation,  K  refers  to  the  polarization  transfer,  M  refers  to  the  scattered 
particle  polarization,  and  N  refers  to  the  recoil  particle  polarization.  Some  of  the 
equations  used  in  the  determination  of  the  invariant  amplitudes  are  presented  in  Table  2.3 
[58].  In  the  table,  subscripts  on  the  observables  label  the  polarization  of  the  scattered  (5), 
recoil  (r),  beam  {b),  and  target  {t)  particles  and  the  primes  on  the  subscripts  distinguish 
between  initial  {k),  scattered  {k'),  and  recoil  ( A: " )  basis.  These  equations  and  others 

^  These  five  amplitudes  are  invariant  with  respect  to  parity  conservation,  time  reversal  invariance,  the  Pauli 
principle  and  isospin  invariance  [31].  Parity  conservation,  time  reversal  invariance  and  isospin  invariance 
are  discussed  in  Section  3.1.3 
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available  in  the  literature  may  be  eombined  to  eliminate  unknowns  and  solve  for  one  of 
the  five  invariant  amplitudes.  It  is  these  invariant  amplitudes  that  define  the  experimental 
Scattering  Matrix,  S 


S{k,k')  =  +  b)  +  {a  - b){<j^»n){<j^»n)  +  {c  +  d){<j^»m){<j2»m)  + 

(c  -  d){(j^  •/  ){(J2  •/ )  +  e(crj  +  cTj  )•«)] 


(2.3) 


Here,  crj(cr2)  refer  to  the  usual  Pauli  Spin  Matrices  and  the  center  of  mass  basis  vectors 


n,  m,  and  I  reflect  the  orientation  of  scattered  and  incident  particles  [58].  Typically, 
results  from  nuclear  scattering  matrix  calculations  are  presented  in  the  form  of  a  phase 
shift.  How  phase  shifts  are  calculated  from  scattering  matrix  elements  is  discussed  in 
Section  4.  Here,  a  Center  of  Nuclear  Studies  (CNS)  fit  (solid  line)  to  experimental  data 
(error  bars)  is  shown  for  the  Pq  proton  neutron  scattering  matrix  element  (Figure  2.6). 


Figure  2.6.  Sample  CNS  Fit  to  Experimental  ^Po  Phase  Shift  Data  for  a  Proton-Neutron  Scattering 
Event.  The  error  bars  are  experimental  results.  The  solid  line  is  a  fit  to  the  observed  data. 
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Table  2.3.  Some  Scattering  Observables  with  0,1,  2  and  3  Spin  Indices.  Here  0  labels  tbe  center  of 
mass  scattering  angle,  0^  labels  tbe  laboratory  angle  of  tbe  scattered  particle,  and  labels  tbe 

laboratory  angle  of  tbe  recoil  particle.  Observables  may  be  combined  to  eliminate  nnknowns,  solve 
for  tbe  invariant  amplitndes  (a,  b,  c,  d,  or  e)  and  constrnct  tbe  Scattering  Matrix,  S  [58] 


.  Designation 

Definition 

I 

Differential  Cross  Section 

1  ll  |2  It  |2  1  |2  1  ,|2  1  |2'| 

cr  =  — ||a  +  o|  +|c  +a  +e  j 

P 

Polarization 

^Pnooo  =  Re(a)  •  e 

A 

Spin  Correlation 

^Ao„o  =|{H"  “H"  “H" +H"} 

^Aoss  -  ^  cos(^)  +  Re(/))  •  c  -  Im(<i)  •  esin(0) 

D 

Depolarization  tensor 

^Dnono =|{h' +H'  “kr  +kr} 

cDs-oso  =  Re(a)  ■bcos{6  -  6^)  +  Re(c)  •  d  sin(0j)  -  Im(/))  ■ecos{6  -  6^) 

K 

Polarization  Transfer 

^Knno =|{h'  “M" +kr  +kr} 

<jKo^.,^o  ~  •  ccos(^  +  O2)  ~  Re(^)  •  d  sin(02)  +  Ini(c)  •  ecos(^  +  O2) 

M 

Scattered  Particle  Polarization  Contribution 

cM^oss  -  Re(<i)  •  ecos(^)  +  Im(a)  •  d  sin(^) 

(jM^osk  -  ~  Re(<i)  •  esin(^)  +  Im(a)  •  d  cos(0)  -  Im(/))  •  c 

N 

Recoil  Particle  Polarization  Contribution 

o'No,..,o  ~  “Re(c)  ■ecos{0  +  62)  -  Im(a)  •  csin(^  +  62)-  Im(/))  •  d sin(02) 

cjNo,„k„  =  Re(c)  •  esm{6  +  O2)-  Im(a)  •  c  cos{6  +  ^2)  +  •  d  005(^2) 
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3.  The  AV1 8  Potential 


The  Argonne  National  Labs  (ANL)  Theoretical  Physics  Division  has  been 
constructing  nucleon-nucleon  models  for  over  three  decades.  ANL’s  latest  eighteen 
operator  model,  known  as  the  AVI 8,  has  been  able  to  accurately  describe  the  binding 
energy  of  light  nuclei  through  the  atomic  number  A=  12  and  faithfully  reproduce  the 
experimental  phase  shifts  associated  with  nucleon-nucleon  scattering.  Embedded  within 
the  AVI 8  are  actually  three  phenomenological  models,  one  for  each  possible  two-body 
interaction,  pp,  pn,  and  nn.  Each  of  these  models  contains  an  electromagnetic  ( ), 

one-pion  exchange  ( )  and  intermediate  short-range  ( )  component  [4] 

viPP)  =  Vem  (pp)  +  (pp)  +  (pp) 

v{np)  =  {pn)  +  {pn)  +  (;?«)  (3.1) 

v{nn)  =  Vem  {nn)  +  {nn)  +  {nn) 

The  sum  of  the  later  two  components  of  each  model,  andv^  ,  constitutes  the  strong 
force  contribution.  The  proton-proton  model,  proton-neutron  model,  and  neutron-neutron 
model  are  summarized  in  Section  3.1,  Section  3.2,  and  Section  3.3  respectively. 

3.1.  The  AV18  Proton-Proton  Model 

The  most  complicated  model  within  the  AV 1 8  is  the  pp  as  it  deals  with  the 
interaction  of  two  charged  particles.  Therefore,  a  majority  of  chapter  3  will  deal  with  the 
development  of  this  model.  This  chapter  has  three  sub-sections,  one  for 
Vem(PP%  v^{pp),  andvE{pp).  . 
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3.1.1.  The  Proton-Proton  EM  Potential,  v^^{pp) 


The  form  of  the  eleetromagnetic  portion  of  the  potential  utilized  by  the  AV 1 8  has 
the  general  strueture, 

^EM  iPP)  =  ^Cl  iPP)  +  ^C2  iPP)  +  ^DF  iPP)  +  ^VP  iPP)  +  ^MM  iPP)  (3  -2) 

where  labels  the  one-photon  exehange,  labels  the  two-photon  exchange,  labels 
the  Darwin-Foldy  term,  Vyp  labels  the  vacuum  polarization,  and  labels  the  magnetic 
moment  contributions. 

The  Coulombic  interaction  is  governed  by  the  exchange  of  virtual  photons,  one 
for  each  charged  particle  involved  in  the  interaction.  Hence,  a  scattering  event  involving 
a  single  proton  would  undergo  an  exchange  of  a  single  virtual  photon,  v^j ,  and  a  two 

proton  event  would  undergo  an  exchange  of  two  virtual  photons,  .  The  one  and  two 
photon  exchange  components  for  pp  scattering  are  represented; 

Vci{PP)  =  (^'— -  (3-3) 

r 


Vc2iPP)  = 


-aa' 


r 


(3.4) 


interaction  between  the  electron  and  the  photon.  The  F  function,  (x)  first  proposed 
by  Auerbach  [43]  reflects  the  finite  size  of  the  nucleon  charge  distribution  and 

compensates  for  the  —  singularity  in  as  r  approaches  zero  [4],  It  has  functional  form, 
r 


F^'’{x)  =  \- 


f 


11 


3  1  ^ 

^  2  ,  ^  3 


1h - X^ - X  H - X 

V  16  16  48 


(3.6) 


withx  =  hr  and  b  referring  to  the  exponential  cutoff  parameter  4.27  fm^ .  Auerbach  F 
functions  appear  in  all  AVI 8  electro-magnetic  components. 

From  special  relativity,  it  is  known  that  as  a  particle’s  velocity  approaches  the 
speed  of  light,  the  mass  of  the  particle  increases.  Since  the  domain  of  the  AVI  8  extends 
up  to  the  pion  production  threshold  of  350  MeV,  (approximately  37  percent  of  a  proton’s 
rest  mass)  relativistic  effects  cannot  be  ignored.  The  Darwin-Foldy  term. 


Vdf{pp)  =  - 


a 


AM 


2 


(3.7) 


is  included  to  compensate  for  the  relativistic  effects  of  a  finite  sized  charged  spin  14 
particle  interacting  with  an  electromagnetic  field.  Here,  Mp  refers  to  the  mass  of  the 
proton  mMeV  and  Fg{x)  refers  to  another  one  of  Auerbach’s  F  functions. 


FAx)  =  b^ 


f  1  1 


1 


4 


- 1 - x-\ - X 

vl6  16  48  j 


(3.8). 


The  vacuum  polarization  term  ( v^p  ), 


^(PP)  = 


2aa '  FA  (r)  ‘ 
2n  r 


00 

f  g(-2»isra) 

1+  ^2 

J 

1 

2x"_ 

X 


-dx , 


(3.9) 
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is  a  quantum  correction  to  classical  electrodynamics.  As  shown  in  Figure  3.1,  virtual 
charged  particle  pairs  form  around  the  positively  charged  proton,  which  reduces  the 
protons  effective  charge. 


Figure  3.1.  Vacuum  Polarization  Shielding  Effect.  Virtual  charged  particle  pairs  form  around  a 
proton  in  a  vacuum.  The  formation  of  the  virtual  charged  particle  pairs  reduces  the  protons  effective 
charge. 

The  effective  charge  increases  as  distance  to  the  charged  particle  decreases  [59].  The 
vacuum  polarization  term  is  required  to  fit  low  energy  scattering  data  [4].  Contributions 
homv yp{pp)  cease  to  be  important  above  about  30  MeV [27].  In  the  above  equation  for 

the  vacuum  polarization  contribution,  is  the  mass  of  the  electron  in  MeV  and  F^'’{r)is 

the  same  F  function  used  for  v^j  back  in  the  one  photon  exchange  term. 

The  last  term  of  equation  (3.2)  is  the  particles  magnetic  moment,  Vf^^{pp) .  For 

the  proton-proton  interaction,  the  magnetic  moment  contribution  takes  the  form, 
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.X  a  2 

^mm(PP)  = -^772 


4M 


3  r 


a  ..  ,xF„(r) 


2M: 


.(4u-1)^^L-S{3.10) 


where  ju  labels  the  magnetie  moment,  (Tj  •  (Tj  labels  the  dot  produet  of  the  Pauli  spin 


matriees  of  the  two  protons,  and  L  ■  S  labels  the  dot  produet  of  the  total  orbital  angular 
momentum  and  total  spin  angular  momentum.  In  addition  to  F^(x)  already  introduced. 


this  expression  also  contains  two  more  of  Auerbach’s  F  functions. 


F^(x)  =  l- 


1  2  1  3  1  4 

+  -X  + — X 

6  24 

f  17  1^ 

1  ,  1^2,'  3  ,  r  4 

1  +  X  +  — X  H - X  H - X 


1  -*-2-*-  3  -*-4  -*-7 

1  +  XH - X  H - X  H - X  - X 

V  2  6  24  144 


exp(-x) 


(3.11) 


48 


48 


exp(-x) 


and  the  tensor  force,  Sn- 

=  3(crj  •  f )(cr2  •  f ) - CTj  •  CTj  (3.12) 

The  tensor  force  term,  Sn,  stems  from  the  slight  asymmetry  associated  with  the 
observed  nature  of  the  deuteron  wave  function.  If  the  deuteron  were  comprised  of  merely 
two  point  particles,  then  the  potential  would  be  entirely  symmetric  (pure  S  state)  and  be 
completely  described  by  spherical  harmonics.  This  is  not  however,  what  is  observed. 

The  interacting  bodies  in  a  nucleon  collision  are  not  actually  the  nucleons  but  the 
constituents  of  the  nucleons,  the  quarks.  Protons  contain  two  up  and  one  down  quark 
whereas  neutrons  contain  two  down  and  one  up  quark.  It  is  the  interaction  between  the 
quarks  that  gives  rise  to  the  slight  4  percent  D  state  component  of  the  deuteron  wave 
function.  [53].  A  list  of  all  AVI  8  constants  is  provided  in  Table  3.1. 


24 


Table  3.1.  Constants  nsed  by  tbe  AV18  code. 


AV  18  Constants 

he  (MeVfm) 

197.32705 

(MeV/c^) 

134.9739 

(MeV/c^) 

139.5675 

Mp  (MeV/c") 

938.27231 

M„  (MeV/c") 

939.56563 

137.03599 

2.79285 

!^n 

-1.91304 

/ 

0.075 

b  (fm‘) 

4.27 

r^(fm) 

0.5 

a  (fm) 

0.2 

2.1 
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3.1.2.  The  Proton-Proton  One  Pion  Exchange  Potential,  v^{pp) 


In  quantum  electrodynamics,  the  force  between  two  particles  can  be  described  as 
being  mediated  by  the  exchange  of  virtual  photons.  Photons  have  no  rest  mass  so  the 
range  of  the  electromagnetic  interaction  is  infinite.  It  was  Yukawa  [26]  that  first 
suggested  that  the  concept  of  virtual  exchange  could  be  extended  to  describe  the  strong 
force  between  two  nucleons.  In  his  theory,  he  proposed  that  since  the  strong  force  has  a 
finite  range,  it  could  be  described  as  being  mediated  by  the  exchange  of  virtual  particles 
of  non-zero  mass  called  mesons  [9,  60].  This  simple  extension  of  quantum 
electrodynamics  is  depicted  as  a  Feynman  diagram  in  Figure  3.2. 


Figure  3.2.  Feynman  Diagram  of  a  Pi  Meson  Exchange  between  Two  Arbitrary  Nucleons  [9] 


26 


Here,  a  single  meson  n  mediates  the  interaetion  between  two  arbitrary  nueleons 
labeled  by  Sinee  the  virtual  intermediate  state  shown  here  is  at  an  energy  higher 
than  the  initial  or  final  states,  the  Heisenberg  uneertainty  prineiple, 

(3.13) 

me 

plaees  a  limit  on  how  long  the  intermediate  meson  state  may  exist.  This  length  of  time 
depends  inversely  on  the  mass  of  the  particle.  Since  a  particle  cannot  travel  farther  than  c 
times  this  lifetime.  At ,  a  limit  is  established  on  the  range  of  the  force  component  the 
particle  can  carry.  Some  calculated  meson  lifetimes  ( At )  and  force  component  ranges 
(Ax)  based  on  the  Standard  Model  of  Fundamental  Particles  and  Interactions  are  listed  in 
Table  3.2. 


Table  3.2  Approximate  Meson  Masses,  Lifetimes  ( At )  and  Ranges  ( Ax )  Calculated  from  the 
Standard  Model  of  Fundamental  Particles 


Name 

Mass  (Me  V) 

At  (s) 

Ax  (fm) 

Pion 

140 

4.70*10-"‘’ 

1.41 

Kaon 

494 

1.33*10-^'' 

0.40 

Rho 

770 

8.54*10'^^ 

0.26 

D+ 

1869 

3.51*10-^’ 

0.11 

Eta-c 

2979 

2.2*10'^'’ 

0.07 

As  shown  in  the  table,  the  lightest  meson,  the  pion,  has  a  range  roughly  equivalent 
to  the  observed  1-2  femtometer  range  of  the  strong  force.  This  is  a  reason  why  the  strong 
force  can  be  viewed  as  an  exchange  of  virtual  pions.  In  the  laboratory,  positive 
( 71^ ),negative  ( 7r_ ),  and  neutral  {tIq)  pions  have  been  observed.  The  negative  pion 
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is  the  anti-particle  for  the  positive  pion  and  both  have  an  approximate  rest  mass  of 
140  MeV.  On  the  other  hand,  the  neutral  pion  {tIq)  has  a  slightly  lower  rest  mass  of  135 
MeV  and  is  its  own  anti-particle  [60]. 

The  charge  independent  Hamiltonian  based  on  the  exchange  of  virtual  pions 
between  two  nucleons  can  be  expressed  as  [9], 

i=\  ' 

where  is  a  three  component  vector  in  isospin  space  which  represents  the  three  different 


pions, 

fl3 

|;r,)  =  |/  =  l,4=+l}=  0 

loj 

"0^ 

|;z- )  =  =  =0)=  1  (3.15) 

lo. 

0 
0 
1 

r.  labels  the  location  of  the  nucleon,/  labels  an  arbitrary  coupling  constant,  cr.  labels 
the  Pauli  spin  matrices,  t-  labels  the  isospin  matrices  and  ju  labels  the  inverse  scattering 
length/cutoff  parameter. 


me 


(3.16) 


The  isospin  quantum  number,  r ,  enables  us  to  treat  the  proton  and  neutron  as  a 
different  state  of  the  same  particle,  the  nucleon  [7].  Possible  values  for  the  isospin 
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projection,  ,  in  the  uncoupled  basis  are  -14  and  14,  with  the  neutron  designated  as  -  14 
and  the  proton  as  +  14.  In  the  coupled  basis,  the  total  isospin  T  reflects  the  total  isospin  of 
the  interaction  pair  analogous  to  the  total  spin  quantum,  S,  either  0  or  1 .  In  the  triplet 
state,  T=  1 ,  three  possible  combinations  exist, 

n)=|i) 

^|T^)  +  |^t)  =  |0)  (3.17) 

V2 

|U)  =  |-i) 

Here,  a  projection  of +1  indicates  a  proton  pair,  -la  neutron  pair,  and  0  a  linear 
combination  of  proton  and  a  neutron.  The  singlet  state  is  also  a  linear  combination  of 
protons  and  neutrons, 

=  (3.18) 

Since  we  are  dealing  with  relativistic  energies  and  spin  0  pions,  the  relativistic 
Klein-Gordon  equation, 

=  source  terms  (3.19) 

is  used  to  construct  the  wave  function.  Taking  the  static  limit,  removes  the  time- 
dependence  from  the  d’Alembertian  such  that, 

(3.20) 

and  we  are  left  with, 

=  (3.21) 

i=\ 

From  equation(3.21),  the  pion  field  generated  by  nucleon  2  is  determined  to  be. 
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An/j.  \r 


(3.22) 


An  assessment  of  the  interaction  energy  of  nucleon  1  with  this  pion  field  is  obtained  via 
substitution  into  the  interaction  Hamiltonian  from  equation  (3.14), 


(3.23) 


Eisenberg  [9]  has  reduced  this  equation  by  performing  the  specified  derivatives  to  find, 


\  J  f  ^  A  e  1  3  3 

v^=—mc  -  r,»r,  <cr,»cr, +0,-  IH - 1 - 

"  3  "  [Anhc]^'  2  ur  (nr 


(^r)  J  t^r 


(3.24) 


(3.25) 


where  T(r)  and  Y(r)  label  the  tensor  and  Yukawa  terms  in  equation  (3.24).  As  indicated 
by  the  isospin  coupling  term,  (fj'fj ) ,  a  slightly  different  solution  exists  for  each  possible 
interaction  combination  (nn,  pn,  or  pp). 


The  inclusion  of  a  pion  mass  ratio  — —  and  an  exponential  cutoff  function, 

m. 


1  -exp(-c^r^)  are  the  only  differences  between  the  Eisenburg  solution  above  and  the 
documented  form  of  the  AVI  8  One  Pion  Exchange  Potential  [4], 


\  fn  n  2~ 

=  <^i*<^2^(0(l“ex/2(-c^r^)l  +  5'j2T(r)(l-eyfi(-c^r^)l  (3.26) 

3  4n:hc  [  in  ^  - 


The  scaling  mass,  ms  non-dimensionalizes  the  coupling  constant,  - and  has  an 

Anhc 


experimentally  determined  value  of  0.075  [4].  The  exponential  cutoff  function. 
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\-exp{-cy) 


(3.27) 


ensures  both  the  Yukawa  ( Y^{r)  )  and  Tensor  ( T^{r)  )  eomponents  of  the  one  pion 

exehange  potential  terminate  smoothly  as  r  ^  0  [4],  In  the  AVI  8,  the  eutoff  parameter 
c^is  set  at 2.1  fm^  [4]. 

3.1.3.  The  Proton-Proton  Short  Range  Potential  v^{pp) 

Data  obtained  from  low  energy  seattering  events  not  only  provides  insight  into  the 
nuelear  potential’s  physieal  eharaeteristics  but  also  aids  in  the  development  of  a  realistie 
parameterizations  ofv^  .  Observations  have  shown  that  any  eonstructed  short-range 
nuelear  potential  must:  deeay  rapidly  around  1-2  fm,  be  strongly  repulsive  at  distanees 
less  than  0.5  fm,  be  strongly  attraetive  between  0.5  fm  and  2  fm,  have  a  depth  on  the  order 
of  40  MeV,  and  depend  on  the  eaeh  nucleon’s  intrinsic  spin/isospin.  With  these  ideas  in 
mind,  we  begin  our  discussion  of  the  theory  behind  phenomenological  potentials  [61]. 

Since  any  nucleon  can  be  labeled  by  its  position  r  ,  momentum  p  ,  spin  a ,  and 
isospin  f ,  any  realistic  parameterization  of  the  two  nucleon  nuclear  potential  must  begin 
with  the  general  functional  form, 

v{r^,r^,p^,p^,B^,B^,f,Tf  (3.28) 

Eisenbud/Wigner  and  later  Okubo/Marshak  proposed  that  this  general  form  could  be 
restricted  since  any  potential  should  adhere  to  several  invariance  requirements  [28,  61]: 

1 .  Translational  Invariance  -  The  concept  of  translational  invariance  simply 
means  that  the  potential  should  depend  only  on  the  relative  separation 
between  the  two  nucleons,  r  =  r^-r^,  not  where  the  two  nucleons  are  in  space. 
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The  system  is  elosed  and  no  external  forces  act  on  the  system  other  than  those 

dp 

between  the  two  particles  so  total  linear  momentum  is  conserved,  F  =  —  =  0  . 

dt 

2.  Rotational  Invariance  -  The  potential  should  be  constructed  such  that  it  is 
invariant  under  proper  rotation.  Since  the  system  is  again  closed,  there  are  no 

external  torques  acting  in  the  system,  r  =  —  =  0  and  the  total  angular 

dt 

momentum  L  =  rx pis  conserved. 

3.  Galilean  Invariance  -  The  concept  of  Galilean  invariance  implies  that  the 
potential  should  depend  only  on  the  relative  momentum  between  the  two 
nucleons,  p  =  p^-  P2  and  be  independent  of  inertial  reference  frames.  A 

change  of  coordinates  involving  a  constant  velocity  does  should  not  influence 
the  potential.  The  same  laws  of  physics  apply  to  all  inertial  reference  frames. 

4.  The  potential  should  be  Hermitian.  This  means  that  any  potential  matrix 
should  be  diagonalizable  by  a  unitary  transformations  and  that  all  the 
eigenvalues  along  the  diagonal  should  be  real  (observable). 

For  the  moment,  we  limit  the  discussion  to  r  ,p  ,5^,  and  ■ 

Given  the  above  invariance  requirements,  the  functional  form  of  the  nucleon- 
nucleon  potential  can  be  restricted  to  [9], 

v{r,p,B^,B^),  (3.29) 

In  addition  to  the  four  vectors  r  ,p  ,B^,  and  B^ ,  six  other  vectors  may  be  constructed, 

(TjXCTj,  B^xr ,  B^xp,  (TjXr,  B2XP,  rxp.  (3.30) 
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Note  that  we  do  not  need  to  eonsider  other  eombinations  of  vectors  or  terms  containing 
more  than  one  cTj  and/or  cr^  since  the  identities, 


Ax[bxC^  =  b(a-C^-C(a-B^ 
and  the  Pauli  matrices  property, 

=  A^»B  +  ia»{^xB^ 


(3.31) 


(3.32) 


can  be  used  to  reduce  them  back  to  a  lower  form. 

From  these  six  vectors  and  r  ,p  ,d^,  and  ,  the  following  scalar  combinations 

are  possible  for  the  nucleon-nucleon  interaction  [9]; 

r-r=r^,  ^I'P’  ^ '  P-^ 

r-(o-iXo-2),  ^-(o-ixo-j),  a^-{rxp),  a2-{rxp), 

((Ti-(T2)(r  •^),  (ai-r)((T2-r),  ((T;  •r)((T2  •^),  {a^-p){a2-r) 

{apr){r-p),  {a2-r){r-p),  {r  ■  p)r -{a^xa^),  (r  •r)((T2 -r) 

As  we  are  merely  trying  to  gamer  an  understanding  of  how  nucleon-nucleon  potentials 
are  developed,  we  have  limited  this  discussion  to  1®*  order  terms  in  momentum.  In 
practice,  higher  order  terms  in  momentum  are  often  included  in  potentials  like  the  AVI  8. 
To  these  first  four  requirements,  four  more  are  typically  added  [9]; 

5.  The  potential  must  commute  with(fj  +  ) »  which  implies  charge  must  be 

conserved.  This  requirement  allows  for  five  possible  operators  in  isospin 
space, 

1,  TpT^,  h+T^,  h-T^,  T,XT2  (3.34) 
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In  addition  to  charge  conservation,  charge  independence  is  also  typically 
imposed  so  that  only  the  scalars,  1  and  Tj  •  Tj  ,  which  are  invariant  under 

rotation  in  isospin  space  may  appear  in  the  nuclear  potential  [9]. 

6.  Parity  Invariance  -  The  potential  should  be  independent  of  reflection  so  only 
even  powers  of  position  and  momentum  may  appear  together  in  the  potential 
expansion  [62], 

v(r , o-j ,  0-2 ,  Tj ,  Tj )  =  v(-r,  (3.35) 

Of  the  scalars  mentioned  in  (3.33)  only 

1,  <Ti-(T2,  r-p, 

B,-[rxp),  ap[rxp), 

-p),  ((Ti-r)((T2-r),  (3.36) 

{a,-r){a^-p),  {a^-p)[a^-r), 

[r  ■p){a^-r){a^-r) 

are  invariant  under  parity  inversion  in  addition  to  the  isospin  scalars  1  and 

TpT^. 

7.  Time  Reversal  Invariance  -  The  potential  should  be  invariant  under  time 
reversal  operations  which  implies  that  spins  and  momentum  can  only  couple 
only  in  even  number  combinations  [62], 

v{r,p,(j^,<j^,T^,T^)  =  v{r,-p,-<j^,-<j^,T^,T^)  (3.37) 

Time  reversal  invariance  then  further  restricts  the  allowable  scalars  from 
(3.36)  to 

1,  <Ti-(T2,  (<?!  •r)((T2 -r 

Bp[rxp),  B^-[rxp), 


(3.38) 
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Since  the  interaetion  always  involves  both  partieles,  the  last  two  terms  of 
(3.38)  ean  be  expressed  as  the  linear  eombination^((Ti  +o'2){rxp^.  Typieally 


this  linear  eombination  is  written  as  Z  •  S'  given  the  definitions  for  total  orbital 
angular  momentum,  L  ,  and  total  spin,  S  , 


L  =  rxp 


(3.39) 


S  =  d'^+ a  2 

Protons  and  neutrons  do  not  change  their  identity  under  time  reversal  so  the 
isospin  sealars  1  and  z)  •  Tj  remain. 

Gammel  and  Thaler  in  1957  utilized  eombinations  of  these  sealars  to  ereate  the 
first  nueleon-nueleon  phenomenological  model  [28], 

^Gammel  ^central  ir)  +  v,^,^Xr)Sn+v,s{r)L-S .  (3.40) 


by  expressing  v^^„,^„,(r),  Vy,^„^^^(r)  and  (r)  in  the  Yukawa  form  v(r)  =  —  and  by 

r 

expressing  the  tensor  foree  Sjj  (Seetion  3.1.1)  as  a  eombination  of  the  sealars 
Sj2  =  3(crj  •r){(j2  •r)  -  Cj  ‘CTj  .  Although  the  Gammel  and  Thaler  model  was  unable  to 

reproduee  all  observed  nn  and  pp  seattering  data,  it  did  serve  as  a  basis  for  further  short- 
range  model  development. 

The  AVI  8  short-range  expression  for  given  by 


ipp)  =  ir)  +  vff  {r)L^  +  {r)S,2  +  vf.’J  {r)L  •  5  +  {r){L  -S  f  (3.41) 


’  The  linear  combination  compliment  (cTj  ~  <72  )  ^  P^  vanishes,  as  it  is  odd  under  interchange  of 


particle  labels  [21]. 
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includes  second  order  terms  in  momentum  [L-S^  and  .  Each  of  the  five  functions  in 
equation  (3.41)  can  be  expressed  as, 

+  +{^rf  R^’^]W{r)  (3.42) 

Here,  the  i  subscript  labels  the  central,  L^,  Tensor,  L»S ,  and  {L»Sy  potentials.  The 
superscripts  S/  T  label  the  total  spin/isospin  of  the  two-nucleon  system  and  /u  is  the 
average  of  the  charged  and  neutral  pion  masses.  The  constants  /,  P,  Q,  and  R  are 
obtained  from  curve  fitting  observed  scattering  data  (Table  3.3). 

Since  T=1  for  the pp  interaction,  only  two  possible  potentials  exist  ioxVj^{pp) , 

one  for  the  singlet  state  S=0,  v”  ' (/»/») 


v''r{PP)  =  T 


,04 

central 


(r)  +  v°’*(r)E 


(3.43) 


and  another  for  triplet  state  S=l,  v^{pp)  . 


Vr\pp)  =  v^ 


14 

central 


i^r)W;l{r)L^+vy_Xr)S, 


.  14 

12  H  s 


{r)L-S  +  vl^^^{r){L-Sr  (3.44) 


Equation  (3.42)  also  includes  a  short-range  repulsive  core  represented  by  the  Wood- 
Saxon  function,  [4] 


W{r)  = 


1 


\  +  Exp 


^r-r^ 
[  a  ) 


(3.45) 


and. 


rr 


T;{r)  = 


1  3  3 

1  +  —  +  - 


vv 


Pr  [pr) 


/  ,x2 


pr 


(3.46) 
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which  has  Yukawa  form, - .  The  parameters  ro  (core  radius),  a  (surfaee  thickness), 

jur 


and  c,^(eutoff)  are  given  in  Table  3.3[4]. 


Table  3.3.  {pp)  Strong  Component  Constant  Parameters  of  the  AV18  model 


Channel 

TYPE 

IMeV 

VMeV 

QMeV 

KMeV 

S=0,  T=1 

Central 

-11.27028 

3346.6878 

1859.5627* 

0 

0.12472 

16.780 

9.0972* 

0 

Central 

-7.62701 

1815.4920 

969.3863* 

1847.8059 

u 

0.06709 

342.0669 

185.4713* 

-615.2339 

Tensor 

1.07985 

0 

-190.0949 

-811.2040 

LS 

-0.62697 

-570.5571 

-309.3605* 

819.1222 

{LSf 

0.74129 

9.3418 

5.0652* 

-376.4384 

3.2.  The  AV18  Neutron-Proton  Model 

Since  the  AVI 8  neutron-proton  model  is  not  as  complicated  as  the  AVI 8  proton- 
proton  model,  this  seetion  will  be  abbreviated.  This  section  also  has  three  sub-seetions, 
one  for  eaeh  component. 

3.2.1.  The  Neutron-Proton  EM  Potential,  Vj,^{pn) 

The  form  of  the  eleetromagnetie  portion  of  the  AV 1 8  model  has  the  general 

structure. 


(pn) 


(3.47) 
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where  again  v^j  labels  a  one-photon  exchange  and  labels  a  magnetic  moment 

contribution.  The  two  photon  exchange,  Darwin-Foldy,  and  vacuum  polarization 
components  are  omitted  since  neutrons  interact  minimally  with  the  electromagnetic  field. 
The  remaining  two  terms  differ  slightly  from  the  proton-proton  counterparts 
Vc^{pp)Mi&VMM{PP), 


Vc,{pn)  =  al3^^ 


(3.48) 


a 


-PnPn 


AMMp  '  n  3 


2  „  ,  _  .  F,{r) 


Fg{r)cr,-cr^  + 


3  ^12 


a 


2MM. 


-Pn 


F,{r) 


f 


L-S  +  L 


v2 


(3.49) 


In  the  above  set  of  equations,  [3^  is  the  form  factor  associated  with  the  electron-neutron 
interaction  at  the  point  of  zero  momentum  [30,  63].  Since  fixed  neutron  targets  do  not 
exist  in  nature,  (3^  at  zero  momentum  is  extrapolated  from  high-energy  (GeV  range) 

electron-deuteron  scattering  experiments  [96].  In  the  AV18,  =  0.0189ym  [14].  The 

only  other  term  not  already  introduced  in  Section  3.1.1.  is  another  one  of  Auerbach  F 
functions. 


F^’'{x)  =  ^^(15x  +  15x  -l-6x  +x  )exp(-x) 


(3.50) 


Here,  again  x  =  br  and  b  refers  to  the  exponential  cutoff  parameter  of  4.27  fm\  is  the 
nucleon-proton  reduced  mass,  and  —  (ctj  -  (Tj  )  represents  a  mixing  of  spin  singlet  and 

v2  J 

triplet  states.  This  mixing  contributes  minimally  to  the  magnetic  neutron-proton 
scattering  amplitude  [4]. 
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3.2.2.  The  Neutron-Proton  OPEP,  v^ipn) 

Since  the  theory  associated  with  v^{pn)  has  already  been  presented  in  Seetion 
3.1.2,  this  discussion  begins  with  the  general  expression  forv^  presented  in  equation 
(3.25), 

1  (  ^ 

— —  (TpT^){<jpcT^Y{r)  +  S,J'{r)Y{r)]  (3.51) 

*5  y  ^TTil C  J 

Here,  a  eomplex  mixing  of  neutral  pion  ^  and  eharged  pions  ^  ean  transform  the 

ineident  neutron  into  a  proton  and  the  ineident  proton  into  a  neutron.  In  the  AVI 8, 

if  ^ ^ t  fV 
\  ^  J 

(-1)"!  ^\^(l-exp(-cy-))(a,-a,Y(r)  +  S„T(f)r(f))  (3.52) 

V  ^  J 
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(3.55) 


and  two  for  S=l. 


0,1  /-  0,1 
0,0  /  X  0,1 


(r)  +  vJ'(r)Z^ 
ir)  +  vlf{r)L^ 


v^nipn)  =  vlLAr)  +  v^^ir)L^ 
^'fipn)  =  vlLAr)W;{r)L^ 


+  VrLsor  ir)S,,  +  vl:',  ir)L  ■  5  +  v;;  (r)(Z  •  5)^ 

+  v^L.  (^)‘^i2  +  ir)L  •  5  +  (r)iL  •  5)^ 


(3.56) 


Table  3.4.  {pn)  Strong  Component  Constant  Parameters  of  the  AV18  Model 


Channel 

TYPE 

IMeV 

PMeV 

QMeV 

RMeV 

S=0,  T=I 

Central 

-10.66788 

3126.5542 

1746.4298* 

0 

0.12472 

16.780 

9.0972* 

0 

S=0,  T=1 

Central 

-2.09971 

1204.4301 

511.9380* 

0 

-0.31452 

217.4559 

117.9063* 

0 

S=l,  T=1 

Central 

-7.62701 

1815.5315 

966.2483* 

1847.8059 

t 

0.06709 

342.0669 

185.4713* 

-615.2339 

Tensor 

1.07985 

0 

-190.0949 

-811.2040 

LS 

-0.62697 

-570.5571 

-309.3605* 

819.1222 

{LSf 

0.74129 

9.3418 

5.0652* 

-376.4384 

S=I,  T=0 

Central 

-8.62770 

2605.2682 

1459.6345* 

441.9733 

t 

-0.13201 

253.4350 

137.4144* 

-1.0076 

Tensor 

1.485601 

0 

-1126.8359 

370.1324 

LS 

0.10180 

86.0658 

46.6655* 

-356.5175 

{LSf 

0.07357 

-217.5791 

-117.9731* 

18.3935 

3.3.  The  AV18  Neutron-Neutron  Model 


Since  the  neutron  has  an  internal  strueture,  it  is  still  subject  to  electromagnetic,  one  pion 
exchange,  and  short-range  effects.  This  section  summarizes  these  components. 
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3.3.1.  The  Neutron-Neutron  EM  Potential, 


The  electromagnetic  portion  of  the  nn  potential  only  has  one  component, 

^EMi^n)  =  v^^{nn).  (3.57) 

which  results  from  the  magnetic  moment  of  the  neutron.  This  single  component  of 
{nn)  may  be  represented  by, 


a  2 


AM 


2 

3 


T;(r)c7,.c72  + 


(3.58) 


As  previously  mentioned,  the  observed  magnetic  moment  of  the  neutron  was  one  of  the 
first  indications  of  an  underlying  particle  structure. 


3.3.2.  The  Neutron-Neutron  OPEP,  v^{nn) 


For  the  nn  OPEP,  the  AVI 8  assumes  this  solution  to  be  the  same  as  v^{pp)  as 
experimental  data  has  shown  little  variation  in  the  coupling  constants  below  350  MeV. 


1 


.2  r 


v(nn)  =  —m„c 

3  “  Anne 


J 


[l-exp{-cy)) 


[apa,Yir)  +  S,,Tir)Yir)]  (3.59) 


3.3.3.  The  Neutron-Neutron  Short  Range  Potential,  v^{nn) 


Analogous  to  (pp) ,  T=1  for  the  nn  interaction.  Therefore,  two  potentials  are 
possible,  one  for  the  S=0, 

v'/(np)  =  vZ,,Jr)  +  vl^ir)L^  (3.60) 

and  one  for  S=7. 


Vr  (np)  =  V, 


{r)W;Ur)L 


central 


.1,1 


.1,1 


+  VT:nsorir)S,^+V^,S 


{r)L-SW’l^^,{r){L-Sy 


(3.61) 
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Here  as  before,  eaeh  of  the  terms  of  equation  (3.61)  can  be  described  by  equation  (3.42) 
with  the  constants  parameters  I,  P,  Q,  and  R  provided  in  Table  3.5. 


Table  3.5.  v^(nn)  Strong  Component  Constant  Parameters  of  the  AV18  model 


Channel 

TYPE 

IMeV 

VMeV 

QMeV 

KMeV 

S=0,  T=1 

Central 

-11.27028 

3342.7664 

1857.4367* 

0 

0.12472 

16.780 

9.0972* 

0 

Central 

-7.62701 

1811.5710 

967.2603* 

1847.8059 

u 

0.06709 

342.0669 

185.4713* 

-615.2339 

Tensor 

1.07985 

0 

-190.0949 

-811.2040 

LS 

-0.62697 

-570.5571 

-309.3605* 

819.1222 

{LSf 

0.74129 

9.3418 

5.0652* 

-376.4384 

3.4.  The  Numerov  Method 


ANL  typically  utilizes  the  Numerov  method  to  compute  AV 1 8  scattering  phase 
shifts  [5].  This  time  independent  approach  is  ideally  suited  to  linear  second  order 
differential  equations  like  the  Schrodinger  equation, 

^  =k^[E-V{r))xi/{r) 

(3.62) 


dr 


k-  ^ 


which  have  no  first  derivative.  Here,  F is  the  AVI 8  potential,  //  is  the  reduced  mass  of 
the  nucleon  pair,  and  h  is  Planck’s  constant  divided  by  Itt 

To  use  this  technique,  we  first  perform  a  Taylor  series  expansion  of  y/  at  the 
point  (r  +  h) 
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(3.63) 


h'^ 

y/{r  +  h^  =  \i/{r)  +  hy/'  {r)  + —y/”  {r)  + —y/"  {r)  +  {r)  + . 

then  perform  a  separate  Taylor  series  expansion  of  y/  at  the  point  (r-h) , 


y/[r-h)  =  y/{r)  -  hy/'{r)  +  ^y/"{r)-^y/"'{r)  +  ^y/""{r)  + .... 

2  6  24 


(3.64) 


The  sum  of  these  two  expansions  eliminates  odd  powers  of  h  and  leaves  an  expression 
for  ^"(r)  which  includes  an  undesirable  fourth  order  derivative, 


(3.65) 


Here,  h  refers  to  the  integration  step  size  and  should  not  be  confused  with  the 
h  contained  within  in  equation  (3.62).  In  order  to  reduce  this  equation  into  something 


(h^  1 

more  tractable,  we  operate  on  the  Schrodinger  equation  (3.62)  withl  + - j  to 

i  12  , 


generate  terms  analogous  to  the  Taylor  Series  Expansion, 


¥"{r)  +  ^y/""{r)  +  k^y/{r)  +  VC'')]  =  0 


(3.66) 


and  then  back  substitute  into  equation  (3.65)  for  y/"{r)  +  —y/""{r) , 


r  H  pz 

y/{r  +  h)  +  y/{r  -h)-  2y/{r)  +  h^k^y/{r)  -\ - Y\k^y/{r)\  +  0(/z'’)  =  0  (3.67) 

12  dr 
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^2 

All  that  remains  in  equation  (3.67)  is  one  lone  seeond  derivative, 

whieh  is  normally  redueed  to  first  order  through  an  elementary  differenee  formula 
approximation’, 


dr 


(3.68) 


After  some  algebra,  we  arrive  at  the  general  form  for  the  Numerov  method. 


y/{r  +  h)  =  - 


\-  —  ¥k 
12 


i//(r)- 


2,2 


l  +  —  h^k 
12 


y/{r-h) 


l  +  —  h^k^ 
12 


■O(h^) 


(3.69) 


The  Numerov  method  above  is  a  three  point  differenee  formula  that  eomputes 

y/{r  +  h)  based  on  the  value  of  y/  at  two  previous  points  y/{r)  and  y/{r  -h)  .  Of  course  to 

start  the  iterative  method,  we  must  first  seed  a  couple  of  values.  The  first  point  at 

y/{r-h)  is  given  by  the  boundary  condition  at  r  =  0 ,  y/{Qi)  =  0 .  The  second  point  at  a 

step  h  from  the  origin  is  considered  to  be  sufficiently  small  so  y/{r)  is  only  slightly 

greater  than  zero.  From  these  two  points,  ^(r  +  h)  is  determined.  Then, 

y/{r)  ^y/{r  -  h)  and  if/(r  +  h)^  y/(r)  and  a  new  y/{r  +  h)  can  be  computed.  The  process 

is  repeated  until  y/{r)  is  completely  determined  for  a  single  energy  value  [64,65]. 


^  Although  this  difference  approximation  has  an  error  of  0\Jl^  j ,  it  has  minimal  impact  on  the  solution 
since  it  is  multiplied  by  the  much  smaller  . 
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Thus  far,  we  have  shown  how  to  obtain  a  solution  to  the  Sehrodinger  equation  via 
the  Numerov  method  but  not  shown  how  a  phase  shift  can  be  obtained  from  this  solution. 
To  help  motivate  the  phase  shift  discussion,  consider  Figure  3.3. 


Figure  3.3.  Centripetal  Potential  and  Effective  Potential  Comparison.  At  some  point  r^ax 

the  centripetal  barrier  begins  to  dominate  at  arbitrary  points  ri  and  >  r^ax  the  solutions  the  radial 
equation  must  be  a  linear  combination  of  free  solutions  [64] 

Here,  a  centripetal  barrier  potential  and  an  effective  potential  including  both  the 
centripetal  potential  and  one  of  the  AVI  8  potentials  are  presented.  Notice  at  a  point 
roughly  equal  to  the  centripetal  potential  begins  to  dominate.  As  r  continues  to 
increase  past  ,  the  solution  to  the  radial  equation  at  a  point  must  be  a  linear 

combination  of  free  solutions. 
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(3.70) 


¥^l ir)  =  Akr,  {kr^ )  cos 5^  - rt^  {kr^ ) sin 5^  ] 

Here,  A  is  an  arbitrary  eonstant,  is  the  phase  shift,  and  and  the  regular  and 
irregular  spherieal  Bessel  funetions  of  order  L.  We  eannot  ealeulate  the  phase  shift  from 
the  solution  at  alone  sinee  the  equation  eontains  two  unknowns.  Caleulation  of  the 
phase  shift  is  faeilitated  by  determining  a  seeond  solution  at  a  point  Tj  >  also  given  by 
the  solution  to  the  free  radial  equation, 

y/li  (r)  =  Akr^  {kr^ )  eos  -  ij^  {kr^ )  sin  S^~\.  (3.71) 

Given  the  two  solutions,  ,  A  ean  be  eliminated  through  some  more 

algebra  to  obtain  an  expression  for  the  phase  shift  at  a  partieular  energy, 

tan S,  =  V\('-)k(h\)-iSkr,)  ^3  ,2) 

Note  that  the  points  and  Tj  are  typieally  ehosen  to  be  mueh  farther  away  from  the 
interaetion  potential  than  shown  in  the  diagram.  These  positions  are  for  illustrative 
purposes  only  [64,  65]. 

In  order  to  get  a  better  understanding  of  how  phase  shifts  are  eurrently  obtained 
from  the  AVI 8,  a  separate  time  independent  eode  was  developed  to  reereate  the 
published  phase  shifts  at  a  2  MeV  energy  resolution  and  eoordinate  step  size,  0.00153  frn. 
A  eomparison  between  a  phase  shift  obtained  from  this  time  independent  eode  and  a 
published  phase  shift  for  the  S^pn  potential  is  presented  in  Figure  3.4. 


*  Plots  of  this  potential  may  be  seen  in  Chapter  5. 
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Figure  3.4.  ANL  and  AFIT  Numerov  *So  Phase  Shift  Comparison.  ANL’s  published  data  is 
represented  hy  circles  and  the  Numerov  calculation  using  the  above  method  is  represented  by  the 
solid  line.  The  AFIT  Numerov  code  exhibited  sensitivity  to  coordinate  step  size  and  required  twice  as 
long  as  the  CPM  technique  to  complete  the  same  calculation. 


Throughout  this  exercise,  it  was  observed  that  the  fidelity  of  the  phase  shift  produced  by 
the  Numerov  method  was  inherently  sensitive  to  coordinate  step  size.  Step  sizes  greater 
than  0.00153  fm  exhibited  far  greater  noise  than  the  present  in  the  plot  above.  Plots  of 
the  same  '  phase  shift  determined  via  the  channel  packet  method  (chapter  6)  had  a 
higher  energy  resolution  (2  orders  of  magnitude  lower)  and  were  produced  in  half  the 
time  of  the  Numerov  method.  These  CPM  characteristics  are  well  suited  for  examining 
phase  shifts  from  many  body  potentials  where  sharp  scattering  resonances  are  likely  to  be 
encountered. 
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4.  Channel  Packet  Method  Theoretical  Background 


In  quantum  scattering  theory,  the  scattering  operator  S  relates  the  incident  reaetant 
asymptotie  state  \y^i„)  to  the  resultant  produet  asymptotie  state 

(4.1) 

This  unitary,  time  independent  operator  eontains  all  information  of  experimental  interest, 
sinee  only  asymptotie  free  motion  is  observed  in  the  laboratory.  Therefore  given  S , 
physieal  insight  into  both  reaetive  and  non-reaetive  seattering  events  ean  be  obtained  for 
verifieation  against  experimental  data  [51]. 

The  seattering  operator,  S  ,  is  defined  as  the  product  of  two  isometric  Moller 
operators  [45,  51,], 

5  (4.2) 


where. 


Cll  =  lim  exp 


iHt 


exp 


-iH^t 

h 


(4.3) 


Notice  that  the  Moller  operator  above  is  essentially  a  product  of  time  evolution  operators; 
one  eorresponding  to  the  asymptotie  Hamiltonian,  and  the  other  the  full 

Hamiltonian,  H .  When  the  Moller  Operator  is  applied  to  some  given  reaetant  state 
\y^i„)  ,  the  wave  paeket  is  propagated  baekward  in  time  under  the  asymptotic 
Hamiltonian,  ,  to  an  intermediate  state  at  t  =  -go  then  forward  in  time  under  the  full 
Hamiltonian  baek  to  its  initial  position  att  =  0 .  The  other  Moller  Operator  Q  has  a 
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similar  effect  on  the  product  state  |  'j ,  first  propagating  the  wave  packet  forward  to 
time  t  =  +CO  then  backwards  in  time  back  to  t  =  0 .  The  action  of  propagating  these  wave 
packets  to  and  from  the  asymptotic  limit  integrates  information  about  how  Hq  is  modified 
by  H  into  each  wave  packet  prior  to  interaction.  These  new  wave  packets  are  designated 
as  the  reactant  and  product  Moller  states,  |  ^  |  y/^^  and  |  ^  |  ^  , 

respectively.  Throughout  the  remainder  of  this  section,  +  subscripts  reference  the 
reactant  wave  packet  and  -  subscripts  reference  to  the  product  wave  packet. 

In  order  to  define  these  Moller  states,  it  is  necessary  to  first  define  a  useful  basis 

set  in  which  to  construct  the  initial  reactant  and  product  states  |  y/^^  ^  and  |  y/^^^  'j  .  A 
convenient  choice  is  the  momentum  basis,  |A:) .  The  usefulness  of  this  choice  becomes 
apparent  when  we  consider  in  the  momentum  basis,  |  k) 


2r  2 


-  fi  k 

Here  we  have  expressed  in  terms  of  the  momentum  eigenvalues -  where  h  is 

2// 


Plank’s  constant,  k  is  the  wave  vector  conjugate  to  r,  and  /u  is  the  reduced  mass.  If  we 
apply  the  Moller  operator  to  these  states 


(4.5) 


and  use  the  intertwining  relation  =  //Q+  [51] 
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(4.6) 


i)2^2  *2r2 

H  n^\k)  =  ^Cl^\k)  =  —\K 


2// 


2// 


we  see  that  not  only  do  we  have  a  mapping  between  the  eigenbasis  \k')  of  Hq  and  the 
eigenbasis  |  A:^)  and  |A:_)  of  H ,  but  we  also  see  that  the  eigenvalues  between  the  basis  are 


the  same. 


H  \k  )  = 


In 

n^k^ 

In 


(4.7) 


k_) 


The  expansion  of  the  reactant  Moller  state  xj/ ^  )  and  product  Moller  y/_  )  state  in  the 


momentum  basis  set  k^'j  have  the  form, 


(4.8) 


x//_)  =  Q_\xnout)  =  ^-\  ri_(k)\k)dk=  ^  ri_{k)\k_)dk 


(4.9) 


where  //(A:^)  and  rj{k_)  represents  the  initial  expansion  coefficients  of  the  product  and 
reactant  wave  packets. 

With  the  two  Moller  states  now  defined,  we  are  now  ready  to  determine  the 
reaction  probability  S  between  these  two  states, 

)  =  {fFout  I  I  fFin )  =  {y^out  I  •  (4.10) 
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To  help  further  develop  a  usable  expression  for  S  in  the  standard  momentum 


representation,  consider  now  the  arbitrary  state  \A^{E)'^  represented  by  the  Fourier 
transform  of  the  time  evolution  of  the  Moller  state, 

|4(^))  =  J^^exp|-^i^|exp  \y/^)  (4.11) 

Since  the  Hamiltonian  can  be  expressed  as  an  energy,  E,  and  both  energies  can  be 

expressed  in  wave  vector  form  E  = -  and  E '  = - ,  we  have 

l/u  IjU 

.  hk'^  _hk^  ^ 

A^{E))  =  ^dt  r+)- 

Here,  the  reactant  Moller  state  is  independent  of  t,  so  the  time  integral  has  the  delta 
function  solution. 

If  we  substitute  in  the  expansion  for  y/_^^  from  equation  (4.8)  and  use  two  delta  function 
identities, 

(4.14) 

=  ^|— j-[^(x  +  «)  +  ^(x-«)]  (4.15) 

we  have. 


AAE))  =  2k  8 
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(4.16) 


\ASE))  =  ^^\dkri,(li,)[s(k: -k)^S(k:  -«:)]|i;) 

Only  two  possible  solutions  for  the  integral  above  exist,  one  for  k'  =  +k  and  one  for 
k'  =  -k .  Therefore, 

1 4  (£»  =  In  ['?.  )  I  )]  ■ 

Here,  it  is  important  to  remember  that  the  +  subseript  refers  to  the  reaetant  state  and  the  - 
subseript  refers  to  the  produet  state  sinee  the  notation  beeomes  more  challenging  as  we 
proceed  further.  The  reaction  probability  then  is  given  by  the  overlap  between  the 
evolving  reactant  state  and  the  stationary  product  Moller  state, 

{yy_  1 4  {E))  =  I  dk'riik')  {k'  \  4  (t))  (4.18) 

{¥-\AiE))  =  ^^^\dk'ri^_{k'){^-\  (+^+)|+^+)  +  ^+ (-^+)|“^+)  (4-19) 

{¥-  1 4  iE))  =  j  [4  {k')ri^  ( +4 )  (4  I  +4  )  +  d  (^')?7+  ( -K )  (^'  I  -K  )]  (4-20) 

Here  there  also  exists  degeneracy  in  A:  analogous  to(4.16)  which  enables  (4.20)  to  be 


written  as  the  sum  of  four  similar  integrals, 

{¥- 1 4  iE))  =  J  dk:rf_  (+4)?;+  (+4 )  (+4 1+4)  + 

j  dk'Tii*_  (+4 )?7+  ( -4 )  (+4 1  -4 )  + 
jjA:'4(-4)/7+(+4)(-4|+4)  + 
jjA:'4(-4)?7+(“4)(-4|-4) 


(4.21) 
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The  integrals  vanish  provided  the  orthogonality  relations  from  Taylor  [51], 


I  1 5 1  ^)  =  (^’  I  I  k)  =  S{li  -  k)  ^ 


(4.22) 


yielding  an  expression  for  the  S  matrix  elements  in  the  momentum  representation 


iTrjLi 

/i|A:| 

2kh 

/i|A:| 

Injii 


V-{+k)ri^{-k)S^,_,  + 
V-{-k)ri^{+k)S_,^,  + 
V-i-k)ri^{-k)S_,^_, 


(4.23) 


Here,  the  sign  of  the  k  subseript  on  S  indieates  the  wave  packets  direction  of  travel,  either 
toward  or  away  from  the  interaction  region.  Therefore,  in  the  case  of  the  symmetric 
square  well,  the  first  two  terms  would  represent  the  probability  of  transmission  and 
reflection  across  the  well  from  an  incident  wave  packet  from  the  left,  whereas  the  last  two 
terms  would  represent  the  probability  of  transmission  and  reflection  from  an  incident 
wave  packet  from  the  right.  Since  only  one  particular  set  of  matrix  elements  can  be 
determined  for  each  momentum  pair,  we  can  express  (4.23)  is  the  compact  form. 


h  ki^E') 


2K/^ri^_{±k{E))ri^[±k{E)) 


{¥-\ME)), 


(4.24) 


This  is  the  CPM  representation  of  the  S-Matrix  [36-38]. 

Solutions  to  equation  (4.24)  hinge  on  the  numerical  evaluation  of  the  inner 
product  I  A^{E)'^ ,  which  is  given  by  the  Fourier  Transform  of  the  correlation  function. 


C{t) ,  between  the  time  evolution  of  the  reactant  and  stationary  product  Moller  states. 
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Scattering  phase  shifts  5 {E)  are  obtained  directly  from  5*+^  {E)  if  we  reeall  that 

solutions  may  be  expressed  in  eomplex  polar  form^  as  exp(z<5')  when  the  scattering  matrix 
is  diagonal  (i.e.  no  coupling  between  states), 


5{E)  = 


(4.26) 


Taylor  mentions  a  few  constraints  to  scattering  calculations  [51].  The  first  states 

_2 

the  potential  should  be  no  more  singular  than  f  at  the  origin,  a  eriterion  easily  met  by  all 
existing  nuelear  models.  The  seeond,  however,  requires  the  potential  deeay  faster  than 
r  at  infinity.  This  is  a  clear  problem  for  the  long-range  Coulombic  tail,  whieh  governs 
charged  nueleon  interaction  at  large  r.  Fortunately,  a  teehnique  for  handling  the  tail 
exists  in  the  literature  and  an  outline  of  the  theory  is  provided  in  Section  5.6 


^  The  scattering  matrix  is  unitary  so  the  r  associated  with  complex  polar  form  is  one. 
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5.  AV18  CPM  Synthesis 


To  couple  the  CPM  methodology  with  the  AVI 8  potential,  an  appropriate 
eoordinate  system,  basis,  Hamiltonian  and  units  system,  must  be  determined.  These 
topics  are  the  subjeets  of  the  first  four  sections  of  Chapter  5.  Seetion  5.5  defines  the 
relationship  between  the  center-of-mass  reference  frame  and  laboratory  frame  from  an 
energy  standpoint.  The  standard  methodology  for  handling  calculation  of  s-matrix 
elements  in  the  presence  of  a  1/r  potential  is  reviewed  in  seetion  5.6.  Section  5.7 
develops  a  useful  system  of  units  for  nuclear  scattering.  Seetion  5.8  gives  an  overview  of 
the  two  leading  parameterizations  for  the  non-diagonal  components  of  the  nuclear  S- 
Matrix.  Finally,  this  seetion  concludes  with  a  sample  problem  designed  to  validate  our 
time-dependent  CPM  teehnique  against  a  time-independent  analytic  solution. 

5.1.  Coordinate  System 

Since  we  are  considering  two  nucleons,  there  are  six  total  degrees  of  freedom  {xj, 
yi,  zi,  X2,  yz,  and  Zi)  whieh  must  be  taken  into  aeeount.  However,  if  we  move  to  the  space 
fixed  center  of  mass  eoordinate  system  (SF-CM)  as  illustrated  in  Figure  5.1,  the  problem 
becomes  more  tractable. 
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Figure  5.1.  Space  Fixed  -  Center  of  Mass  Coordinates. 

In  this  coordinate  system,  the  aforementioned  six  degrees  of  freedom  are  exehanged  for 
six  others,  three  pertaining  to  the  motion  of  the  center  of  mass  {Xcm,  yxm,  Zcm)  and  three 
others  that  define  the  orientation  of  the  two  bodies  around  the  eenter  of  mass  {R,  0,  (fi) . 

The  resulting  wave  funetion  can  thus  be  treated  as  the  produet  of  two  different  wave 
functions, 

¥=¥r{R,0,(I))  (5-1) 
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One  which  specifies  the  motion  of  the  center  of  mass  through  space  W  cu^^cm^ycm^'^cm) 
and  another  which  specifies  the  dynamics  of  the  particles  about  that  center  of  mass 
{R,  6,  (j)) .  This  construct  provides  an  advantage.  Since  we  are  only  interested  in  the 
physics  associated  with  scattering  dynamics  between  the  nucleons  not  where  in  space  the 
interaction  takes  place;  can  be  effectively  ignored.  Only  three  coordinates  therefore 

need  to  be  considered,  R  which  labels  the  separation  between  the  nucleons  and  0  and  (j) 
which  describe  nucleons  orientation  about  the  center  of  mass. 


5.2.  The  Hamiltonian 

The  Hamiltonian  utilized  by  the  CPM  is  expressed  using  SF-CM  coordinates  as 

[5], 


H  = 


+v  = 

2  AVIS 


2// 


-  +  H 


eff 


(5.2) 


AB 


Here,  represents  the  reduced  mass  of  the  two  body  system  and  is  the  effective 


potential  that  includes  together  with  the  kinetic  energy  tumbling  term 


2  ' 


The  radial  kinetic  energy  is  given  by 


2// 


AB 


5.3.  Basis 

A  convenient  basis  well  suited  for  representing  this  Hamiltonian  in  equation  (5.2) 
is  a  mixed  radial  coordinate,  angular  momentum  representation  labeled  with  R,  the 
distance  between  the  nucleon  pair,  the  total  angular  momentum  J,  the  total  orbital  angular 
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momentum  Z,  the  total  spin  S ,  the  total  isospin  T,  the  intrinsie  spin  a ,  and  isospinr  of 
eaeh  particle 

\y/)  =  \j  LS  (Tj  cTj  T  Tj  T2  R)  (5.3) 

Since  we  are  dealing  with  two  bodies  in  the  SF-CM  coordinate  system,  quantities  are 
independent  of  M  the  projection  of  J  onto  the  space-fixed  axis  [9].  In  this  section,  we 
examine  the  scalar  products  ( Cj  •  cTj  ,  rj  •  rj , X  •  5*,  ,  S'jj )  introduced  in  Section  3.1.3 

in  the  specified  basis  above 

{R'  T'2  t[  r  cr'  cj[  S'L'J'\V\J  LS  cj.cj^T  T,  R)  (5-4) 

to  develop  a  useful  abbreviated  basis  =  L  S  T  R'j .  Matrix  elements  are  then 

computed  to  determine  the  form  of  the  potential  in  the  abbreviated  basis  set.  From  this 
point  forward,  R ’s  presence  throughout  the  remainder  of  this  discussion  is  implied. 

A  convenient  way  to  evaluate  the  spin  scalar  product  cTj  •  (J2  is  through  the 


(5.5) 


(5.6) 


and  the  requirement  that  must  also  satisfy 

;S"|5}  =  5(;S  +  l);i"|5}  (5.7) 

to  express  the  spin  scalar  product  in  a  more  convenient  form. 
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0-1  -0-2  =2S{S  +  \)-3. 


(5.8) 


Now  as  we  have  already  seen  in  Seetion  3.1.2,  the  total  spin  S  of  the  eomposite  system  is 
either  1  for  the  symmetric  triplet  state  or  0  for  the  anti-symmetric  singlet  state.  So,  the 
spin  scalar  product  <j^  ■  <7^  in  the  coupled  basis  has  two  possible  eigenvalues,  [66,  67] 


o-j  •  o-j  I  triplet)  =  1 1  triplet) 
o-j  •  o-j  I  singlet)  =  -3 1  singlet) 


(5.9) 


The  isospin  scalar  product  has  exactly  the  same  eigenvalues  for  the  coupled  triplet 
and  singlet  states,  since  the  isospin  matrices 


= 


"0 

0 

= 

^0 

= 

n 

0^ 

vl 

0; 

y 

0; 

z 

vO 

-1. 

(5.10) 


are  analogous  to  the  Pauli  matrices  [74] 


= 


"0 

1^ 

^0 

= 

fl 

0^ 

vl 

07 

y 

07 

z 

lo 

-1. 

(5.11) 


and  the  total  isospin  T  for  the  two  nucleon  interaction  is  also  0  ox  1. 

The  total  angular  momentum  J  is  a  conserved  quantity  and  is  related  to  the  total 
orbital  angular  momentum,  L,  and  the  total  spin,  S,  via  the  relation. 


J"  =(p  +  =  +2L-S 


(5.12) 


which  may  be  manipulated  to  obtain  an  expression  T  •  S  as, 

L-S  =  '^  ^  ^ 


(5.13) 


Since  the  values  for  S  are  restricted  to  0  and  1  and 


J=L+S , 


(5.14) 
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L  can  only  have  possible  values  of  J+1,J  ox,  J-1  in  the  eoupled  basis.  The  evaluation  of 

2 

L  in  the  eoupled  basis  follows  direetly  from  the  eigenvalue  equation, 

L^\l)  =  L{L  +  \)n^\l)  (5.15) 


These  limitations  on  L  impose  restrietions  on  the  size  of  the  potential  matrix  under 
eonsideration  and  are  useful  in  the  evaluation  of  the  tensor  oomponent5'j2 . 

The  tensor  pieee  represents  the  only  non-central,  off  diagonal  eomponent  of 
the  potential.  As  discussed  in  Seetion  3.1.3,  ^  combination  of  the  sealar 

quantities  (Cj  •r){<72  •r)  and  cTj  'CTj  whieh  are  utilized  to  help  deseribe  the  observed 
asymmetry  of  the  deuteron  wave  funetion.  Wiringa  [5]  defines  the  8^2  matrix  to  be. 


R'  T'2  t[  T'  (T'2  (t[  S' L'  J'\Sy2\J  L  S  (Ty  (T2  T  Tj  Rj  = 

(-1)^'^^  [30(2T  + 1)(25  +  1)(2T '+ 1)(25 '+ 1)] 


X 


\J  S'  L' 


\2  L  S 


^  L'  2  L 
^0  0  Oy 


$2  S2 


S  2 


i  Fi  Fi 


(5.16) 


Here,  9-j  Wigner  symbol  ensures  that  there  is  no  off-diagonal  contribution  when  S'  =  0  ,  so 


here  we  ean  foeus  our  attention  on  the  3-J  eoeffieient 


^  L'  2  L 
vO  0  Oy 


and  6-J  symbol 


[J  S'  L 
[2  L  S 


when  S  =  1 .  The  possible  matrix  elements  for  Sj2  are  also  limited  when  3-J 


eoeffieient  is. 


^Xi  X2 

,0  0 


=  0 


if  Xj  +  X2  +  X3  is  odd 


(5.17) 
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In  this  case,  the  non-repetitive  eombinations  of  L  andZ '  are  redueed  to  four, 


L  =  J  +  \,  L'  =  J'+\ 
L  =  J  +  l,  L'  =  J'-\ 
L  =  J,  L'  =  J' 

L  =  J-\,  L'  =  J'-\ 


(5.18) 


The  3-J  eoeffieient  for  L  =  J  +  \,  L'  =  J'+  \  may  be  eomputed  using  the  3-J  formula. 


r/ 


/  2^ 
-M  0^ 


_ 2  [3  -/(/  +  !); _ 

[(2/-l)(2/)(2/  +  l)(2/  +  2)(2/  +  3)]‘'' 


(5.19) 


by  making  substitutions  J  =  L  =  L' ,  J€  =  0 ,  S  =  S'  =  \,  and  L  =  L'  =  J  +  1.  Similarly, 


the  6-J  symbol  may  be  eomputed  with  the  6-J  formula, 

[q  h  o'!  .  .  \  a+Z)+c 

[2  c  " 

^ _ 2[3(7-1)-4Z>(Z>  +  1)c(c  +  1)] _ 

[{2b  -  \){2b){2b  +  \){2b  +  2)(2^  +  3)(2c  -  l)(2c)(2c  +  l)(2c  +  2)(2c  +  3)]‘'" 

(5.20) 


Y  =  /)(/)  + 1) +  c(c  + 1) -a(a  + 1) 

by  making  the  substitutions  a=J,  b  =  S  =  S' ,  and  c  =  L  =  L'  =  J  +  l .  The  remaining  3-J 
eoeffieients  and  6-J  symbols  ean  be  determined  in  a  similar  fashion  by  using  these  two 
formulas  or  others  found  in  the  appendiees  of  Edmonds  [68]. 

Now  given, 

'1/2  1/2  r 
<1/2  1/2  1 
1  1  2 

(l/2||ai||l/2)  =  (l/2||a2||l/2)  =  V6 


=  1/9 


(5.21) 
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and  the  formulas,  the  matrix  elements  for  S^2  [9,  13], 


i^R'  T2  t[  T'  (j 2  (j[  S' L'  J'\S^2\J  L  S  (j^  (J2  T  Tj  T2  R^  = 


L=J+1 

L=J 

L=J-1  ^ 

L=J+1 

-2{J  +  2) 

0 

2J  +  \ 

2J  +  \ 

L=J 

0 

2 

0 

L=J-1 

0 

1 

rT 

1 

2J  +  \  2J  +  \  J 


(5.22) 

Since  we  are  dealing  with  fermions,  the  wave  function  must  be  anti-symmetric 
under  the  exchange  of  like  particles  to  satisfy  the  Pauli  Exclusion  Principle. 

y/  ( particle  1,  particle  2)  =  -y/  ( particle  2,  particle  l)  (5 .23) 

In  spherical  polar  coordinates,  is  a  function  of  r,  0,  and  (p  where  r  reflects  the  relative 
separation  between  the  two  nucleons  and  9  and  (p  define  the  orientation  of  the  two 
nucleons  in  space.  Since  the  radial  distance  between  the  two  particles  is  independent  of 
particle  exchange,  the  requirement  that  the  wave  function  be  anti-symmetric  must  be 
satisfied  by  6  and  (p .  Typically,  the  angular  dependence  is  described  by  the  spherical 

harmonics  Y^  m  >  which  have  a  parity  of  (-1)^ .  So,  even  angular  momentum 

eigenvalues  will  contribute  to  a  symmetric  wave  function  whereas  odd  angular 
momentum  eigenvalues  will  contribute  to  an  anti-symmetric  wave  function. 

For  nucleons  and  other  particles  with  spin  however,  symmetry  of  the  wave 
function  is  not  completely  determined  by  the  orbital  angular  momentum  alone,  the 
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intrinsic  spins  of  the  particles  also  contribute.  The  5*  =  1  state  where  the  nucleon  intrinsic 
spins  are  aligned, 

tt 

+  (5.24) 

is  a  symmetric  state  whereas  the  S'  =  0  state  where  the  nueleon  spins  are  aligned 

(5.25) 

is  an  anti-symmetric  state.  This  implies  that  in  order  to  produce  an  anti-symmetric  wave 
function  either,  an  even  L  state  must  be  paired  with  a  S  =  0  state  or  an  odd  L  state  must  be 
paired  with  a  S  =  1  state.  No  such  restrictions  are  imposed  on  the  wave  funetion  if  the 
interacting  pair  is  a  proton  and  a  neutron. 

Having  evaluated  all  of  the  potential  terms,  we  ean  now  examine  the  full 
Hamiltonian  in  the  abbreviated  |  T  S  J  T  ^  basis  for  a  given  value  of  J, 

{L'S'J'  T'\h\L  S  J  T)^{L'  S'  J'  T'\  —  \LSJT)  +  {L'S'J'  T'\V^j^\L  S  J  T)  (5.26) 

The  radial  kinetic  energy  operator  matrix  elements  are  diagonal  in  all  of  the  angular 
momentum  quantum  numbers  but  not  R  due  to  the  derivative, 

(!■  S'  J'  T'\^\L  S  J  T)  =  ^-!--^RS(R-R')S,.A.,S„  (5.27) 

ZfX  K  oK 

while  the  effective  potential  is  diagonal  except  when  L  =  J±\  and  T '  =  J  + 1  as  shown  in 
Figure  5.2. 
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Figure  5.2.  Evaluation  of  the  Potential  in  the  LSTJR  Basis 


Here,  non-zero  matrix  elements  on  the  diagonal  are  labeled  using  speetroseopie  notation, 
where  the  elements  have  the  form,  .  The  off  diagonal  elements  are  labeled  and 

refer  to  the  mixing  parameter  that  eouples  the  L=J-1  and  the  L=J+1  states.  The  matrix  is 
only  fully  populated  as  depicted  above  if  the  interacting  pair  is  a  proton  and  a  neutron.  If 
in  the  event  the  interacting  pair  is  two  protons  or  two  neutrons,  then  only  those  elements 
highlighted  in  red  will  be  allowed  due  to  the  Pauli  Exclusion  Principle. 

Figure  5.3  and  Figure  5.4  illustrate  how  some  of  the  potentials  are  modified 

i} 

with  and  without  inclusion  of  a  kinetic  energy  tumbling  term - ^  in  the  selected 


basis. 
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Figure  5.3.  Sample  AV18  Potentials  without  Centrifugal  Correetion 


Figure  5.4.  Sample  AV18  Potentials  with  Centrifugal  Correction 
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From  a  simple  analysis  of  the  potentials,  we  ean  gain  some  physieal  insight  about  the 
problem.  As  shown,  the  angular  momentum  eomponent  quiekly  dominates  at  small  r. 
This  implies  three  things.  One,  resonant  features  sueh  as  multiple  peaks  in  the  seattered 
wave  paeket  should  not  be  expeeted  sinee  there  is  little  evidenee  of  seeondary  wells  in 
any  of  the  plots.  Two,  phase  shifts  should  have  negative  slopes  at  low  L  and  beeome 
positive  as  L  inereases.  A  partiele  aeeelerates  entering  a  well  yielding  a  negative  phase 
shift  whereas  a  barrier  retards  it,  produeing  a  positive  phase  shift.  Finally,  sinee  the 
ability  to  aeeess  regions  defined  by  the  potential  diminishes  as  the  angular  momentum 
inereases,  values  of  J  >  5  should  eontribute  minimally  to  the  over  determination  of  the 
seattering  eross  seetion.  As  a  result,  values  of  J  >  5  do  not  need  to  be  ealeulated. 

5.4.  Coupled  Basis  Digitalization 

As  shown  in  the  previous  seetion,  the  Hamiltonian  kinetie  energy  matrix  is 
diagonal  in  the  ehosen  \L  S  J  basis  whereas  the  AVI 8  potential  is  bloek  diagonal 

(Figure  5.2).  This  mixing  of  states  in  the  potential,  sueh  as  between  |0,1,1,0) 
and|2,l,l,0)  eouples  the  wave  funetion  from  one  surfaee  to  the  other.  However,  it 
eomplieates  the  Fourier  transform  of  the  eorrelation  funetion  C{t)  from  equation  (4.25), 

T 

F(E)  =  {i//_\A^(t))=  ^  Exp 

-T 

sinee  exponentiation  of  a  non-diagonal  matrix  requires  additional  eomputational  effort. 
A  useful  teehnique  for  dealing  with  an  unwieldy  adiabatie  matrix  is  to  perform  a  unitary 


-iEt 

~Y~ 


C{t)dt 


(5.28) 
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transformation  to  an  appropriate  diabatic^^  basis  where  the  potential  is  diagonal. 
Although,  this  operation  greatly  simplifies  the  exponentiation,  there  is  a  problem.  The 
diagonal  representation  of  the  potential  operator  is  not  the  diagonal  representation  for  the 
kinetic  energy  operator.  So,  the  same  unitary  matrix  must  be  applied  to  recover  the 
original  potential’s  adiabatic  representation  prior  to  performing  the  DFT  and  the  kinetic 


energy  operator,  T .  Since,  no  coupling  will  occur  between  surfaces  in  the  diabatic  basis, 
we  must  transform  back  to  the  adiabatic  basis  and  then  repeat  the  process.  Equation 
(5.29)  summarizes  the  process  and  identifies  when  a  unitary  transformations,  U,  are 
required  [36-40]. 


(5.29) 


In  quantum  mechanics,  “adiabatic”  refers  to  a  situation  where  a  system 
undergoes  a  gradual  change,  which  leaves  the  eigenstates  unaffected.  As  an  example, 
consider  the  simple  harmonic  oscillations  of  a  pendulum.  An  “adiabatic”  change  to  the 
oscillation  would  occur  if  the  support  is  moved  slowly  enough  to  allow  the  system  to 
adapt  leaving  the  motion  unperturbed  [25]. 
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5.5.  Energy  in  the  Center  of  Mass  Frame  and  Lab  Frame 

Our  phase  shift  calculations  are  performed  in  the  SF-CM  frame  whereas 
published  data  are  generally  presented  in  the  laboratory  frame.  In  order  to  compare  our 
phase  shifts  to  published  data,  we  need  to  understand  how  to  convert  from  one  frame  to 
the  other.  This  section  investigates  the  relationship  between  the  two  reference  frames. 

Consider  the  following  pictorial  diagram  of  the  center  of  mass  and  laboratory  for  a 
simple  1-D  scattering  experiment. 


Figure  5.5.  Center  of  Mass  vs.  Laboratory  Frame 


Here,  primes  label  the  center  of  mass  frame,  v  label  pre-collision  velocities,  and  u  labels 
post  collision  velocities.  Vectors  and  (not  shown)  connect  the  positions  of  the 


particles  to  the  origin  whereas  the  vector. 


^CM  ~  ■ 


m,  +^2 


(5.30) 


connects  the  origin  to  the  center  of  mass  of  the  system.  Using  definition  (5.30)  above,  we 
can  write  the  momentum  between  the  two  systems  as, 

(5-31) 
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where  M  has  been  used  to  label  the  total  system  massM  =  mj  +  .  Sinee  the  veloeity  of 

^2  in  equation  (5.31)  prior  to  collision  in  the  lab  frame  is  zero,  a  relationship  between 
and  Vj  in  terms  of  masses  of  the  particles  can  be  established 


(5.32) 


(5.33) 


we  can  substitute  equation  (5.32)  into  (5.33)  for  and  arrive  at  a  relationship  between 
the  center  of  mass  energy  and  laboratory  frame  energy  given  by, 


^Lab  ^CM 


(5.34) 


Since  the  mass  of  the  proton  is  approximately  the  mass  of  the  neutron,  about  a  factor  of 
two  in  energy  exists  between  the  two  reference  frames.  As  the  mass  of  the  target  particle 
increases,  approaches  [70]. 


5.6.  The  1/r  Potential 

The  Scattering  Matrix  developed  in  Section  4  was  only  defined  for  potentials, 

which  decay  at  a  rate  of  or  better.  This  presents  a  problem  for  numerical  calculations 

r 

since  a  large  number  of  scattering  processes,  including  the  pp  interaction  of  interest 
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here*',  involve  the  Coulombic  potential  whieh  decays  at  a  rate  of  — .  As  r  approaches 

r 

infinity,  the  potential  still  influences  particle  motion  violating  the  boundary  condition  that 
the  wave  function  be  stationary  in  the  asymptotic  limit.  An  outline  of  Taylor’s  [51] 
methodology  for  dealing  with  the  Coulombic  tail  quandary  is  summarized  here  [30,  71- 
73], 


The  premise  of  the  solution  rests  on  the  awareness  that  pure  Coulombic  potentials 
do  not  exist  in  nature.  At  some  distance,  other  charged  particles  effectively  “screen”  the 
original  Coulombic  interaction  between  two  particles,  suggesting  an  appropriate  way  to 
handle  the  “well-behaved”  potential  is  to  impose  some  form  of  cutoff  condition.  For  this 
discussion,  we  begin  with  the  radial  component  of  Schrodinger’s  equation. 


^  L{L  +  \) 


dr 


u{r)  =  0 


(5.35) 


J 


where  we  have  used  by  using  the  standard  convention 


2juE 

~ir 


Here,  h  is  Plank’s 


constant,  E  is  the  total  energy  of  the  two-particle  system,  L  refers  to  a  particular  value  of 
angular  momentum,  and  ju  is  the  reduced  mass.  Before  proceeding  to  discuss  the 
specifics  surrounding  the  phase  shifts  induced  by  the  Coulombic  or  strong  force,  it  is 
useful  to  consider  asymptotic  free  motion. 


"  "^EMiPP)  of  equation  (2.4) 
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5.6.1.  Asymptotic  Free  Motion 


Under  the  asymptotic  free  motion  assumption,  V(r)  and 


L{L  +  \) 


are  both  set  to 


zero  and  the  radial  equation, 


dr^ 


u{r)  =  0 


(5.36) 


has  a  solution  which  is  a  linear  combination  of  sin(A:r)  and  cos(A:r) 


u{r)  =  A  sm{kr)  +  B  cos{kr) 


(5.37) 


where  A  and  B  are  arbitrary  constants.  Given  that  the  boundary  conditions  at  r  =  0 
requires  5  =  0,  the  free  solution  at  large  r  must  be  proportional  to, 

u{r)  ~  Asm{kr)  (5.38) 


5.6.2.  The  Method  of  Partial  Waves 

The  potentials  of  interest  here  are  all  central  in  nature  and  depend  on  the  relative 
distance  r  and  the  angular  momentum  L  that  is  a  constant  of  the  motion.  In  these  cases,  it 
is  useful  to  express  the  wave  function  as  a  product  of  radial  and  angular  parts  and  then 

sum  over  the  contributions  of  angular  momentum  L, 

00 

=  (5-39) 

L=0 

Here,  coefficient  refers  to  the  amplitude  of  each  partial  wave  and  refers  to  the 

spherical  harmonics  which  are  independent  of  azimuthal  angle.  To  help  understand  why 
the  use  of  partial  waves  is  particularly  useful  to  nuclear  scattering,  we  will  examine  the 
impact  parameter  b. 
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From  classical  mechanics,  we  know  that  the  impaet  parameter  b  is  defined  as  the 
perpendieular  distanee  from  the  eenter  of  foree  to  the  veloeity  veetor  of  the  ineident 
partiele.  The  impaet  parameter  b  ean  be  expressed  in  terms  of  the  ineident  partiele’s 
momentum  fjv  and  the  elassieal  angular  momentum  /  through  the  relation, 

b  =  —  (5.40) 

//V 

Scatter  for  higher  values  of  angular  momentum  are  only  possible  if  the  impaet  parameter 
is  less  than  the  range  of  the  potential  energy  between  the  two  partieles  [70]. 

This  same  eoneept  ean  be  made  applieable  to  quantum  meehanies  by  exehanging 
the  angular  momentum  /  with  the  eigenvalues  of  the  momentum  operator  L  and 
expressing  the  momentum  in  wave  veetor  form  as  p  =  tik , 


b  = 


L{L  +  \) 
IpE 
V  b" 


(5.41) 


Here  we  have  ehosen  to  express  the  wave  veetor  k  in  terms  of  energy 


to 


make  elear  the  relationship  between  impaet  parameter  and  the  energy  of  the  system. 

Since  angular  momentum  in  quantum  mechanics  can  only  assume  diserete  values,  we  can 
readily  determine  the  number  of  partial  waves,  whieh  may  interaet  with  a  finite  range 
potential.  In  Figure  5.6,  the  horizontal  dotted  line  represents  the  range  of  nuelear 
potential  and  the  eurved  lines  represent  the  impact  parameter  b{E)  for  the  first  five 
values  of  angular  momentum.  An  angular  momentum  value  of  zero  represents  a  head  on 
eollision  between  the  two  particles  is  not  shown.  From  the  diagram,  it  is  elear  that  partial 
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waves  of  Z  >  4  contribute  minimally  to  wave  function  in  equation(5.39)  and  do  not  need 
to  be  calculated  [30,  66,  74,  75]. 


Figure  5.6.  Impact  Parameters  and  the  Range  of  the  Nuclear  Force  Range  (dotted  line) 


from 


So  now,  let  us  relax  our  initial  asymptotic  assumption  and  allow  contributions 

L(L  +  \)  .  ,  , 

- r —  into  the  radial  equation. 


^  L{L  +  \) 
dr^ 


+r 


M^(r)  =  0 


(5.42) 


If  we  non-dimensionalize  the  radial  equation  by  performing  a  change  of  variables 
to  p  =  A:r ,  we  obtain  the  familiar  form  of  Bessel’s  differential  equation. 
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+p^-L{L  +  \)  u,{p)  =  Q 
dp  J 


(5.43) 


This  solution  is  a  linear  combination  of  regular  spherical  {p)  and  irregular  spherical 
functions  {p) . 


UL{p)  =  AJ^{p)  +  B^n^{p) 


(5.44) 


As  in  the  asymptotic  free  case,  the  irregular  solution  (p)  is  eliminated  through 


application  of  the  boundary  condition  m^(0)  =  0  leaving, 


Wi(p)  =  Aa(p) 


(5.45) 


and  an  asymptotically  form. 


Wi(p)~— sin  p — — 


(5.46) 


The  inclusion  of  additional  terms  in  the  argument  of  the  sine  function  suggests  that  the 
effect  of  introducing  a  potential  into  the  radial  equation  is  to  induce  a  shift  in  phase  in  the 
solution  for  large  r. 


5.6.3.  The  Coulomb  Potential 


Having  discussed  the  asymptotic  free  solution  and  the  method  of  partial  waves, 
we  now  consider  inclusion  of  the  Coulombic  potential. 


F(r)  = 


Ansp- 


(5.47) 


into  the  radial  equation. 


d'"  L{L  +  \)  2p 


h  Ansp 


+  k^  u^{r)  =  0 


(5.48) 
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Here,  labels  the  permittivity  of  free  spaee  and  e  labels  the  fundamental  charge  of  the 

electron.  As  we  desire  to  present  the  Coulombic  wave  equation  in  a  more  compact  form, 
it  is  useful  to  perform  the  same  change  of  variables  from  the  previous  section  p  =  kr  and 

2 


make  the  substitution  y  = 


pe 


Ans^hk 


^  il  L{L  +  \) 


2\ 


dp^ 


P  P 


Wi(p)  =  0 


(5.49) 


Abramowitz  and  Stegun  [76]  give  the  solution  to  this  equation  for  positive  p  ,  positive 
integer  values  of  L,  and  -oo  <  /;  <  go  as  a  linear  combination  of  the  regular  and  irregular 


logarithmic  hyper-geometric  functions  F^^{y,p)  and  G^{y,p), 


Ul  ip)  =  (r,  p) + (r,  p) 


(5.50) 


which  have  the  asymptotic  property. 


(  T  TT  ' 

p — —  +  5f^-y\n{2kp) 

V  ^  / 

(  T  71  ^ 

Gi(r,Pp^»)~cos  p — —  +  5'[-y\n{2kp) 


(5.51) 


For  the  pure  Coulombic  potential,  the  boundary  condition  at  r  =  0  requires  B^=0 .  In 
equation  (5.51)  we  see  two  new  terms,  5‘l  andyln(2^/»)  as  arguments  of  the  sine  and 


cosine  functions.  The  first,  is  the  actual  Coulomb  phase  shift  for  the  E  partial  wave. 


5^  =  argr(C  +  l-l-zy) 


(5.52) 


and  the  other  y  \n{2kp^  is  a  Coulomb  cutoff  phase  correction  designed  to  compensate  for 


the  infinite  range  of  the  nuclear  force  where  p=r  is  the  cutoff  radius.  The  rationale 
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behind  imposing  a  cutoff  condition  becomes  clear  when  we  consider  that  pure  Coulombic 
potentials  do  not  exist  in  nature.  At  some  point,  the  Coulombic  potential  between  two 
bodies  is  effectively  screened  by  a  Coulombic  potential  of  another  charged  particle.  For 
example  in  the  classic  Rutherford  scattering  experiment,  the  Coulomb  field  from  gold 
nuclei  was  completely  shielded  by  atomic  electrons  after  a  few  angstroms  [51].  Since 
now  imposing  a  cutoff  condition  seems  reasonable,  the  question  remains  at  which  point 
do  we  impose  the  cutoff  The  answer  is  that  the  cutoff  point  is  somewhat  arbitrary  as 
Figure  5.7  shows. 


Figure  5.7.  Coulombic  Phase  Corrections  at  140  fm,  1000  fin,  10000  fin,  and  \nt.  Coulombic 
corrections  to  the  phase  shifts  do  not  vary  greatly  over  small  variations  in  r. 


Coulombic  phase  corrections  simply  do  not  vary  greatly  with  small  variations  inr .  So 
here,  we  make  the  convenient  choice  of  the  cutoff  point  oip=  140  fm. 
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5.6.4.  The  Coulomb  Plus  Short  Range  Potential 


In  the  last  section,  we  considered  scattering  from  a  pure  Coulombic  potential. 

Now  we  wish  to  examine  the  implications  of  adding  a  short-range  potential  V^g(r)  to  the 
Coulombic  term, 

r(r)  =  -^  +  V,,(r).  (5.53) 

Since  we  are  adding  a  short-range  potential  to  the  Coulombic  potential,  we  would  expect 
the  solution  would  appear  as  a  small  perturbation  to  the  Coulombic  solution.  A  useful 
way  of  representing  this  perturbation  is  to  admit  contributions  from  G^(Y,Pp^^)  into  the 
asymptotic  solution, 

m(^)  ~  y4sin(^)  +  5cos(^)  (5.54) 


Ltt 

Here  we  have  made  the  substitution  6  =  +  -y\n{2kp) .  Equation  (5.54)  can  be 

written  in  terms  of  exponentials 


u(e) 


li 


(5.55) 


^  B'^ 


converted  to  polar  form  with  magnitude  A  ±  iB^^^  and  phase  v  =  tan  '  —  to  yield 


u(0) 


A  +  iB\ 


^\A-iB\ 


e  e 


li 


li 


(5.56) 


Since  the  scattering  matrix  is  unitary,  the  magnitudes  above  are  one  and  we  have, 

i(0^v)  _  -H0+V) 

u(e)~- - ^ - sin(0  +  v)  (5.57) 


li 
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This  implies  that  the  asymptotie  solution  for  a  Coulombic  potential  modified  by  a  short- 
range  potential  has  the  asymptotic  form, 


(  Lk 

Ul  (r,  )  ~  sin  p  -  —  +  cJf  -  y  In  (2^;?)  + 

V  ^  y 


(5.58) 


and  overall  phase  shift  5^  is  given  by 

-rln(2^p)  +  v^  (5.59) 

The  term  represents  the  additional  phase  shifts  introduced  by  all  short-range  forces 
and  should  not  be  confused  with  the  phase  shift  due  to  short-range  forces  alone  [51].  For 
pp  scattering,  ,  the  phase  shift  from  the  short  range  forces  alone,  is  reported.  In  order 

to  obtain  the  computed  pp  short-range  phase  shift,  ,  for  comparison  to  published  data,  a 
correction  must  be  applied  to  the  total  phase  shift  obtained  from  our  nuclear  CPM  pp 
calculations, 

v,=5,-Sl  +  rH2kp)  (5.60) 


5. 7.  Nuclear  Units 

In  order  to  perform  channel  packet  method  calculations,  a  system  of  units  must  be 
defined.  As  all  previous  calculations  utilizing  the  CPM  have  been  performed  on  the 
Atomic  Level,  the  traditional  choice  of  units  were  the  atomic  units  (See  Table  5.1) 
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Table  5.1.  Selected  Atomic  Units 


Designation 

MKS  Equivalent 

Length 

ao 

5.291772x10'"  [meters] 

Mass 

trie 

9.105953x10'"  [kg] 

Angular  Momentum 

h 

1.054571x10'"  [J]  [sec] 

Energy 

En 

4.359744x10''''  [J] 

Time 

niE„ 

2.418884x10'"  [sec] 

where  the  unit  of  length  is  based  on  the  Bohr  Radius,  the  unit  of  mass  on  the  rest  mass  of 
the  eleetron,  angular  momenta  are  measured  in  units  of  h  and  the  unit  of  energy  is  the 
hartree  (27.21 1  eF);  twiee  the  ionization  energy  of  the  hydrogen  ground  state  eleetron. 
The  atomie  unit  of  time  is  a  derived  quantity  based  on  h  and  the  hartree.  These 
traditional  units  are  ill  suited  for  nuelear  ealeulations  sinee  both  the  masses  and  energies 
involved  are  several  orders  of  magnitude  larger  than  those  eneountered  in  an  atomie 
ealculation. 

In  order  to  form  a  set  of  units,  three  base  units  must  be  eonsidered,  mass,  length 
and  time.  The  dimensions  of  the  nuclear  potential  give  insight  to  a  natural  unit  selection. 
Width  is  measured  in  femtometers  (fm)  and  the  depth  is  measured  in  MeV.  If  we  trade 
the  mass  for  energy,  we  have  only  one  more  parameter  to  fix.  Employing  the  useful 
atomic  convention  of  setting  h  to  one,  the  nuclear  unit  of  time  may  be  derived  from  the 
uncertainty  relation, 

EAt<h=l  (5.61) 

The  nuclear  unit  of  mass  may  then  be  obtained  from  our  definition  time,  length  and 
energy.  Conversion  factors  between  MKS  and  our  derived  units  are  included  in  Table  5.2 
and  values  for  nuclear  masses  are  included  in  Table  5.3. 
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Table  5.2.  Nuclear  Units 


Designation 

MKS  Equivalent 

Length 

jm 

Mass 

Angular  Momentum 

n 

Energy 

MeV 

Time 

Table  5.3.  Nuclear  Masses 


Designation 

Value 

Neutron  Mass 

Mn 

0.02412982  m 

Proton  Mass 

trip 

0.02409661  m 

Neutron-Proton  Reduced  Mass 

l^pn 

0.01205660 

f^nn 

0.01206491 

Proton  -  Proton  Reduced  Mass 

Mpp 

5.8.  S-Matrix  Parameterizations 

As  we  saw  in  the  matrix  representation  of  F^^baek  in  Figure  5.2,  a  mixing 
between  two  states  of  different  orbital  angular  momenta  is  possible  in  the  |L,S,J,T,R^ 
basis.  The  first  mixing  between  the^5'j  and  ^ states  is  shown  in  Figure  5.8. 


(L',S',J',T'|F^^|L,S,J,T)  |( 

3,1, 1,0) 

|2,1,1,0) 

r 

■N 

(o,i,i,o| 

X 

(2, 1,1,0 

V. 

J 

Figure  5.8.  Mixing  Between  ^  and  ^  States 
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The  S-Matrix  associated  with  this  2x2  block  will  also  be  a  2x2  since  there  is  a  possibility 
of  transition  from  ^  to  and  from  to  ^5*;  in  addition  to  the  possibility  of 


reflection  back  into  original  states.  Although  many  parameterizations  of  this  2x2  S- 
matrix  block  could  exist,  only  two  are  repeatedly  found  in  the  literature;  the  Stapp, 
Ypsilatis,  and  Metropolis  “bar”  [18]  and  the  Blatt/Biedenham  “eigen”  [14] 
parameterizations. 

Blatt  and  Biedenham  proposed  viewing  the  symmetric  2x2  matrix  in  its  most 
natural  eigen-basis  form, 


S  = 


in(^.)  cos(^,) 


cos 

sin 


^  ^  US, 

e 


0 


r 


2iS 


•i.j-i 


COS 

-sin 


V 


in(^J  cos(^J 


(5.62) 


'  e  '  cos 


S  = 


(  \  ,  20,.,  1  .  2 

iSjj  +  e  '  sm 


( cos  ( sin  ( ^ 


-cos 


V 


(^Sj  ]  sin  ^  ‘  -  j 


20,. 


''0-1 


COS 


(^y) 


+  sin^  i^Sj  ^ 


(5.63) 


This  unitary  transformation  can  be  thought  of  as  a  series  of  rotations  where  U  is  a 
function  of  the  mixing  parameter,  Sj  and  S  are  the  phase  shifts  of  the  eigen-state 


associated  with  a  specific  /  state.  Although  on  the  surface  it  would  seem  to  be  the  most 
natural  basis  for  presenting  the  phase  shift  information,  difficulties  exist  in  separating  the 
Coulombic  and  strong  contributions  for  the  pp  triplet  state.  Thus,  most  phase-shift 


analyses  are  performed  with  the  “barred”  parameterization. 
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(5.65) 


cos  (  I  sin  2s  j  ^  ^ 

z  sin  ( 2s  j  ^  ^  cos  ( 2 


where  the  three  parameters,  and  Sj ,  are  the  parameters  that  are  generally 

reported.  Since  the  unbarred  phase  shifts  are  required  in  the  determination  of  the  total 
cross  section,  the  following  transformations  between  representations  are  often  useful. 


^i=j-i  +  ^i=j+i  -  ^i=}-\  +  ^i=}+\ 

.  tan  (2^) 

s.n(4,-,-4=,.,)  =  ^^ 

.  .  sin  (2^) 

s.n(4,-,-4=,.,)  =  ;^ 


(5.66) 


5.9.  1-D  Sample  Calculation  -  Square  Well  Potential 

It  is  often  useful  to  first  calculate  and  validate  a  new  technique  against  a  problem 
that  has  a  known  analytic  solution.  In  this  case  we  use  an  asymmetric  square  well  of 
roughly  the  same  dimensions  as  the  actual  nuclear  potentials  under  investigation  (Figure 
5.9).  This  test  verifies  operation  of  all  CPM  modules  exeept  the  2x2  module.  Section 
5.9.1  will  first  develop  the  general  form  for  an  analytic  solution  for  a  square  well  of 
arbitrary  dimensions  then  apply  the  formalism  to  the  test  case.  Section  4.9.2  will  develop 
the  analogous  CPM  solution.  Results  are  eompared  in  section  4.9.3. 

3000  MeV  if  X  <  0.65  fim 

F(x)  =  Jo  MeV  ifx  >  1.65  fin  [  (5.67) 

-100  MeV  otherwise 
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Figure  5.9.  Square  Well  Approximatiou  to  *So  Poteutial.  The  dimeusious  of  the  well  were  adjusted 
to  the  *So  poteutial  to  provide  a  method  of  validatiug  the  Chauuel  Packet  Method  algorithm  agaiust  a 
known  solution  prior  to  propagating  with  actual  AV18  model. 


5.9.1.  Analytic  Solution 


Consider  the  following  generic  square  well  problem, 


V{x)  =  \ 


V,  if  X  (  a 
V„  if  a  <  X  < 


V,„  if  b  <  X 


b 


Region  I 
>  Region  II 
Region  III 


(5.68) 
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Here,  solutions  to  the  time  independent  Sehrodinger’s  equation,  Hy/  =  {f  +  V^y/  =  Ey/  , 
for  each  region  are  given  by, 

{x\y/')  =  y/j^{x)  =  AjfCxp +  Bj^exp ,  (5.69) 


where, 

_^2p(E-V,) 

^  P  -  -) 


(5.70) 


and  R  refers  to  regions  I,  II,  or  III  Since  the  Schrodinger’s  equation  is  a  second  order 
differential  equation,  continuity  must  exist  at  the  boundary  between  regions  I/II  and 
regions  II/III  for  both  y/  and  its  first  derivative.  Thus,  we  have  the  pair  of  transfer 
matrices. 
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(5.72) 


where 


a=—  1h - Qxpi^ia  (k2-k\)) 
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(5.76) 


These  matriees  ean  be  eombined;  eliminating  two  of  the  six  unknown  eoeffieients  to 
ereate  a  seattering  matrix  between  regions  I/III, 


bA  s,, 


)  W  2 1  “^22  /  V  ^ III 


(5.77) 


where  the  eomponents  of  S  are  expressed  in  terms  of  a,/3,c,  and  d  and^l/  {Bm)  refers  to 
the  ineident  amplitude  from  the  left  (right)  and  Am  (Bj)  refers  to  the  outgoing  amplitude 
Am  to  the  left  (right). 

Now  if  we  select  our  incident  wave  packet  to  be  incident  from  the  right,  Ajis  zero 
and  only  two  scattering  matrix  elements  must  be  considered;  the  transmitted  amplitude  to 
the  left,  S'i2  and  reflected  amplitude  to  the  right  S22, 


{ac  +  /3^d){ac  +  /3d*) -{ad*  +  /3*c){ad  +  /3c) 
{ac  +  pd*) 


(5.78) 
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{ad  +  /3c  ) 
{ac  +  f3dy 


(5.79). 


5 


22 


So  far,  the  discussion  has  been  limited  to  the  generic  square  well.  Our  asymmetric  square 
well  approximation  for  the  nuclear  well  from  Figure  5.9,  however,  has  a  barrier.  Here, 
the  potential,  F, ,  at  x  =  a  is  considered  to  be  sufficiently  large  enough  to  minimize 
transmission.  is  zero.  With  only  one  remaining  element  to  consider,  the 
determination  of  the  reflected  analytical  phase  shift  is  straightforward  via  application  of 
equation  (5.79), 
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5.9.2.  CPM  Solution 


For  CPM  calculations,  a  useful  choice  is  a  complex  Gaussian  for  y/. 
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(5.82) 


because  it  is  simple  in  form  and  can  be  analytically  propagated.  Since  CPM  calculations 
involve  both  the  coordinate  and  momentum  spaces,  choosing  this  form  provides  an 
analytic  means  of  verifying  numerical  propagation  techniques  in  the  absence  of  an 
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interaction  potential.  In  the  above  equation,  t  labels  to  time,  h  labels  Plank’s  constant, 
fj.  labels  the  reduced  mass  of  the  two  nucleons,  k  labels  the  wave  vector,  labels  the 

initial  packet  position,  ko  labels  the  initial  momentum,  and  Axg  labels  the  wave-packet 
spread. 

Although  any  position  on  the  coordinate  grid  may  be  chosen  for  the  initial 
complex  Gaussian  product  (reactant)  wave-packets,  we  choose  to  center  both  wave- 
packets  peaks,  designated xio  andx2o  respectively,  at  30  fm  as  shown  in  Figure  5.11  and 
proceed  to  find  optimal  values  for  the  wave-packets  propagation. 
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Figure  5.11.  Scaled  Wave  Packets  at  t=0  in  Coordinate  Representation. 


If  we  consider  the  time-independent  form  of  equation  (5.82), 
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(5.83) 


(5.84) 


it  is  evident  that  a  deerease  in  eoordinate  spaee  wave-paeket’s  spread,  Ax , corresponds  to 
an  increase  in  momentum  space  wave-packets  spread.  Equations  (5.83)  and  (5.84)  were 
plotted  using  Ax  values  of  2.474 fm,  1.814  fm,  and  1.154  frn  to  help  identify  an 
appropriate  choice  for  Ax  given  our  chosen  1600 fm  coordinate  space  and  our  desired  0- 
400  MeV energy  range.  These  are  shown  in  Figure  5.12. 
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Figure  5.12.  The  Coordinate  Representation  (top  frame)  and  Momentum  Representation  (bottom 
frame)  of  the  Same  Three  Wave  Packets.  A  decrease  in  coordinate  space  wave  packet  spread 
constitutes  an  increase  in  the  wave  packet  spread  in  the  momentum  space  complementary  wave 
packet.  Both  representations  must  be  monitored  to  ensure  the  wave  packets  fit  on  the  available  grid 
space. 
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For  a  Ax  of  2.474 fm,  the  wave-packet  amplitude  tends  to  zero  at  k  values  of 
approximately  -2.0  [l//m]  and  -0.5  [l//m].  The  effective  energy  range  may  be  estimated 
from  these  k  values  and  the  reduced  mass,  /u  ,  of  the  interacting  pair  via, 

£  =  (5.85) 

l/U 

So  for  this  particular  choice  of  Ax  ,  an  effective  energy  range  between  20  MeV  and  330 
MeV  is  obtained.  Although  this  is  a  perfectly  acceptable  selection,  what  is  desired  is  a 
slightly  wider  energy  range.  So,  a  narrower  Ax  would  be  a  better  choice. 

Now,  consider  Ax  of  1.154  fm.  Figure  5.12  illustrates  two  problems  with  this 
choice.  The  first,  rather  obvious,  problem  is  that  the  wave-packet  does  not  fit  on  the 
selected  grid  space.  The  second  is  a  little  more  subtle.  Note  that  the  right  side  of  the 
wave-packet  tends  to  zero  past  the  point  of  zero  momentum  in  the  positive  momentum 
region.  Since  this  wave-packet  contains  both  positive  and  negative  momentum 
components,  the  wave-packet  will  split.  The  negative  components  will  probe  the 
potential  whereas  the  positive  components  will  leave  the  interaction  region.  In  other 
words,  the  reflection  amplitude  of  the  S-matrix  elements  will  no  longer  be  one  and  the 
information  necessary  to  compute  the  correct  phase  shift  will  be  lost.  Here,  the  most 
useful  choice  for  Ax  is  1.814ym  as  it  both  fits  on  the  grid  space  and  has  an  energy  range 
from  about  5  MeV  to  about  420  MeV.  The  remaining  parameters  pertaining  to  this 
calculation  are  summarized  in  Table  5.4. 
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Table  5.4.  1-D  Scattering  Parameters 


Parameter 

Quantity 

Units 

Mass 

0.01201 

X  max 

1600 

fm 

k  max 

2.56 

1/ fm 

Coordinate  Points 

213 

Coordinate  Step  Size 

0.195 

fin 

T 

max 

20 

A 

Temporal  Points 

218 

Time  Step 

7.63x10-5 

A 

Initial  Coordinate  Reactant  Wave-Packet  Position  (  xIq  ) 

30 

fin 

Initial  Coordinate  Product  Wave-Packet  Position  {x7,q  ) 

30 

fin 

Initial  Reactant  Momentum  (  kl^ ) 

-1.253 

1/ fm 

Initial  Product  Momentum  (  klg  ) 

1.253 

1/ fm 

Packet  Spread  (  Ax  ) 

1.814 

fin 

Since  no  portion  of  the  wave-paeket  resides  within  the  potential,  the  requirement  to  first 
propagate  the  wave-paekets  to  the  asymptotie  limit  under  the  ehannel  asymptotie 
Hamiltonian  and  baek  under  the  full  ehannel  Hamiltonian  is  eliminated.  For  this 
problem,  these  are  our  Moller  states  of  equation  (4.8). 

To  ealeulate  S-Matrix  elements,  a  FORTRAN  wave-packet  propagation  eode  was 
developed.  Snapshots  of  the  reaetant  wave-paeket’s  time  evolution  are  provided  in 
Figure  5.13-  Figure  5.16  with  absolute  values  shown  for  elarity. 
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Figure  5.13.  Scaled  Evolving  Wave  Packet  at  0.122  .  The  evolving  wave  packet  has  left  its  initial 

position  and  has  traversed  half  the  distance  to  square  well.  Absolute  Values  are  Plotted 


Figure  5.14.  Wave  Packet  Propagation  at  0.244  .  The  evolving  wave  packet  has  begun  interacting 

with  the  square  well.  No  evidence  of  bifurcation  is  present. 
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Figure  5.15.  Wave  Packet  Propagation  at  0.366  .  The  higher  energy  components  of  the  wave 

packet  have  left  the  interaction  region  and  are  returning  to  the  initial  position  so  that  the  correlation 
with  the  product  wave  packet  may  he  computed. 


Figure  5.16.  Wave  Packet  Propagation  at  0.977.  With  the  exception  of  the  lowest  energy 
components,  the  wave  packet  has  all  hut  left  the  interaction  region. 
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In  the  first  frame,  the  wave-paeket  peak  has  advaneed  under  the  full  Hamiltonian 
from  t=0  approximately  \3fm  toward  the  potential.  At  this  point,  no  useful  data  has 
been  eolleeted  as  the  wave-paeket  has  not  yet  entered  the  interaetion  region.  By  the 
seeond  frame  (0.244  ),  the  higher  momentum  eomponents  of  the  wave  paeket  have 

reaehed  the  interaetion  region  and  eollided  with  the  barrier.  Here,  the  total  energy  within 
the  well  area  elimbs  as  kinetie  energy  is  exehanged  for  potential  energy.  No  evidenee  of 
bifureation  is  apparent.  Later  at  0.366  ,  the  wave-paeket  after  eollision  is  shown 

beginning  to  exit  the  area  having  eolleeted  information  about  the  potential.  At  eaeh  time 
step,  the  eorrelation  funetion  C{t)  is  determined  by  taking  the  sealar  produet  between  the 

evolving  wave-paeket  and  the  produet  Moller  state.  By  0.977  (Figure  5.16),  the 
ealeulation  is  essentially  eomplete  as  all  but  the  lower  momentum  eomponents  have 
exited. 


From  the  eorrelation  funetion  C(t)  between  the  evolving  reaetant  state  and  the 
stationary  produet  state  reeorded  at  eaeh  time  step,  we  obtain  the  eorrelation  funetion  as  a 
funetion  of  energy  by  taking  the  Fourier  Transform 

T 

{¥-\ME))=\  Exp 


-iEt 

~Y~ 


C{t)dt 


(5.86) 


S-matrix  elements. 


-k^+k 


{E)  = 


h\k{E) 


2npri*_[-k{E))ri^[+k{E)) 


{^M^E)) 


(5.87) 


are  then  eomputed  by  expressing  both  k  and  the  reaetant/produet  initial  expansion 
eoeffieients  77+  [k]  in  the  momentum  representation  as  a  funetion  of  energy. 
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(5.88) 


k{E) 


v^ 


h 


r]{k{E))=  — Ax^  exp[^-(A:(£')-A:o)^»Ax^ +zA:(£')Xq]  (5.89) 

y 

The  ?7+  [A:]  s  are  determined  analytieally  from  the  Gaussian  wave  function  which  defined 
our  initial  states  at  t  =  0  (equation(5.82)).  Here,  k^  and  x,,  refer  to  the  respective  initial 
conditions  for  each  wave-packet  from  Table  5.4. 

The  overlap  between  initial  expansion  coefficients  77+  ,  the  wave 

vector  (£') ,  and  the  Fourier  transform  of  the  correlation  function 
shown  in  Figure  5.17. 


Figure  5.17.  Scaled  Overlap  of  equation  (5.87)  components  as  a  function  of  Energy.  Since  the 
correlation  function  can  rise  much  faster  than  the  product  of  the  expansion  coefficients  a  division 
error  is  possible  at  low  energies.  How  to  deal  with  low  energy  calculations  is  discussed  in  Section  6. 
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Clearly,  a  division  error  resulting  from  dividing  a  larger  number  by  a  smaller  number 
exists  at  low  energies  sinee  the  Fourier  Transform  of  the  eorrelation  function  can  rise 
much  faster  than  the  product  of  the  initial  expansion  coefficients  in  the  denominator  of 
equation  (5.87).  Methods  for  dealing  with  low  energy  resolution  will  be  discussed  in 
Chapter  6  where  low  energy  phase  shifts  are  investigated. 

An  issue  may  arise  during  the  calculation  of  the  correlation  function  between  the 
dispersing  wave-packet  and  the  periodic  nature  of  the  Fourier  transform.  When  the 
dispersing  wave-packet  crosses  the  grid  boundary,  the  transform  will  cause  the  exiting 
wave-packet  to  reemerge  from  the  opposite  boundary  and  continue  propagating  across 
the  same  grid  space.  If  this  occurs  before  the  dot  product  of  the  reactant  wave-packet’s 
trailing  edge  and  product  Moller  state  is  negligible,  the  leading  edge  will  introduce  non¬ 
physical  anomalies  into  the  correlation  calculation.  This  non-desirable  situation  can 
precipitate  construction  of  larger  and  larger  coordinate  grid  spaces  to  compensate,  which 
is  not  computationally  advantageous.  Fortunately,  a  proven  workaround  does  exist,  the 
inclusion  of  an  absorbing  boundary  condition. 

An  absorbing  boundary  condition  is  a  useful  tool  often  employed  to  compensate 
for  this  phenomenon.  To  employ  this  technique,  the  interaction  potential  operator,  F, ,  of 
the  Hamiltonian  is  augmented  by  an  additional  imaginary  piece, 

F  =  K±it,  (5.90) 

which  attenuates  the  wave  function  as  it  approaches  the  edge  of  the  grid  space.  Inclusion 
of  this  non-physical  condition  does  not  have  an  adverse  effect  on  the  calculation  unless 
the  absorbing  boundary  overlaps  the  product  Moller  state.  If  the  absorbing  boundary  did 
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overlap  the  produet  Moller  state,  the  absorbing  boundary  would  artificially  reduce  the 
evolving  wave-packet/product  Moller  state  correlation  function.  For  this  research,  we 


utilize  the  same  exponential  functional  form  successfully  employed  by  Calfas  and  Weeks, 
[39]  and  later  Niday  [41]  for  three  body  atomic  problems, 

-(x-Xo)" 


V^=  A  exp 


B 


(5.91) 


where  the  amplitude  (A)  and  width  (B)  parameters  can  be  tweaked  to  maximize 
attenuation  and  minimize  reflection  at  the  boundary.  Here,  Xg  was  set  to  x^^^  =1600 fm 

whereas  the  values  for  3000  MeV  and  B=27195 fm  were  determined  by  trial  and 

error. 


5.9.3.  Analytic  -  Numerical  Comparison 

Given  the  information  from  the  previous  section,  determining  scattering  matrix 
elements  from  equation  (5.87)  is  uncomplicated.  Figure  5.18  compares  the  scattering 
amplitude  analytic  to  the  numerical  solution  whereas  Figure  5.19  contrasts  the  phase  shift 
solutions.  The  analytic  solutions  are  represented  as  circles. 
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Figure  5.18.  Analytic  vs.  Numerical  S-Matrix  Amplitude.  A  scattering  amplitude  of  unity  is 
expected  since  the  entire  wave  form  should  he  reflected  from  the  potential  harrier  at  an  r  ~  0.5  fin . 
The  apparent  ringing  low  energy  is  associated  with  a  division  error  and  can  he  remedied  hy  an 
appropriate  choice  of  a  low  energy  wave  packet.  The  lower  sampling  rate  failed  to  converge  to  one 
due  to  a  poor  choice  in  coordinate  step  size. 


Intuitively  we  would  expect  the  scattering  amplitude  to  be  one  since  the  entire 
waveform  should  be  reflected  from  the  potential  barrier.  The  failure  of  the  lower 
sampling  rate  to  meet  this  expectation  stems  from  a  poor  choice  of  step  size.  Here  the 
lower  sampling  rate  corresponds  to  about  two  samples  per  fm  which  implies  that  the 
square  well  appears  to  have  more  of  a  trapezoid  shape  than  that  of  a  square.  If  we 
increase  the  sampling  rate  by  four  orders  of  magnitude,  the  expected  scattering  amplitude 
of  one  is  obtained.  Note  that  both  solutions  also  exhibit  a  ringing  effect  at  the  upper  and 
lower  energy  limits,  which  is  independent  of  sampling  rate.  This  ringing  is  a  result  of  the 
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division  by  the  product  of  the  expansion  coefficients  [^(-E)]  in  equation(5.87),  which 


are  small. 
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Figure  5.19.  Analytic  vs.  Numerical  Phase  Shift.  The  phase  shift  is  more  susceptihle  to  a  poor  choice 
in  coordinate  step  size  than  the  scattering  matrix  amplitude.  At  the  lower  sampling  rate,  the  square 
well  appears  to  he  more  of  a  trapezoid  than  a  square  well  and  fails  to  converge  to  the  analytic 
solution.  Since  the  actual  potentials  are  free  of  sharp  discontinuities,  lower  sampling  rates  should  he 
possible  when  calculations  with  the  AV18  are  performed. 


Shifting  the  focus  now  to  Figure  5.19,  the  asymmetric  square  well  phase  shift 
exhibits  a  greater  sensitivity  to  the  trapezoid  appearance  of  the  well  than  the  scattering 
matrix  amplitude.  At  the  low  2  rate,  the  analytic  and  calculated  differ  by  about  15 
degrees.  As  the  sampling  rate  is  increased  to  2'^,  the  phase  shift  aligns  with  the  analytic 
solution.  These  issues  should  not  pose  a  problem  with  the  actual  AV 1 8  potential  since  all 
surfaces  are  smooth  and  free  of  sharp  discontinuities  within  the  range  of  the  nuclear 
strong  force. 
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6.  Presentation  of  Results 


The  results  of  this  researeh  are  presented  in  three  main  seetions.  The  first 
examines  the  'So  potential.  The  'So  potential  is  the  simplest  two-nueleon  case  where  all 
three  nucleon-  nucleon  combinations  are  possible  (pp,  pn,  and  nn).  The  lack  of  a 
centrifugal  component  makes  it  possible  to  position  the  wave-packets  close  to  the 
potential  well  on  a  small  grid  space.  As  a  result,  computational  time  is  significantly 
reduced.  Subsequent  calculations,  however,  become  more  complex  as  the  effect  of 
angular  momentum  can  no  longer  be  ignored.  Coordinate  grid  spaces  must  be  expanded 


to  compensate  for  the  longer  range  of  the  centrifugal  term  ■  Section  6.2 

2juR 

assesses  the  computational  impact  of  angular  momentum  through  J=4.  Finally,  a  more 
efficient  method  for  calculating  nuclear  phase  shifts  that  exploits  the  innate 
characteristics  of  the  potential  is  developed  in  Section  6.3. 


6.1.  The  ^So  Potential 

For  the  'So  potentials,  there  exists  one  effective  potential  surface  for  each  possible 
nucleon  combination,  one  for  neutron-neutron,  one  for  neutron-proton,  and  one  for 
proton-proton.  These  are  shown  in  Figure  6. 1 . 
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Figure  6.1.  The  *So  AV18  potentials.  The  first  frame  shows  the  dimensions  of  the  potentials.  The 
second  frame  highlights  the  0.3  MeV  Coulomhic  harrier  on  approach  to  the  well.  The  third  frame 
illustrates  the  slight  differences  in  maximum  well  depth  between  the  three  potentials. 
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A  noticeable  surface  similarity  between  the  plots  is  immediately  apparent.  In  fact,  the 
only  distinguishing  feature  between  them  is  the  slight  inerease  in  the  pp  potential 
approaehing  the  nuelear  well.  The  small  differenees  in  the  last  frame  of  Figure  6.1  are 
due  to  the  charge  dependence  of  the  strong  foree  whereas  the  slight  inerease  in  the  middle 
frame  is  a  consequence  of  Coulombie  repulsion.  Despite  the  surfaee  similarity,  some 
general  information  about  the  phase  shift’s  qualitative  aspeets  may  be  inferred: 

1 .  Secondary  wells  are  not  evident  so  portions  of  the  wave  packet  eannot  beeome 
delayed  in  exiting  the  interaetion  region  and  generate  resonant  features  sueh 
as  seeondary  peaks  in  the  amplitude  of  the  seattered  wave  paeket. 

2.  The  phase  shifts  should  be  similar  sinee  only  minor  differenees  exist  between 
the  three  'S'g  potentials. 

3.  All  are  attraetive  (negative)  therefore,  the  resulting  phase  shifts  should  be 
primarily  positive  .[77]. 

As  shown,  the  AV 1 8  potentials  are  smooth  and  free  of  any  diseontinuities  within  the 
interaction  region.  This  is  not  the  ease  outside  the  interaction  region  where  only  the  nn 
and  np  potentials  tend  to  zero  well  before  the  \A0fm  eut-off.  Here,  the  pp  potential  is 

still  slightly  positive  at  termination  due  the  slow  —  Coulombie  tail  deeay  (Figure  6.2). 

r 
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Figure  6.2.  The  *So  potentials  at  cut-off.  The  0.01  MeV  step  in  the  proton-proton  potential  should 
have  minimal  impact  on  the  phase  shift  calculation. 

Although  the  potential’s  step  at  140  fm  has  minimal  impact  on  the  phase  for 

the  ^So  potential,  we  still  choose  the  initial  coordinate  position  of  the  wave-packets  to  be 
well  outside  the  interaction  region  at  400 fm  to  minimize  any  possible  low  energy 
interference.  Subsection  6.1.1  discusses  the  initial  guess  for  coordinate/temporal  step 
size.  Refinement  of  parameters  and  test  for  convergence  of  these  step  sizes  are  the 
subject  of  Subsection  6.1 .2.  Subsection  6. 1 .3  describes  a  methodology  for  determining 
an  appropriate  set  of  absorbing  boundary  conditions.  Finally,  Subsection  6.1.4  presents 
the  results  of  the  CPM  for  the  ^So  potential. 
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6.1.1.  Selection  of  Initial  Coordinate  and  Temporal  Step  Size 

Obviously,  the  choice  of  coordinate  and  temporal  grid  parameters  size  can 
significantly  influence  the  computational  cost  of  this  technique.  A  poor  choice  may 
never  provide  any  useful  information  nor  prove  to  be  more  advantageous  than  traditional 
time-independent  methods.  In  this  section,  an  outline  on  how  to  estimate  an  appropriate 
step  size  is  presented. 

Nyquist’s  sampling  theorem  provides  a  useful  prescription  for  establishing  an 
initial  guess  for  an  appropriate  coordinate  step  size.  The  theorem  stipulates  that  in  order 
to  reconstruct  an  original  analog  signal  digitally,  the  sampling  rate  must  be  twice  the 
highest  frequency  observed  in  the  original  analog  waveform.  Here,  the  interesting 
features  of  the  potential  well  illustrated  in  Figure  6.1  only  exist  over  the  first  few 
femtometers  so  at  least  two  samples  per  fm  are  required  to  provide  some  reasonable  phase 
shift  information.  Given  our  initial  choice  of  a  1600 fm  coordinate  grid  and  the  power  of 
2  FFT  condition, 

^™  =  0391[>,]  (6.1) 

is  the  first  possible  step  size  candidate  that  meets  the  criterion. 

An  initial  estimate  for  the  temporal  step  size  can  just  as  easily  be  determined 
using  our  nuclear  convention  form  of  the  energy-time  uncertainty  relation, 

M  =  -^,  (6.2) 

E 

max 

Here,  E^ax  is  the  sum  of  the  potential  energy  given  by  the  depth  of  the  nuclear  well  (115 
MeV)  and  the  kinetic  energy  computed  from  the  highest  value  of  k  contained  within  the 
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wave  packets  (543  MeV).  For  convenience,  we  set  =  lOOOMeF  which  establishes 


0.001  as  an  upper  limit  for  Ar  .  Given  that  it  takes  about  80  for  the  amplitude  of  the 

correlation  function  between  the  evolving/product  wave-packet  to  fall  below  10'^,  we 
have, 

80  T 

0.001  r,  >  =  0.0006104  ,  (6.3) 

for  the  temporal  step  initial  guess.  These  choices  again  represent  initial  guesses  and 
should  not  be  assumed  to  be  optimal.  Refinement  of  these  parameters  is  the  subject  of 
the  next  section. 

6.1 .2.  Parameter  Refinement  and  Convergence  Tests 

In  this  section,  we  seek  to  ascertain  the  most  appropriate  coordinate  and  temporal 
step  size  to  achieve  suitable  results  while  minimizing  computational  effort.  Beginning 
with  the  initial  guesses  identified  in  the  last  subsection,  systematic  convergence  tests 
were  performed  using  the  potential.  To  test  for  coordinate  step  size  convergence, 
the  time  step  was  held  constant  while  Ar  was  successively  reduced  by  factors  of  two  until 
the  phase  shift  converged  on  a  stable  phase  value.  As  shown  in  Figure  6.3,  the  threshold 
for  coordinate  step  convergence  is  obtained  at  a  Ar  of  0.195  fm. 
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Figure  6.3.  Delta  r  Convergence  for  Fixed  Delta  t.  A  step  size  smaller  0,195  fin  is  unnecessary  since 
the  phase  shift  has  already  converged  to  a  stable  value. 


Here  we  have  chosen  to  zoom  into  a  small  energy  range  since  small  variations  in 
the  phase  shift  would  not  be  readily  apparent  on  the  0-350  MeV  scale.  This  choice  for 


Ar  fixes  the  maximum  momentum  supported  by  the  grid  space  to  be  ±2.56  —  since 

fm 


—  =  .  An  optimal  value  for  At  (0.0003 1  )  was  determined  similarly  by  holding 

Ar 

this  new  Ar  constant  while  reducing  At  by  factors  of  two. 
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6.1.3.  Absorbing  Boundary  Conditions 


Obviously,  the  amount  of  time  required  to  perform  CPM  ealculations  is  a  funetion 
of  grid  size;  the  larger  the  grid,  the  more  time  that  is  required  to  eomplete  the  ealeulation. 
For  these  and  subsequent  higher  order  J  ealeulations,  insertion  of  an  appropriate 
boundary  eondition  is  eritieal  as  it  not  only  reduees  grid  size  and  eomputation  effort  but 
also  prevents  faster  moving  wave-paeket  eomponents  from  reaehing  the  boundary  before 
the  slower  eomponents  have  probed  the  potential. 

Seleetion  of  an  appropriate  boundary  hinges  on  proper  seleetion  of  the  boundary 
width  B  and  amplitude  A  parameters  of  equation  4.59, 


V^=  A  exp 


B 


(6.4) 


An  appropriate  selection  was  made  by  computing  the  correlation  function  between  the 
reflected  portion  of  an  evolving  wave-packet  incident  on  the  boundary  and  negative 
momentum  product  state.  A  negative  momentum  product  state  refers  to  a  wave-packet 
with  the  same  exact  parameters  as  the  initial  reactant  state  save  one,  the  sign  of  the  initial 
momentum.  The  dissipation  of  a  propagating  wave-packet  by  an  absorbing  boundary  is 
shown  in  Figure  6.4  . 
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Figure  6.4.  Dissipation  of  an  Evolving  Wave-Packet  by  an  Absorbing  Boundary.  The  choice  of  an 
appropriate  absorbing  barrier  quickly  dissipates  the  propagating  wave  packet.  An  absorbing 
boundary  improves  computational  efficiency  by  reducing  the  coordinate  grid  space  necessary  to 
perform  a  time  dependent  scattering  calculations. 


By  150  most  of  the  incident  wave  packet  has  been  dissipated.  The  absorption  of  the 
boundary  was  maximized  by  varying  the  parameters  by  trial  and  error  until  the 
correlation  function  between  the  wave  packet  reflected  off  the  barrier  and  the  negative 
momentum  product  state  was  minimized.  Using  this  approach,  an  amplitude,  A,  of  3000 
MeV and  a  width,  B,  of  329.063  fm  were  found  to  be  acceptable.  All  'So  CPM  parameters 
are  summarized  in  Table  6.1. 
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Table  6.1.  Initial  Parameters  used  for  nuclear  CPM  calculations 


Parameter 

Quantity 

Units 

Xmax 

1600 

fin 

kmax 

2.56 

1/ fin 

Neutron-Proton  Reduced  Mass,  jUpn 

0.01205660 

niy 

Neutron-  Neutron  Reduced  Mass, 

0.01206491 

niy 

Proton  -Proton  Reduced  Mass,  Upp 

0.01204830 

niy 

Coordinate  Points 

Coordinate  Step 

0.195 

fin 

Time 

80 

A 

Temporal  Points 

218 

Time  Step 

0.000305 

A 

Initial  Reactant  Coordinate,  xlg 

400 

fin 

Initial  Product  Coordinate,  x2q 

400 

fin 

Initial  Reactant  Momentum,  kl^ 

-1.253 

1/ fm 

Initial  Product  Momentum,  kl^ 

1.253 

1/ fm 

Packet  Spread,  Ax 

1.814 

fin 

Barrier  Amplitude  (A) 

3000 

MeV 

Barrier  Width  (B) 

329.063 

fin 

6.1.4.  The  ^So  Phase  Shifts 

The  absorbing  boundaries  and  wave-packets  were  configured  based  upon  the 
Table  6.1  parameters  for  propagating  on  the  three  'So  potentials.  For  each  case,  a 
correlation  function,  S-matrix,  and  phase  shift  were  computed  for  comparison  to  the  ANL 
published  phase  shift  data.  The  CPM/ ANL  '  So  phase  shift  comparisons  are  shown  in 
Figure  6.5  through  Figure  6.7.  ANL’s  data  was  only  available  at  the  discrete  energies  of 
1,  5,  10,  25,  50,  100,  150,  200,  250,  300,  and  350  MeV. 
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Figure  6.5.  *So  CPM  vs.  ANL  neutron-neutron  phase  shift. 


Figure  6.6.  *So  CPM  vs.  ANL  proton-neutron  phase  shift 
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All  of  these  plots  are  smooth,  free  from  resonant  features,  and  exhibit  a  positive  slope 
that  is  indicative  of  an  attractive  potential.  Here,  the  neutron-neutron  and  neutron-proton 
phase  shifts  achieve  extremely  good  agreement  with  ANL’s  published  results,  whereas 
the  proton-proton  case  does  not  due  to  the  long-range  Coulombic  tail.  Application  of  the 
Coulombic  correction  term  from  section  5.6.4  is  discussed  in  the  next  subsection. 

6.1 .5.  The  ^So  pp  Phase  Shift  Coulombic  Correction 

As  we  saw  in  the  last  section,  phase  shift  did  not  agree  well  with  the 
published  ANL  results.  Here,  we  apply  the  Coulombic  correction  developed  in  Section 
5.6  to  obtain  dL^S^pp  solution  similar  to  ANL.  Agreement  with  ANL  may  be  obtained  by 
correct  application  of  the  Coulombic 
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cr^  =  arg  r(l  +  iy  +  L)  , 


(6.5) 


and  the  cut-off  correction, 


where. 


y\n{2kp). 


(6.6) 


7  =  ^, 


AnSohk 


r 


(6.7) 


Here,  p  is  the pp  reduced  mass,  E  is  Energy,  ti  is  Plank’s  constant, is  the  coordinate 
grid  Coulombic  cutoff  point  of  140  fm,  is  the  permittivity  of  free  space  and  e  is  the 

electron  charge  magnitude  (Table  6.2).  Note  that  only  the  Coulombic  term  depends  on 
the  angular  momentum  L  whereas  both  the  Coulombic  and  the  cutoff  correction  terms 
depend  on  the  momentum,  hk  .  It  is  apparent  from  the  plot  of  these  two  terms  (Figure 
6.8)  that  the  cut-off  correction  has  little  impact  on  the  solution  and  that  the  Coulombic 
term  dominates.  Given  our  calculated  'S'g  pp  phase  shift  <7^ ,  a  solution  that  more  closely 


aligns  with  ANL’s  data  may  be  obtained  by  using  equation  (5.60), 

=SL-cr^+rH2kp) 

The  corrected  CPM  solution  is  also  shown  in  Figure  6.8. 


(6.8) 


Table  6.2.  *So  Coulombic  and  Cut-off  Correction  Parameters 


Parameter 

Quantity 

Units 

Angular  Momentum,  L 

0 

Electron  Charge  Magnitude,  e 

1.602x10'''' 

C 

Permittivity  of  Free  Space, 

8.854x10'^' 

\farad/fm] 

Coulombic  Cut-off  Point,/? 

140 

fm 

Proton  -Proton  Reduced  Mass,  Ppp 

0.01204830 
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Figure  6.8.  The  Coulombic  and  Cut-off  Correetions  (upper  diagram)  to  the  *SoPp  Phase  Shift  (lower 
diagram).  The  cutoff  correction  is  small  with  respect  to  the  magnitude  of  the  Coulomhic  correction. 
When  these  corrections  are  applied,  the  *So  CPM  result  aligns  with  the  published  results  from  ANL. 
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6.2.  Wave-Packets  and  the  Asymptotic  Limit 


We  now  look  at  the  non  2x2  isospin  triplet  for  eaeh  value  of  total  momentum,  J,  through 
J  =  3.  In  the  previous  section,  the  angular  momentum,  L,  was  zero  and  had  no  impact  on 
the  calculation.  However,  this  is  not  the  case  for  both  the  remaining  1x1  and  2x2  isospin 
triplets.  As  L  increases,  so  does  the  centrifugal  barrier  contribution  to  the  potential 
(Figure  6.9). 


Figure  6.9.  Centrifugal  Barrier  Potentials  for  i,=l-3.  As  the  centrifugal  barrier  increases,  the 
nuclear  potential  represents  only  a  minor  perturbation  to  the  barrier.  Phase  shifts  for  Z  >  2  are 
expected  to  be  small. 


In  fact  for  Z  >  2  ,  the  short  range  nuclear  well  represents  only  a  minor  perturbation  to  the 
barrier  potential. 

In  the  previous  'So  phase  calculations,  the  wave-packets  could  be  placed  relatively 

close  to  the  interaction  region  at  400  fin.  For  that  problem,  the  wave  packets  could  be 
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considered  in  the  asymptotie  limit  since  the  AVI 8  potential  terminated  at  and  the 

centrifugal  barrier  provided  no  contribution  (See  Figure  6.2).  Here  howeverZ  0 ,  so  a 
choice  for  the  initial  position  of  the  wave  paekets  must  be  made  further  away  from  the 
origin.  Ideally,  a  point  will  be  chosen  close  enough  to  the  interaction  region  to  minimize 
eomputational  effort  but  still  far  enough  away  from  the  area  to  minimize  the  barrier 
potentials  impact  on  the  correlation  function.  As  an  initial  guess,  points  are  chosen  where 
each  potential  falls  to  roughly  1x10'^ Me F.  For  the  ^Pi,  this  oeeurs  at  a  distanee  of  800 fm 
(Figure  6.10). 


800  1000  1200  1400  1600  1800  2000  2200  2400 
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Figure  6.10.  Decay  of  the  Barrier  Potential  at  Higher  Values  of  r.  The  asymptotic  limit  was  initially 
assumed  to  he  where  the  harrier  potential  fell  helow  1x10  “*  Me  F  Moller  states  were  established  to 
begin  wave  packet  propagation  based  on  this  assumption  at  800  fm,  1600  fm,  and  3200  fm  for 

L  =  \,  2,  and  3  respectively. 
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For  each  higher  order  L  calculation,  the  initial  positions  were  doubled  such  that 
the  ^T>2  and  cases  started  at  1600  fin  and  3200  fim  respectively.  To  accommodate 
these  changes,  the  coordinate  grid  size  and  number  sampling  points  were  increased  by 
factors  of  two  to  maintain  the  0.195  fim  coordinate  step  size.  The  same  principle  was 
applied  to  the  number  of  time  steps  and  .  The  values  for  Ax  ,  k,  and  packet  spread 

were  the  same  as  those  previously  given  by  Table  6.1.  The  results  of  the  pn  CPM 
calculations  are  presented  in  the  next  three  figures. 


Figure  6.11.  ^Pi  pn  Phase  Shift.  Calculated  from  Wave  Packets  Initially  Placed  at  800 fm. 
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Figure  6.12.  'Di  pn  Phase  Shift  Calculated  from  Wave  Packets  luitially  Placed  at  1600  fin. 


Figure  6.13.  ^Fj  pu  Phase  Shift  Calculated  from  Wave  Packets  luitially  Placed  at  3200  fm.  The  low 
euergy  uoise  in  this  plot  and  the  previous  two  plots  results  from  a  poor  choice  of  wave  packet.  The 
disagreement  with  the  ANL,  which  increases  as  a  function  of  L,  is  a  result  of  a  had  asymptotic  limit 
assumption. 
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Although  these  solutions  track  well  with  ANL’s  published  results,  two 
deficiencies  are  readily  apparent.  First,  there  is  a  high  degree  of  noise  in  the  low  energy 
regime  (below  50  MeV),  which  suggests  an  inappropriate  wave-packet  choice  was  made 
to  resolve  this  specific  energy  range  adequately.  Our  original  wave  packet  amplitude  was 
small  at  low  MeV  energies.  How  to  increase  resolution  at  low  MeV  energies  is  discussed 
in  Section  6.3.  The  other  more  noteworthy  deficiency  is  that  the  error  increases  with  each 
increase  with  L.  This  suggests  a  poor  choice  was  made  for  the  asymptotic  assumption.  In 
an  attempt  to  correct  this  problem,  the  initial  start  positions  and  grid  sizes  were  doubled 

-5 

yet  again.  The  F3  phase  shifts  from  this  trial  were  still  not  satisfactory.  Crude 
calculations  showed  (See  Figure  6.14)  that  for  the  F3  an  initial  position  of  12,800  fin,  a 

17  7^ 

coordinate  grid  of  2  points,  and  temporal  grid  of  2  points  would  be  more  suitable. 


Figure  6.14.  Crude  ^Fs+Phase  Shift  Calculated  from  Wave  Packets  Initially  Placed  at  12,800  fm.  The 
results  of  this  crude  calculation  indicate  that  placing  the  wave  packets  at  12,800  fm  from  the  nuclear 
well  would  he  a  good  asymptotic  limit  approximation. 
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Projections  for  both  MSRC  and  Desktop  processing  based  on  previous  optimized  run 
times  indicated  these  grid  sizes  would  require  a  minimum  14-15  days  per  calculation  for 
the  F3  and  high  L  cases  independent  of  platform  (Table  5.2).  Parallelization  provides 
little  or  no  additional  benefit  as  the  CPM  is  a  serial  algorithm. 


Table  6.3.  Projected  MSRC  and  PC  Run-Times 


Coord/Time  Points 

MSRC-HPCll  RunTime  (Days) 

P4  2400  Desktop  RunTime  (Days) 

17,17 

0.22 

0.24 

17,18 

0.44 

0.47 

17,19 

0.88 

0.94 

17,20 

1.75 

1.89 

17,21 

3.50 

3.78 

17,22 

7.00 

7.56 

17,23 

14.0 

15.11 

17,24 

28.0 

30.22 

Dispersion  further  compounds  the  problem.  As  we  move  farther  from  the  interaction 
region,  the  lower  energy  wave-packet  components  require  more  time  to  probe  the 
potential  and  return  to  the  initial  position;  making  the  computational  cost  of  this 
technique  become  less  attractive.  The  next  section  presents  an  alternative  methodology 
for  handling  these  two  issues. 

6.3.  An  Intermediate  State  Calculation 

As  we  saw  in  the  last  section,  the  computational  cost  of  placing  the  wave-packet 
in  the  asymptotic  limit  can  be  extremely  high.  The  following  alternative  approach 
leverages  the  specific  problem  characteristics  to  our  benefit.  This  section  is  divided  into 
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two  subsections.  The  first  deals  with  the  1x1  matrix  elements  of  and  the  second 
(Section  6.3.2)  deals  with  the  2x2  elements  matrix  elements  ofV^ . 


6.3.1.  The  Single  Element  Case 


To  minimize  dispersion  and  accelerate  subsequent  calculation  of  the  correlation 
function,  it  would  be  useful  to  place  the  wave-packets  as  close  as  possible  to  the 
interaction  region.  A  useful  choice  would  be  to  center  the  wave  packets  at  400 /m  as  we 
did  in  Section  6.1  since  it  is  well  outside  the  AVI  8  pp  cutoff  point  and  far  enough  away 
from  the  absorbing  boundaries  for  them  to  pose  any  problems.  In  order  to  start  the  wave- 
packet  calculations  in  closer  proximity  to  the  well,  we  return  to  the  definition  of  the 
Moller  operator 


=  lim  Exp 


iHt 

~Y 


Exp 


-iH^t 

h 


(6.9) 


Previous  calculations  assumed  that  the  propagation  began  in  the  asymptotic  limit  solely 
under  the  full  Hamiltonian,  H.  Since  acceptable  results  were  not  obtained  under  this 
assumption,  the  full  form  must  be  employed. 

If  we  recall  that  the  first  operator.  Ho,  contains  only  kinetic  energy  information, 
we  can  use  equation  (5.82)  to  analytically  determine  the  reactant  (product)  wave  packet’s 
form  at  some  position  where  the  centripetal  barrier  is  negligible.  Here,  r  was  increased 
until  the  wave  packet  peaks  were  well  outside  our  asymptotic  estimate  of  12,800 fm 
determined  in  section  5.2.  (Figure  6.15)  These  intermediate  states  were  then  propagated 
back  to  r  =  0  under  the  full  Hamiltonian,  H,  using  the  propagation  algorithm.  This 
process  established  new  Moller  states  for  the  phase  shift  calculations. 
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Since  it  requires  eonsiderable  time  to  propagate  the  wave-paekets  baek  from 
-r  ( r ),  it  would  be  useful  to  know  if  the  time  step  eould  be  inereased  given  that  over  this 
region  the  potential  remains  negligible.  So,  eomparisons  were  made  between  a  wave 
paeket  determined  analytieally  at  235  and  a  wave  packet  propagated  to  that  same 

point  using  the  propagation  algorithm  using  Ar ’s  of  0.00022  ,  0.00045  z^ ,  0.00090  z^ , 

and  0.0018  .  Absolute  error  tests  indicated  the  time  step  eould  be  inereased  by  a  factor 

of  2  without  significant  waveform  degradation. 


Figure  6.15.  Absolute  Value  of  Dispersed  luitial  Reactaut  Wave-Packet.  The  dispersed  wave  packet 
was  calculated  usiug  the  aualytic  free  space  propagator  equatiou  (5.82).  The  peak  of  the  wave  packet 
was  ceutered  at  twice  the  distauce  iudicated  by  the  crude  calculatiou  showu  iu  Figure  6.14  to  eusure 
mioimal  barrier  iuterfereuce.  This  dispersed  wave  packet  was  theu  propagated  back  to  the  luitial 
positiou  to  iutegrate  barrier  ioformatiou  iuto  the  Mailer  state  for  scatteriug  matrix  calculatious. 
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There  is  an  inherent  simplieity  to  this  approach.  Consider  the  kinetic  energy 


tumbling  term, 


L{L  +  \) 


(6.10) 


Notice  that  the  only  difference  between  potentials  outside  the  interaction  region  for  a 
given  choice  of  wave-packet  is  due  to  the  reduced  mass,  /u  ,  and  angular  momentum,  L. 
This  implies  that  if  the  reactant  and  product  Moller  states  differ  only  by  the  sign  of  ko,  the 
intermediate  Moller  state  need  only  be  calculated  once  for  a  specific  L  and  ju 
combination.  In  other  words,  these  Moller  states  can  stored  and  used  repeatedly.  This 
capability  facilitates  rapid  calculation  of  improved  phase  shifts  as  both  our  understanding 
of  the  strong  force  and  our  models  that  describe  this  understanding  improve. 

The  usefulness  of  this  approach  is  apparent  if  we  consider  the  calculation  of  the 
neutron-proton  L=  3  phase  shifts,  3F2,  1F3,  3F3  and  3F4.  From  Table  6.3,  it  would 
require  at  least  28  days  of  computational  time  to  calculate  each  phase  shift  on  a  grid  the 
size  shown  in  Figure  6.15  (1 12  days  total).  However  if  the  Moller  state  is  pre-calculated 
as  described  in  the  preceding  paragraphs,  the  same  four  calculations  can  be  performed  in 
less  than  12  hours.  This  is  a  99%  reduction  in  computer  runtime.  Subsequent  neutron- 
proton  L  =  3  phase  shifts  can  be  obtained  from  these  intermediate  Moller  state  in  a  matter 
of  minutes.  This  capability  provides  the  nuclear  community  with  a  powerful  new  tool  to 
calculate  detailed  phase  shift  information  rapidly  as  new  improved  nuclear  models 
become  available.  For  these  and  all  remaining  phase  calculations,  three  wave-packets 
were  employed  to  reduce  the  distortion  prevalent  below  50  MeV  in  the  Section  5.2 
outputs  (Figure  6.16  and  Figure  6.17) 
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Figure  6.16.  Initial  Positions  of  the  Three  Wave-Packets 
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To  achieve  low  energy  resolution,  the  momentum  wave  packet  must  be  narrow  in 
k  spaee  and  elose  to  zero  as  shown  in  Figure  6.17.  However  the  narrower  in  k  spaee  we 
make  the  wave  paeket  the  broader  the  eoordinate  wave  paeket  beeomes  (Figure  6.16). 
This  neeessitates  initial  eoordinate  wave  paeket  positions  farther  from  the  interaetion 
region.  The  eoordinate  initial  peak  positions  used  here  are  eentered  at  200 fm,  300  fm, 
and  400 fm  whieh  eorrespond  to  k  spaee  peaks  eentered  ±  1.253  [Hfin],  ±  0.450  [Hfm], 
and  ±0.150  [l//m]  respeetively.  The  parameters  used  to  ealeulate  the  Moller  states  from 
these  wave  paekets  are  given  in  Table  6.4. 


Table  6.4.  Parameters  used  to  calculate  the  Moller  states  for  each  wave-packet 


Broad  Energy  Wave  Packet 

^max 

51200 

fm 

Time 

235 

Coordinate/Temporal  Points 

2‘®,2'’ 

Coordinate  Packet  Peak 

200 

fm 

Momentum  Packet  Peak 

-  1.253 

1/  fm 

Packet  Spread 

1.814 

fm 

Mid  Energy  Wave  Packet 

^max 

102400 

fm 

Time 

1320 

Coordinate/Temporal  Points 

2*^2'^ 

Coordinate  Packet  Peak 

300 

fm 

Momentum  Packet  Peak 

-  0.450 

1/  fm 

Packet  Spread 

5.277 

fm 

Low  Energy  Wave  Packet 

^max 

102400 

fm 

Time 

3950 

Coordinate/Temporal  Points 

2*^2'^ 

Coordinate  Packet  Peak 

400 

fm 

Momentum  Packet  Peak 

-  0.150 

1/  fm 

Packet  Spread 

16.491 

fm 
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For  each  wave-packet,  an  intermediate  state  was  calculated  using  equation  (5.82) 
under  the  assumption  that  a  dispersed  peak  location  of  at  least  25,000  fin  from  the 
interaction  region  would  be  an  acceptable  asymptotic  limit  approximation  (Figure  6.18). 
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Figure  6.18.  Absolute  Value  of  the  Three  Wave  Functions  in  the  Intermediate  Reactant  States  in  the 
Asymptotic  Limit.  The  peak  positions  of  the  lower  energy  wave  packets  were  doubled  to  further 
reduce  the  possibility  of  barrier  interference. 

As  an  extra  precaution,  the  intermediate  peak  positions  of  the  low  energy  wave  packets 
were  doubled  to  further  reduce  the  possibility  of  barrier  interference.  These  asymptotic 
states  were  then  propagated  back  toward  the  interaction  region  under  their  respective  full 
Hamiltonians  until  r  =  0  using  our  propagation  algorithm;  reconstituting  our  non- 
dispersed  wave  packet  with  the  desired  centripetal  barrier  information.  Phase  shifts  were 
then  calculated  with  the  three  wave-packets  as  before  but  with  different  propagation 
times  since  lower  k  wave  packets  are  slower  and  require  more  time  to  interact  with  the 
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potential  and  return  to  their  initial  positions.  The  results  of  these  ealeulations  were  then 
splieed  together  to  present  a  eomplete  phase  shift  picture.  The  parameters  used  to 
calculate  the  phase  shifts  using  these  Moller  states  are  provided  are  Table  6.5.  On  the 
next  series  of  pages,  the  Pi,  Di,  and  F3  pn  and  pp  phase  shifts  are  shown.  From  this 
series  of  plots,  it  is  apparent  that  very  good  agreement  exists  between  the  ANL  values 
and  our  CPM  technique  down  to  a  few  MeV.  Additional  phase  shifts  combinations  are 
included  in  Appendix  1 . 

Table  6.5.  Parameters  Used  to  Calculate  the  Phase  Shifts  from  Three  Moller  States. 


Broad  Energy  Wave  Packet 

Xmax 

1600 

fin 

Time 

80 

A 

Coordinate/Temporal  Points 

Coordinate  Packet  Peak 

200 

fin 

Momentum  Packet  Peak 

-/+  1.253 

1/ fin 

Packet  Spread 

1.814 

fin 

Mid  Energy  Wave  Packet 

Xmax 

1600 

fin 

Time 

160 

A 

Coordinate/Temporal  Points 

2'^2'^ 

Coordinate  Packet  Peak 

300 

fin 

Momentum  Packet  Peak 

-/+  0.450 

1/ fin 

Packet  Spread 

5.277 

fin 

Low  Energy  Wave  Packet 

Xmax 

1600 

fm 

Time 

320 

A 

Coordinate/Temporal  Points 

2'^2^'’ 

Coordinate  Packet  Peak 

400 

fin 

Momentum  Packet  Peak 

-/+  0.150 

1/ fin 

Packet  Spread 

16.491 

fin 
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Figure  6.19.  Composite  CPM  Calculation  for  the  ^Pi  pn  and  the  ^Pi  pp  Phase  Shift. 
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Figure  6.20.  Composite  CPM  Calculation  for  the  pn  and  the  *D2  pp  Phase  Shift 
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Figure  6.21.  Composite  CPM  Calculation  for  the  pn  and  the  ^Fj  pp  Phase  Shift 
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6.3.2.  The  2x2  calculation 


The  first  2x2  block  of  (Si  and  Di)  represents  our  most  difficult  application 

of  the  CPM  thus  far.  S-Matrix  element  calculations  for  the  2x2  involve  propagation  on 
two  coupled  surfaces  and  subsequent  computation  of  four-correlation  functions.  Since 
the  potential  in  the  \L,S,J,T'^  basis  is  not  diagonal,  it  is  more  efficient  to  transform  to  a 

diagonal  basis  to  simplify  exponentiation.  However,  the  kinetic  energy  matrix  is 
undiagonalized  as  a  consequence  of  this  transformation.  Therefore,  propagation  requires 
two  extra  steps  to  transform  between  representations  to  the  basis  where  the  respective 
matrix  is  diagonal  before  operating  on  the  wave  vector.  These  unitary  transformations 
were  identified  by  U  and  back  in  equation  (5.29). 

Although  the  same  three  initial  wave-packets  from  the  last  section  will  be  used  for 
this  calculation,  to  help  illustrate  the  process  here  we  only  discuss  the  wave  packet 
starting  at  200 fm.  A  sequence  of  four  frames  (Figure  6.22-Figure  6.25)  from  the  multi¬ 
surface  propagation  is  provided  to  guide  the  discussion.  In  the  first  frame,  the  reactant 
wave-packet  is  started  on  the  lower  Si  surface  and  two  identical  product  wave-packets 
are  positioned  at  the  same  location  to  determine  the  scattering  probability  to  each 
channel.  The  lower  product  state  measures  the  reflected  amplitude  and  phase  in  the  Si 
channel  and  the  upper  wave-packet  measures  the  magnitude  and  phase  that  is  transferred 
to  the  coupled  channel.  By  the  third  frame,  the  incident  wave-packet  has  entered  the 
interaction  region  of  the  potential  and  has  become  mixed,  as  amplitude  is  apparent  of  the 
upper  surface.  Both  waveforms  have  started  leaving  the  interaction  region  by  the  fourth 
frame.  The  propagation  is  allowed  to  continue  until  the  correlation  function  between  the 
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incident  and  product  wave-packets  on  both  ehannels  drops  to  the  order  of  10'^, 
completing  the  first  part  of  the  calculation.  The  same  eomputation  is  then  eomputed 
again  except  the  wave-packet  is  started  on  the  upper  Di  surface  to  complete  the  2x2 
bloek  (not  shown). 


Figure  6.22.  Initial  2x2  Wave-Packet  Positions  (Scaled)..  The  reactant  wave  packet  is  superimposed 
over  the  product  wave  packet  on  the  ^Si  channel.  In  this  diagram,  the  ^Di  surface  is  not  shown  and 
the  first  off  diagonal  coupling  surface  from  Figure  5.2  is  artificially  displaced  for  clarity  along  with 
the  ^Di  product  Mailer  state.  Absolute  values  of  the  wave  packet  are  plotted. 
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Figure  6.23.  Evolving  Packet  Converges  on  Potential  (Scaled) 
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Figure  6.25.  Mixed  Packets  Leave  Interaction  Region  On  Both  Channels  (Scaled).  The  lower  energy 
components  of  the  wave  packet  have  yet  to  leave  the  interaction  region. 

The  phase  outputs  obtained  from  this  2x2  calculation  are  highlighted  Figure  6.26 
and  Figure  6.27.  Again,  good  agreement  exists  between  the  CPM  and  the  data  provided 
by  ANL.  The  only  slight  abnormality  is  in  the  low  energy  region  of  the  which  is  on 
the  order  of  a  tenth  of  a  degree  and  considered  to  be  negligible.  Additional  computation 
can  reduce  this  error  further  but  was  deemed  not  necessary  since  it  has  already  been 
shown  that  enhanced  low  energy  results  can  be  obtained  with  the  selection  of  an 
appropriate  wave  packet.  Here,  L  has  even  parity  {L=0  or  2)  and  S  is  symmetric  ( S'  =  1 ). 
So,  the  wave  function  is  symmetric  and  the  pp  and  nn  phase  shifts  do  not  need  to  be 
computed.  However,  in  the  next  2x2  block,  ^  and^Fj ,  all  three  will  need  to  be 
determined  since  L  and  S  are  anti-symmetric. 
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Figure  6.26.  Composite  CPM  Calculation  for  the  ^Si  pn  (top  frame)  and  the  s  hari  pn  (bottom 
frame)  Phase  Shifts.  The  e  hari  pn  phase  shift  is  associated  with  the  off  diagonal  component  of  the 
first  2x2  block  in  the  potential  matrix  given  in  Figure  5.2 
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Figure  6.27.  Composite  CPM  Calculation  for  ^Di  pn  phase  shift  (top  frame).  As  the  lower  frame 
illustrates,  the  three  wave  packets  chosen  to  perform  these  calculations  provide  good  agreement  with 
the  published  results  down  to  about  1  MeV. 
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For  the  ^  and^F2  2x2  block,  phase  shifts  were  computed  for  the  isospin  triplet. 
These  are  shown  in  Figure  6.28  through  Figure  6.3 1 .  Of  the  three  possible  combinations 
only  the  pn  and  pp  are  given  here,  the  remaining  nn  phase  shifts  are  included  in  the 
Appendices. 


Figure  6.28.  Composite  CPM  Calculation  for  the  ^P2  pn  Phase  Shift . 
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Figure  6.29.  Composite  CPM  Calculation  for  the  E_bar2  pn  Phase  Shift  (upper  frame)  and  the  ^F2 
pn  Phase  Shift  (lower  frame)..  The  slight  difference  in  phase  at  about  40  MeV is  a  splicing  artifact . 
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Figure  6.30.  Composite  CPM  Calculation  for  the  ^PiPP  Phase  Shift  (upper  frame)  and  the  e_har2 
pp  Phase  Shift  (lower  frame). 
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Figure  6.31.  Composite  CPM  Calculation  for  the  ^F2  pp  Phase  Shift.  The  deviation  from  the 
published  results  stems  from  how  the  Coulomhic  correction  is  applied  to  the  CPM  results.  Higher  L 
calculations  of  either  1x1  pp  or  2x2  pn  phase  shifts  do  not  exhibit  this  deviation 

-5  -5 

Here,  there  is  again  good  agreement  with  the  published  results  for  the  P2,  8_bar2,  and  F2 

•V 

pn  phase  shifts.  There  is  only  a  small  (tenth  of  a  degree)  step  around  40  MeV  in  the  F2 
pn  phase  shift  whieh  is  merely  a  splieing  artifact.  The  same  step  is  evident  in  the  P2 pn 
phase  shift  when  placed  on  the  same  scale.  What  immediately  stands  out  about  the 
complementary  2x2  pp  phase  shifts  is  the  deviation  of  our  F2  phase  shift  from  the 
published  results  (Figure  6.31).  This  same  error  is  present  in  the  P2  pp  phase  shift,  it  is 
just  not  as  apparent  because  the  scale  of  the  plot  is  an  order  of  magnitude  larger.  Since 
the  error  is  not  apparent  in  the  P2,  s_bar2,  and  F2  pn  shifts  or  1x1  pp  phase  shifts  of 
higher  order  L  (Figure  6.21),  the  source  of  the  error  must  be  associated  with  how  we 
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applied  the  Coulombie  eorreetion  to  the  8_bar2,  and  ^F2  2x2  pp  bloek.  Stapp  [18]  has 
a  possible  explanation. 

The  Coulombie  eorreetion  developed  in  Seetion  5.6.4  eorresponds  to  what  Stapp 
[18]  ealls  a  “pure”  state  ealeulation  or  a  1x1  pp  matrix  element  of  .  In  the  2x2  pp 

ealeulation,  the  evolving  wave  paekets  enter  the  interaetion  region  and  are  mixed  yielding 
s-matrix  elements  and  phase  shifts  that  are  linear  eombinations  of  P  and  F.  These  final  P 
and  F  states  are  no  longer  “pure.”  Ideally,  we  eould  determine  the  unitary  matrix  whieh 
would  diagonalize  the  P2,  8_bar2,  and  F2 pp  2x2  bloek  and  use  the  transformation  to 
eompute  a  mixed  basis  Coulombie  eorreetion.  Here  we  ehoose  however  to  follow 
Stapp’s  adviee.  In  his  1957  paper,  he  states  that  the  “pure”  Coulombie  eorreetion  may 
be  applied  direetly  to  the  “bar”  phase  shifts  without  introdueing  a  signifieant  amount  of 
error.  The  tenth  of  a  degree  phase  differenee  between  the  published  results  and  our  F2 
pp  CPM  solution  is  not  signifieant  and  as  Figure  6.32  illustrates  well  within  the  range  of 

-5 

F2  results  obtained  from  other  models  [23-25,  32]. 
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Figure  6.32.  Comparison  of  the  Calculated  CPM  ^F2  Phase  Shift  to  Published  Results  from  Other 
Time-Independent  Calculations.  The  tenth  of  a  degree  difference  between  the  CPM  result  and  the 
ANL  results  is  not  significant  when  compared  to  the  differences  between  the  AV18  results  and  the 
other  models  results. 
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6.3.3.  The  Cross  Section 

Nuclear  calculations  rely  heavily  on  the  determination  of  accurate  cross  sections 
to  ascertain  the  probability  of  reaction  between  nuclides.  The  phase  shifts  calculated  in 
the  previous  section  lend  themselves  readily  to  the  determination  of  the  nn,  pn,  and  pp 
cross-sections.  Cross  sections  may  be  calculated  by  taking  a  weighted  sum  of  the  triplet 
and  singlet  phase  shifts  for  a  specified  value  of  angular  momentum  L  given  by  [7,  78], 

(3  1  ^ 

^  +-«singie,"  J  •  (6-11) 

with, 

1 

4„.,=TYl;(2i  +  l)sm^<5, 

(6.12) 

1  ^  1 

^triplet  +  1) (sin'  +  sin"  +  sin"  ) 

where  k  is  the  wave  vector.  In  these  calculations,  partial  waves  above  L  =  2  were  not 
considered  in  the  summation  since  sin^  (angle  <  1  degree)  is  very  small  and  does  not 
contribute  significantly  to  the  total  cross  section.  The  calculated  pn,  pp,  and  nn  total 
cross  sections  in  units  of  bam  (1  barn  =  10  )  are  presented  in  Figure  6.33  and  Figure 

6.34  respectively. 

It  was  our  original  intention  to  compare  these  cross  sections  to  the  Evaluated 
Nuclear  Data  Files  (ENDF)  posted  on  the  Eos  Alamos  National  Eab  website.  Although 
there  was  ample  pn  cross  section  data  on  the  website,  complementary  pp  and  nn  data  was 
not  available.  Attempts  were  also  made  to  locate  pp  and  nn  total  cross  sections  in  the 
literature  without  success.  After  an  extensive  search,  we  were  only  able  to  obtain 
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differential  cross  section  data  in  our  0-400  MeV  energy  range.  Nonetheless,  we  can  still 
make  a  couple  of  qualitative  assessments  about  these  plots.  One,  the  total  cross  section 
was  determined  correctly  from  the  phase  shifts  otherwise;  the  pn  phase  shift  in  Figure 
6.33  would  not  have  such  a  tight  correlation  with  the  published  ENDF  files.  Two,  the pp 
total  cross  section  (Figure  6.34)  is  the  total  cross  section  due  to  short-range  forces  alone. 
This  is  not  what  would  be  observed  in  the  laboratory  as  the  long  range  Coulomb  force 
between  the  protons  would  combine  with  the  nuclear  potential  and  produce  a  different 
result.  Indeed,  Blatt  [7]  mentions  that  below  \0  MeV,  the  only  partial  wave  appreciably 
changed  by  the  nuclear  well  is  phase  shift.  All  other  partial  waves  with  T  >  0  would 

be  scattered  by  the  Coulomb  force  long  before  the  particles  could  become  close  enough 
to  experience  nuclear  effects. 


Figure  6.33.  Total  up  Scattering  Cross  Section  Compared  to  the  LANL  ENDF  Data. 
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Figure  6.34.  Total  pp  Scattering  Cross  Section  (top  frame)  and  the  Total  nn  Scattering  Cross  Section 
(bottom  frame).  It  is  important  to  remember  that  the  pp  cross  section  presented  here  is  the  total 
cross  section  due  to  short  range  forces  alone.  No  published  results  could  be  located  for  comparison 
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7.  Conclusion  and  Recommendations  for  Future  Work 


This  research  effort  is  the  first  to  demonstrate  complete  determination  of  nuclear 
scattering  matrix  elements  for  proton-proton,  neutron-proton,  and  neutron-neutron 
scattering  events  utilizing  the  time-dependent  CPM  technique.  Wave-packet  interaction 
with  a  potential  in  this  technique  is  analogous  to  a  laboratory  experiment  and  provides 
intuitive  insight  into  scattering  dynamics.  This  is  a  distinct  advantage  over  the  current 
time-independent  scattering  methods  used  by  the  nuclear  community,  which  retain  no 
additional  collision  data. 

This  work  was  also  the  first  to  demonstrate  how  Moller  states  can  be  leveraged  to 
rapidly  calculate  high-energy  resolution  scattering  matrix  elements  and  phase  shifts  for 
an  entire  block  of  angular  momentum,  J.  This  provides  an  ancillary  utility.  Since  the 
nuclear  potential  has  such  a  short  range,  these  same  Moller  states  can  be  stored  as  reused 
to  compute  improved  phase  shifts  as  higher  precision  nuclear  models  become  available 
with  minimal  computational  effort. 

The  CPM  determines  the  entire  phase  shift  for  a  range  of  energies,  whereas  time- 
independent  techniques  typically  obtain  a  range  of  phase  values  for  a  specified  energy. 
Thus  depending  on  the  desired  information,  CPM  may  be  more  advantageous.  CPM  may 
also  be  advantageous  as  the  problems  become  more  complex  since  computational  cost  of 
this  matrix  multiplication  technique  scales  as  N  or  better  versus  the  time-independent 
approaches  which  scale  as  N  [50].  This  section  details  the  conclusions  and  outlines  the 
recommendations  for  future  research  initiatives. 
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When  we  initially  set  up  our  ehannel  paeket  ealeulations,  we  attempted  to  plaee 
the  Moller  states  within  few  thousand  femtometers  of  the  nuelear  well  and  ealeulate  the 
nuelear  seattering  matrix  elements  under  the  assumption  that  the  eentripetal  barrier’s 
amplitude  was  small  enough  as  to  not  interfere  with  the  ealeulation.  This  was  an 
ineorreet  assumption.  Eaeh  inerease  in  angular  momentum  eorresponded  to  a  signifieant 
inerease  in  wave-paeket  start  position.  From  test  ealeulations  for  F  =  5  and  4,  it  was 
shown  that  the  eomputational  eost  of  positioning  the  wave-paekets  in  a  good  asymptotie 
limit  approximation  was  not  praotieal.  As  many  as  14-15  days  would  be  required  for 
eaeh  phase  shift  ealeulation  for  Z  >  3  to  aehieve  suitable  eonvergenee  to  the  published 
data. 

In  an  alternative  approaeh,  an  analytie  equation  was  used  to  obtain  an 
intermediate  state  in  a  suitable  asymptotie  limit  approximation.  Then,  the  time  dependent 
algorithm  was  used  to  ineorporate  information  about  the  eentripetal  barrier  into  eaeh 
wave  paeket  as  the  wave  paekets  were  propagated  baek  to  their  original  positions  elose  to 
the  nuelear  well.  It  was  shown  that  not  only  eould  these  Moller  states  be  ereated  at  a  time 
step  2  greater  than  utilized  in  the  original  approaeh  but  that  these  Moller  states  need  only 
be  ealeulated  onee  for  eaeh  nueleon-nueleon  and  angular  momentum  pairing.  This  shift 
in  approaeh  resulted  in  a  dramatie  25-50  fold  inerease  in  eomputation  effieieney  and 
reprodueed  published  results  from  ANL  with  added  energy  resolution. 

The  work  presented  in  this  dissertation  provides  a  foundation  for  future  time- 
dependent  nuelear  seattering  ealeulations.  There  are  two  primary  proposals  for  future 
work.  One  eenters  on  applieation  of  this  teehnique  to  a  three  nueleon  seattering  problem. 
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The  other  focuses  on  the  development  of  an  analytic  technique  for  propagating  a  wave- 
packet  on  the  centrifugal  barrier.  These  two  research  initiatives  are  discussed  below. 

7.1.  The  Three-Body  Problem 

The  three-body  problem  is  the  next  logical  extension  of  this  research.  Here,  there 
are  three  interacting  partners  that  require  nine  coordinates  to  completely  describe  the 
motion.  Fortunately,  nature  helps  limit  the  sheer  number  of  possible  scattering  partners 
to  two  since  the  deuteron  is  the  only  two-nucleon  system  to  exhibit  binding.  For  that 
reason,  an  incident  proton  and  an  incident  neutron  on  a  deuteron  target  are  the  only  two 
possible  three-body  combinations.  Here,  as  in  the  two-body  problem,  the  coordinate 
system,  Hamiltonian,  and  basis  must  be  determined  before  any  three-body  CPM 
calculations  can  be  performed. 

Although  the  choice  could  be  made  to  perform  this  calculation  in  the  Space-Fixed 
Center  of  Mass  frame  as  we  did  in  the  two  nucleon  scattering  problem,  a  more  suitable 
coordinate  system  for  the  three-nucleon  problem  is  the  body  fixed  coordinate  system.  In 
this  coordinate  system,  the  motion  of  the  center  of  mass  is  eliminated  by  selecting  Jacobi 
coordinates  to  describe  the  interaction.  These  are  shown  in  Figure  7.1  . 
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Figure  7.1.  Jacobi  Coordinates.  The  Jacobi  coordinate  system  is  a  useful  choice  for  performing 
three-body  time-dependent  scattering  calculations. 


Here,  R  is  aligned  along  the  body  fixed  z-axis  and  represents  the  separation  between  the 
deuteron  and  the  ineident  particle,  r  is  the  distance  between  the  deuteron  constituents,  (f) 

is  the  azimuthal  angle  between  the  R  and  r  ^  and  0  is  orientation  of  the  di-nucleon  about 
body  z  axis. 

The  Full  Hamiltonian  well  suited  for  three-body  problem  is. 


H  = 


-  +  - 


'^R-a,bc 


+v 

2  Interaction 


liu 


P 

^  +V 

A,BC 


(7.1) 
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where  refers  to  the  tumbling  of  the  eomposite  system  around  the  eenter  of  mass  and 
refers  to  the  rotation  of  the  di-nueleon.  These  terms  are  ineorporated  into  the  potential 
ereating  an  effeetive  potential  of  the  form, 


(7.2) 


where  the  AVI 8  model  represented  here  as  v.j  is  augmented  by  a  three  body  eorrection, 

.  Neither  the  AVI  8  Model,  nor  any  of  the  other  NN  models,  serves  as  a  suitable 

deseription  of  the  three-body  potential  by  themselves  sinee  all  fail  to  adequately  deseribe 
the  eoupling  between  three  interaeting  partners.  Bound  state  energies  of  the  Tritium  atom 
in  the  AVI 8  are  underestimated  by  around  by  1  MeV eompared  to  the  measured  value  of 
8.48  MeV  [8].  As  shown  in  Figure  7.2,  the  agreement  between  the  AV18  binding 
predietion  (blue  bars)  and  observation  (green  bars)  eontinues  to  deteriorate  rapidly  with 
inereasing  N  [80].  To  eompensate  for  the  underbinding  in  heavier  nuelei,  ANL  and  the 
University  of  Illinois  at  Urbana  have  developed  a  new  three-body  eorreetion  designated 
as  the  IL2.  The  IL2  (red  bars  in  Figure  7.2)  has  been  shown  to  be  able  closely 
approximate  observed  binding  energies  of  heavier  nuclei  through  ^4 =72.  At  this  time, 
however,  distributable  code  for  the  IL2  is  not  available  so  applying  the  CPM  technique  to 
the  three-body  problem  will  require  development  of  an  AFIT  IL2  FORTRAN  code  [79] 
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Figure  7.2.  AV18,  AV18/IL2,  and  Experimental  Binding  Energy.  [8,  80]  Note  the  improvement  in 
the  binding  energy  when  the  IL2  correction  is  applied  to  the  AV18.  The  IL2  correction  to  the  AV18 
has  been  able  to  reproduce  the  binding  energies  of  nuclei  containing  up  to  twelve  nucleons. 


A  well-suited  basis  set  in  which  to  represent  the  Hamiltonian  is  always  important. 
Since  the  core  of  the  three-body  potential  is  the  AVI  8,  we  expect  that  the  basis  set  would 
resemble  the  basis  set  utilized  in  the  two-body  analysis, 

\\l/^  =  \j  L  S  Sy  $2  T  R  (7.3) 

and  that  the  same  spin  (isospin)  coupling  encountered  in  Chapter  5.4  should  help  reduce 
the  basis  set  to  the  more  manageable  form, 
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\y/)  =  \j  LS  T  Rr)  (7.4) 

Here,  it  is  reeommended  that  the  researeh  be  broken  into  at  least  two  separate  projeets. 
The  first  should  foeus  on  developing  a  working  IL2  eode  and  ealeulating  phase  shifts  for 
the  elastie  problem  below  2.42  MeV.  The  seeond  projeet  should  foeus  on  the  inelastie 
seattering  problem  above  2.42  MeV.  For  this  seeond  problem,  the  CPM  teehnique  will 
need  to  be  extended  for  a  ehange  of  ehannel  from  A+  (BC)  to  A+B+C  before  inelastic 
phase  shifts  can  be  computed.  Both  projects  will  be  challenging. 


7.2.  Analytic  Determination  of  the  Moller  State 

Back  in  Section  6.3,  we  used  equation  (5.82)  to  analytically  determine  an 
intermediate  state  in  a  suitable  asymptotic  limit  approximation  then  used  the  time 
dependent  algorithm  to  propagate  the  wave  packets  back  to  their  original  positions  to 
establish  Moller  states.  This  approach  incorporated  centripetal  barrier  information  into 
each  wave  packet  prior  to  computing  scattering  matrix  elements  and  reduced  the  amount 
of  time  required  compute  scattering  matrix  elements  99%.  Although  we  were  able  to 
make  huge  gains  in  computational  efficiency,  it  still  took  almost  twelve  hours  to  compute 
the  initial  Moller  states  for  the  scattering  matrix  calculations. 

An  even  more  efficient  and  useful  approach  would  be  to  develop  an  analytic 
solution  to  the  time-dependent  radial  form  of  Schrodinger’s  equation. 


dt 


5"  L»{L  +  \)h^  ,,,  , 
-r  +  — ^ —  +  F(r) 


2//r  dr 


2/ur 


W 


(7.5) 


where  the  potential  V{r)  is  set  to  zero.  This  equation  has  the  inhomogeneous  PDF  form 


u,=ku^+f{x,t) 
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(7.6) 


and  can  be  expressed  in  terms  of  a  Green’s  function,  K  [81] 


y/{r,t)  =  j  K{r\r,t)y/{r\  Q)dr' 


which  relates  the  initial  eondition,  ’,0)  ,  to  some  final  state  at  some  later  time,  . 


Andrews  [44]  has  shown  the  solution  to  this  Green’s  funetion  kernel  for  the  repulsive  —r 


potential  as. 


[  n  J  t  2nt  ’ 


where  //  denotes  the  redueed  mass,  denotes  Bessel  funetions  of  the  first  kind,  and  the 
orbital  momentum  dependenee  is  eontained  within  dimensionless  parameter  v  , 


V-  —  +  L»(L  +  V) 

vu 


Sinee  solutions  to  v  are  in  n  +  —  powers  of  angular  momentum,  the  Bessel  functions  can 
be  expanded  in  terms  of  eosines  and  sines  [76] 


i^)  =  \  —  sin(x) 

^  \  TlX 


(7.10) 


j^(x)=  —  ^^^^-eos(x) 

:::  \  TTX  X 


(7.11) 


2(3  )  3 

J^{x)  =  .  — <  — -1  sin(x) — ^eos(x) 

o  \  71 X  \  \  X  J  X 


(7.12) 


Given  the  above  information,  a  test  ealeulation  was  performed  using  the  broad  energy 
wave  packet  parameters  given  in  Table  6.4  and  a  trapezoid  integration  algorithm. 


Although  similar  results  were  obtained  to  those  obtained  via  our  CPM  ealeulations,  an 
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additional  12  hours  was  required  to  complete  the  calculation.  Andrew’s  Green’s  function 
solution  may  be  a  good  starting  point  for  subsequent  research. 
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Appendix 


NN  PHASE  SHIFTS 
FOR/=<?-5 


Data  from  Nijmegen  (PWA)  experimental  data  used  when  ANL  data  not  available 
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NP  PHASE  SHIFTS 
FOR/=<?-5 


Data  Compared  To  Results  Provided  From  ANL 
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PP  PHASE  SHIFTS 
FOR/=<?-5 


Data  from  Nijmegen  (PWA)  experimental  data  used  when  ANL  data  not  available 


175 


Phase  Shift  (deg)  Phase  Shift  (deg) 


\ 


176 


Phase  Shift  (deg)  Phase  Shift  (deg) 


Q  /  I  I  I  I  I  I 

0  50  100  150  200  250  300  350 

MeV 


177 


Phase  Shift  (deg)  Phase  Shift  (deg) 


MeV 


178 


Phase  Shift  (deg)  Phase  Shift  (deg) 


MeV 


179 


Phase  Shift  (deg)  Phase  Shift  (deg) 


180 


Phase  Shift  (deg)  Phase  Shift  (deg) 


181 


182 


Bibliography 


1 .  Bethe,  Hans  A.,  “What  Holds  the  Nucleus  Together?”  Scientific  American,  V 1 89, 
58-63 (1953) 

2.  Machleidt,  R.  and  I.  Slaus.  “Topical  Review:  The  Nucleon-Nucleon  Interaction,” 
Journal  of  Physics  G,  27:  R69-R108  (2001) 

3.  Machleidt,  R.  “The  Theory  of  Nuclear  Forces:  Is  the  Never-Ending  Story  Coming 
to  an  End?,”  Nuclear  Physics  A,  790:  17c-23c  (2007) 

4.  Wiringa,  R.B.,  V.GJ.Stoks,  and  R.  Schiavilla.  “Accurate  Nucleon-Nucleon 
Potential  with  Charge-Independent  Breaking,”  Physical  Review  C,  V5I,NI:38-5I 
(January  1995) 

5.  Wiringa,  R.B.,  R.A.  Smith,  and  T.E.  Ainsworth.  “Nucleon-Nucleon  Potentials  With 
and  Without  A  (1232)  Degrees  of  Freedom  f  Physical  Review  C,  V29,N4:1207- 
1221  (April  1984) 

6.  Schiavilla,  R.,  R.B  Wiringa.,  V.G.J.Stoks,  W.  Glockle,  H.  Kamada,  A.  Nogga,  J. 
Carlson,  R.  Machleidt,  V.R.  Pandharipande,  and  A.  Kievsky,  “Weak  Capture  Of 
Protons  by  Protons,”  Physical  Review  C,  V58,N2: 1263-1277  (August  1998) 

7.  Blatt,  John  M.  and  V.  E.  Weisskopf  Theoretical  Nuclear  Physics.  New  York: 
Dover  Publications,  Inc.,  1991.  pp.  48-190 

8.  Pieper,  Steven  C.  V.R.  Pandharipande,  R.B  Wiringa,  and  J  Carlson.  “Realistic 
Models  of  Pion-Exchange  Three  Nucleon  Reactions,”  Physical  Review  C,  V64:l- 
21 (2001) 

9.  Eisenberg,  Judah  M.  and  W.  Greiner.  Microscopic  Theory  of  the  Nucleus  (Nuclear 
Theory  Volume  3).  ElservierPub  Co,  1976.  pp.  13-90. 

10.  Bethe,  H.A.  “The  Meson  Theory  of  Nuclear  Eorces,  General  Theory  I,”  Physical 
Review,  N 51:260-212  (Eeb  1940) 

11.  Bethe,  H.A.  “The  Meson  Theory  of  Nuclear  Eorces,  General  Theory  II,  Theory  of 
the  Deuteron”  Physical  Review  C,  V57:390-4I3  (Mar  1940) 

12.  Holz,  J.  and  W.  Glockle.  “Nucleon-Nucleon  Scattering  in  a  Time-Dependent 
Treatment,”  Physical  Review  C,  V37N4: 1386-1402  (April  1988) 

13.  Phillips,  R.J.N.  “The  Two-Nucleon  Interaction,”  Reports  Progress  Physics,  22:562- 
634(1959) 


183 


14.  J.M.  Blatt  and  L.C.  Biedenharn,  “Neutron-proton  scattering  with  spin-orbit 
coupling  (I),  General  expressions,”  Physical  Review,  86,  399  (1952) 

15.  Pieper,  Steven  C.,  K.Varga,  and  R.B.  Wiringa.  “Quantum  Monte  Carlo  calculations 
of  A=9,10  nuclei,”  Physical  Review  C,  V66,  2002. 

16.  Frank,  W.  M.,  D.  J.  Land,  and  R.  M.  Spector.  “Singular  Potentials”  Reviews  of 
Modern  Physics,  V43  Nl,  1971. 

17.  Schweiger  W,  W  Plessas,  L.  P.  Kok  and  H.  van  Haeringen.  “Coulomb  Plus 
Separable  Potential  in  Coupled  Channel,”  Physical  Review  C,  V28  N3,  1983. 

18.  H.P.  Stapp,  T.J.  Ypsilantis,  andN.  Metropolis.  “Phase-shift  analysis  of310MeF 
proton-proton  scattering  QycpQvimQnis,"  Physical  Review,.  105,  302  (1957). 

19.  Wiringa,  R.  B.  and  S.  Pieper.  “Evolution  of  Nuclear  Spectra  with  Nuclear  Forces,” 
Physical  Review  Letters,  V89,  N18,  2002. 

20.  Carlson,  J.,  V.R.  Pandharipande,  and  R.B.  Wiringa.  “Three-Nucleon  Interaction  in 
3-  and  4-  and  go  -Body  Systems.”  Nuclear  Physics  A,  401,  59-85,  1983. 

21 .  Pandharipande,  V.R.  “Nuclear  Forces  and  Nuclear  Structure,”  Nuclear  Physics  A, 
V738;66-72  (2004) 

22.  Pieper,  Steven  C.  and  R.B  Wiringa.  “Quantum  Monte  Carlo  Calculations  of  Fight 
Nuclei,”  Annual  Review  of  Particle  Science,  V5 1 :53-90  (2001) 

23.  V.G.J.  Stoks,  R.A.M.  Klomp,  M.C.M.  Rentmeester,  and  J.J.  de  Swart.  “Partial 
Wave  Analysis  of  All  Nucleon-Nucleon  Scattering  Data  Below  350  MeV,” 
Physical  Review  C  48  (1993),  792 

24.  Th.A.  Rijken  and  V.G.J.  Stoks.  “Soft  Two-Meson-Exchange  Nucleon-Nucleon 
Potentials.  I.  One  and  Two  Pair  Diagrams,”  Physical  Review  C,  54  (1996),  2851. 

25.  Th.A.  Rijken  and  V.G.J.  Stoks.  “Soft  Two-Meson-Exchange  Nucleon-Nucleon 
Potentials.  II.  One  and  Two  Pair  Diagrams,”  Physical  Review  C,  54  (1996),  2869 

26.  Yukawa,  H.  “On  the  Interaction  of  Elementary  Particles  I,”  Proceedings  of  the 
Physics-Math  Society  Japan,  17  (1935),  48-57. 

27.  J.R.  Bergervoet,  P.C.  van  Campen,  Klomp,  R.,  de  Kok,  J.,  Rijken,  T.A.,  Stoks, 

V.G.  J.,  and  J.J.  de  Swart.  “Phase  Shift  Analysis  of  all  proton-proton  scattering  data 
below  TLab=350MeV.,”R/2yVca/ Review  C,  V41,  1435-1451  (1990) 

28.  Roy,  R.R.and  B.P.  Nigam.  Nuclear  Physics:  Theory  and  Experiment.  New  York: 
John  Wiley  and  Sons,  1967.  pp.  46-140,  563-583. 


184 


29.  Preston,  M.A.  and  R.K.  Bhaduri.  Structure  of  the  Nucleus.  Massachusetts: 
Addison  Wesley,  1982.  pp.  3-53,  158-166. 

30.  Wong,  Samuel  S.  M.  Introductory  Nuclear  Physics.  New  York:  John  Wiley  and 
Sons,  1998.  pp.  21-102,  409-434. 

31.  Hans,  H.  S.  Nuclear  Physics:  Experimental  and  Theoretical.  New  Delhi:  New 
Age,  2001.  pp.  175-189 

32.  Pashinin,  V  V.G.J.  Stoks,  R.A.M.  Klomp,  C.P.F.  Terheggen,  and  J.J.  de  Swart. 
“Construction  of  high-quality  NN  potential  models,”  Physical  Review  C,  49,  2950 
(1994). 

33.  Machleidt,  R.  “High  Precision,  Charge  Dependent  Bonn  Nucleon-Nucleon 
yiod^V  Physical  Review  C,  V63,  1-32  (2001) 

34.  Kroger,  H.  “Time  Dependent  Methods  in  Scattering  Theory,”  Physics  Reports, 
V210  45-109,  1992. 

35.  Eisenbud,  Leonard.  “Time  Dependent  Formulation  of  the  Collision  Problem  in 
Quantum  Mechanics,”  Journal  of  the  Franklin  Institute,  261,  409-420  (1956) 

36.  Weeks,  David  E.  and  David  J.  Tannor.  “A  time-dependent  formulation  of  the 
scattering  matrix  using  Moller  operators,”  Chemical  Physics  Letters,  207:  301-308 
(May  1993) 

37.  Tannor,  David  J.  and  David  E.  Weeks.  “Wave  packet  correlation  function 
formulation  of  scattering  theory:  The  quantum  analog  of  classical  S-matrix 
theory,”  The  Journal  of  Chemical  Physics,  98:  3884-3893  (March  1993) 

38.  Weeks,  David  E.  and  David  J.  Tannor.  “A  time-dependent  formulation  of  the 
scattering  matrix  form  the  collinear  reaction  H  +  Hfu)  ^  Hfu')  + H  ,”  Chemical 
Physics  Letters,  224:  451-458  (July  1994) 

39.  Calfas,  R.S.  and  David  E.  Weeks.  “A  new  application  of  absorbing  boundary 
conditions  for  computing  collinear  quantum  reactive  scattering  matrix  elements,” 
Chemical  Physics  Letters,  263:  292-296  (December  1996) 

40.  Weeks,  David  E.  and  Roy  S.  Calfas,  Air  Force  Technical  Report,  in  press. 

41.  Niday,  Thomas  A.  “Inelastic  Scattering  Matrix  Elements  for  the  Collision 

+  (7)  ^^(^^3/2) +  -^2(70’”  Air  Force  Institute  of  Technology  (AU), 

Wright-Patterson  AFB  OH,  1999.  Thesis  AFIT/GAP/ENP/99M  06 

42.  Leforestier,  C.,  R.  H.  Bisseling,  C.  Cerjan,  M.D.  Feit,  R.  Friesner,  A.  Guldberg,  A. 
Hammerisch,  G.  Jolicard,  W.  Karrlein,  H.D.  Meyer,  N.  Lipkin,  O.  Roncero  and  R. 

185 


Kosloff.  “A  Comparison  of  Different  Propagation  Sehemes  for  the  Time 
Dependent  Sehrodinger  Equation,”  The  Journal  of  Computational  Physics,  94:  59- 
80(1991) 

43.  Auerbaeh,  N,  J  Huffner,  A.K.  Kerman,  and  C.M.  Shaken.  “A  Theory  of  Analog 
Resonanees,”  Review  Modern  Physics,  V44,  48  (1972) 

44.  Andrews,  Mark.  “Total  Time  Derivatives  in  Quantum  Mechanies,”  American 
Journal  of  Physics,  V71  (4),  2003. 

45.  Althorpe,  S.  C.  and  D.  C.  Clary.  “Quantum  Seattering  Calculations  on  Chemical 
Reactions,”  Annual  Review  of  Physical  Chemistry,  54,  493-529,  2003. 

46.  Jolicard,  G.  “Effective  Hamiltonian  Theory  and  Molecular  Dynamics,”  Annual 
Review  of  Physical  Chemistry,  V46,  83-108,  1995. 

47.  Bowman,  J.  M.  and  G.  Schatz.  “Theoretical  Studies  of  Polyatomic  Bimolecular 
Reaction  Dynamics,”  Annual  Review  of  Physical  Chemistry, \A6,  169-195,  1995. 

48.  Jackson,  B.  “Time-Dependent  Wave-Packet  Approach  to  Quantum  Reactive 
Scattering,”  Hnnwa/ Review  of  Physical  Chemistry,  V46,  251-274,  1995. 

49.  Tannor,  David  J.  Introduction  to  Quantum  Mechanics:  A  Time  Dependent 
Perspective.  California:  University  Science,  2007.  pp.  559-600. 

50.  Kosloff,  Ronnie.  “Time  Dependent  Quantum-Mechanical  Methods  for  Molecular 
Dynamics,”  Journal  of  Chemical  Physics.  V92,  2087-2100,  1988 

51.  Taylor,  J.R.  Scattering  Theory:  The  Quantum  Theory  ofNonrelativistic  Collisions. 
New  York:  Krieger,  1987.  pp.  7-36,  180-185,  259-270. 

52.  Calfas,  Roy  S.  Reactive  Quantum  Scattering  in  Two  Dimensions.  AirEorce 
Institute  of  Technology  (AU),  Wright-Patterson  APB  OH,  1997.  Dissertation 
APlT/DSP/ENP/97  02 

53.  Krane,  Kenneth  S.  Introductory  Nuclear  Physics.  New  York:  John  Wiley  and 
Sons,  Inc.,  1988.  pp.  80-160. 

54.  Turner,  James  E.  Atoms,  Radiation,  and  Radiation  Protection.  New  York:  John 
Wiley  and  Sons,  1995 

55.  Shultis,  J.  Kenneth  and  Richard  E.  Paw.  Fundamentals  of  Nuclear  Science  and 
Engineering.  Plorida,  CRC  Press,  2002.  pp.  27-35. 

56.  Bridgman,  Charles  J.  Introduction  to  Nuclear  Weapons  Effects.  Defense  Threat 
Reduction  Agency,  2001 


186 


57.  Wilson,  Edmund.  An  Introduction  to  Particle  Accelerators.  New  York:  Oxford, 
2001 

58.  Lechanoine-LeLuc,  C  and  F.  Lehar.  “Nucleon-Nucleon  Elastic  Scattering  and  Total 
Cross  Sections,”  Reviews  of  Modern  Physics,  V65,  No.l  (1993) 

59.  Eisberg,  R.  and  R.  Resnick.  Quantum  Physics  of  Atoms,  Molecules,  Solids,  Nuclei, 
and  Particles.  New  York:  John  Wiley  and  Sons,  1985.  pp.  289-296,  667-707. 

60.  Frizsch,  H.  Quarks:  The  Stuff  of  Matter.  lAom  York.:  Basic,  1983 

61.  Eisenbud,  L.  and  E.P.  Wigner.  “Invariant  Forms  of  Interaction  between  Nuclear 
Va.d.iclQS,''’  Proceedings  of  National  Academy  of  Sciences,  27  (1941),  281-289 

62.  Grenier,  W.  and  J. A.  Maruhn.  Nuclear  Models .  New  York:  Springer- Verlag, 

1996. 

63.  Krohn,  V.E.  and  G.R.Ringo.  “Measurement  of  the  Electron-Neutron  Interaction 
by  the  Asymmetrical  Scattering  of  Thermal  Neutrons  by  Noble  Gases,”  Physical 
Review.  V148,  Number  4,  1303-1311,  1966 

64.  Koonin,  Stephen  E.  Computational  Physics.  California,  Benjamin  Cummings: 
1986.  pp.  55-105. 

65.  Glockle,  W.  The  Quantum  Mechanical  Few-Body  Problem.  New  York,  Springer- 
Verlag,  1983.  pp.  40-76. 

66.  Townsend,  John  S.  A  Modern  Approach  to  Quantum  Mechanics.  New  York: 
McGraw  Hill,  1992.  pp.  120-146. 

67.  Mathews,  P.M.  and  K.  Venkatesan.  A  Textbook  of  Quantum  Mechanics.  New 
Delhi:  McGraw-Hill,  1976.  pp.  237-257. 

68.  Edmonds,  A.R.  Angular  Momentum  in  Quantum  Mechanics.  New  Jersey: 
Princeton  University  Press,  1957.  pp.  122-132. 

69.  Griffiths,  David  J.  Introduction  to  Quantum  Mechanics.  New  Jersey:  Prentice 
Hall,  1982.  pp.  323-325,  352-373. 

70.  Marion,  Jerry  B.  and  Stephen  T.  Thornton.  Classical  Dynamics  of  Particles  and 
Systems.  California:  Thomson,  2004.  pp.  345-350. 

71.  Farina,  J.  E.  G.  Quantum  Theory  of  Scattering  Processes .  New  York,  Pergamon, 
1973.  pp.  29-45. 

72.  Messiah,  Albert.  Quantum  Mechanics.  New  York,  John  Wiley  and  Sons,  1970.  pp. 
303-432. 


187 


73.  Mott,  N.  and  H.  Massey.  The  Theory  of  Atomic  Collisions,  3’^'^  Edition.  London; 
Oxford,  1965.  pp.  53-68. 

74.  Devanathan,  Varadaraian.  Quantum  Mechanics.  Middlesex:  Alpha  Science,  2005. 
pp.  221-224. 

75.  Levine,  Raphael  D.  Molecular  Reaction  Dynamics.  New  York:  Cambridge,  2005. 
pp.  49-51. 

76.  Abramowitz  ,  M.,  and  1.  Stegun.  Handbook  of  Mathematical  Functions  with 
Formulas,  Graphs,  and  Mathematical  Tables.  New  York,  Dover,  1964 

77.  Hulthen,  L.  and  M.  Sugawara.  Handbuch  der  Physik,  39,  1957. 

78.  Cohen-Tannoudji,  C.,  B.  Din.,  and  F.  Laloe.  Quantum  Mechanics  (2  Volumes). 
New  York:  John  Wiley  and  Sons,  1977.  pp.  903-964. 

79.  Private  Communication  with  R.  B.  Wiringa 

80.  Pieper,  Steven  C.  “Quantum  Monte  Carlo  Calculations  of  Light  Nuclei,”  Nuclear 
Physics  A,  V75 1:5 16-532  (2005) 

81.  Merzbacher,  Eugene.  Quantum  Mechanics.  New  York:  John  Wiley  and  Sons,  Inc., 
1997.  pp.  35-53,  139-191,  215-251. 

82.  Hamming,  R.W.  Numerical  Methods  for  Scientists  and  Engineers.  New  York: 
Dover  Publications,  Inc.,  1973.  pp.  382-387,  503-563. 

83.  Burden,  Richard  L.  and  J.  Douglas  Faires.  Numerical  Analysis .  Pacific  Grove: 
Brooks/Cole  Publishing  Company,  1997.  pp.  188-201,  249-342. 

84.  Kameyama,  H.,  M.  Kamimura.,  Y.  Fukushima.  “Coupled-Rearrangement-Channel 
Gaussian-Basis  Variational  Method  for  Tri-Nucleon  States,”  Physical  Review  C, 
V40,N2:974-987  (August  1989) 

85.  Taylor,  J.R.  An  Introduction  to  Error  Analysis.  Mill  Valley,  Ca:  University 
Science  Books,  1982.  pp.  45-79. 

86.  Valderrama,  M  Pavon  and  E.Ruis  Arrioloa,  “Fow  Energy  NN  Scattering  at  Next- 
to-Next-to-Next-to-Next-to-leading  order  for  partial  waves  with  j  <5 ,”  Physical 
Review  C,  72,  2005 

87.  Stoks,  V.G.  J.  and  J.J.  de  Swart.  “Comparison  of  Potential  Models  with  PP 
Scattering  Data  Below  350  MeVf  Physical  Review  C,  V47,  761-767  (1993) 

88.  Sitenko,  A  and  V  Tartakoovskii,  Theory  of  Nucleus,  Nuclear  Structure  and 
Interaction,  Boston,  Kluwer,  1997.  pp.  1-89. 

188 


89.  Rodberg,  Leonard  S.  and  R.  M.  Thaler.  Introduction  to  the  Quantum  Theory  of 
Scattering,  New  York,  Aeademie  Press,  1967.  pp.  63-73 

90.  Austen,  G.J.M.  and  J.J.  de  Swart.  “Improved  Coulomb  potential”.  Physical  Review 
Letters,  50,  2039  (1983). 

91.  J.R.  Bergervoet,  P.C.  van  Campen,  W.A.  van  der  Sanden,  and  J.J.  de  Swart. 

“Phase  shift  analysis  of  0-30  MeV  pp  scattering  data,” 

Physical  Review  C,  38,  15  (1988). 

92.  Suslov,  V.M.  and  B.  Vlahovic.,  “Generalization  of  the  Numerov  method  for 
solution  of  Nd  breakup  problem  in  configuration  space”.  Physical  Review  C,  V69, 
2004 

93.  Burcham,  W.  E.  Elements  of  Nuclear  Physics,  New  York,  Longman,  1979. 

94.  Breit,  G.  and  R.D.  Haracz.  Nucleon-Nucleon  Scattering,  High  Energy  Physics 
Volume  1.  New  York,  Academic  Press,  1967. 

95.  MacGregor,  M.  H.,  M.  J.  Moravcsik,  and  H.P.  Stapp.  “Nucleon-Nucleon  Scattering 
Experiments  and  their  Phenomenological  Analysis,”  Annual  Review  of  Nuclear 
Science,  VIO,  1960. 

96.  Glockle,  W.  H.  Witala,  D.  Hiiber,  H.  Kamada,  and  J.  Golak.  “The  Three  Nucleon 
Continuum:  Achievements,  Challenges,  and  Applications,”  Physics  Reports,  274, 
107-285,  1996. 

97.  Carlson  J.  and  R.  Schiavilla.  “Structure  and  Dynamics  of  Pew  Nucleon  Systems,” 
Reviews  of  Modern  Physics.,  V70  N3,  1998. 

98.  Paddeev,  L.  D.  Sov.  Phys.  JETP,  12,  1014,  1961. 

99.  Chen,  Z.  M.,  W.  Tornow,  and  A.  Kievsky.  “Extension  of  the  Proton-Deuteron 
Phase-Shift  Analysis  to  E'p=22.7  MeV  and  "^Pj  Phase  Shifts,”  Few  Body  Systems, 
V35,  15-31,2004 

100.  Tomow,  W.,  A.  Kievsky.,  and  H.  Witala.  “Improved  Proton-Deuteron  Phase-Shift 
Analysis  Above  the  Deuteron  Breakup  Threshold  and  the  Three-Nucleon 
Analyzing-Power  Puzzle,”  Few  Body  Systems,  V32,  53-81,  2002 

101.  Huber,  D.,  W.  Glockle,.  J.  Golak.,  H.  Witala.,  H.  Kamada.,  A.  Kievsky.,  S. 

Rosati.,  and  M.  Viviani.  “Realistic  Phase  Shift  and  Mixing  Parameters  for  Elastic 
neuteron-deuteron  scattering:  Comparison  of  Momentum  Space  and  Configuration 
Space  Methods,”  Physical  Review  C,  V51,  Number  3,  1100-1107,  1995 


189 


102.  Kievsky,  A.,  M.  Viviani S.  Rosati.,  D.  Hiiber.,  W.  Glockle,.  H.  Kamada.,  H. 
Witala.,  and  J.  Golak.  “Benchmark  Calculations  for  Polarization  Observables  in 
Three-Nucleon  Scattering,”  Physical  Review  C,  V58,  Number  6,  3085-3092,  1998 

103.  Breit,  G.  “Relativistic  Corrections  for  High-Energy  p-p  Scattering,”  Physical 
Review,  99,  1581  (1955). 

104  Fujita,  J.  and  H.  Miyazawa.  “Pion  Theory  of  Three-Body  Forces,”  Progress  of 
Theoretical  Physics,  VI 7:360-365  (1957) 

105.  Feit  M.D.,  J.A.  Fleck,  Jr.,  and  A.  Steiger.  “Solution  of  the  Schrodinger  equation  by 
a  spectral  method,”  Journal  of  Computational  Physics,  47,  412  (1982) 

106.  Foveland,  Walter  D.,  David  J.  Morrissey,  and  Glenn  T.  Seaborg  Modern  Nuclear 
Chemistry.  New  Jersey:  John  Wiley  and  Sons,  2006 


190 


Vita 


Major  Brian  S.  Davis  graduated  from  Ware  Shoals  High  School  in  Ware  Shoals, 
South  Carolina  in  1987.  After  having  served  in  the  United  States  Air  Force  as  a  Photo- 
Sensor  Maintenance  Specialist  at  RAF  Bentwaters,  UK,  Major  Davis  entered 
undergraduate  study  at  Clemson  University  where  he  majored  in  Physics  and  graduated 
in  1996.  He  enrolled  in  the  graduate  physics  program  in  1999  and  was  awarded  a  Master 
of  Science  in  Physics  from  Wright  State  University  in  2002.  Major  Davis  was 
commissioned  through  AFROTC  Detachment  770  at  Clemson  University  in  1996. 

His  first  officer  assignment  in  1996  was  at  Holloman  AFB  where  he  performed 
the  duties  of  Propulsion  and  Munitions  Flight  Commander.  After  a  three-year  tour,  he 
worked  as  a  bio-acoustic  research  physicist  assigned  to  the  Human  Effectiveness 
Directorate,  Air  Force  Research  Laboratory,  Wright-Patterson  AFB  Ohio.  In  2002,  he 
entered  the  Nuclear  Physics  program  of  the  Department  of  Engineering  Physics,  Graduate 
School  of  Engineering  Management,  Air  Eorce  Institute  of  Technology.  Erom  2006  to 
present,  he  has  served  as  Sciences  Support  Division  Deputy,  Materials  Directorate,  Air 
Eorce  Technical  Application  Center  at  Patrick  APB. 


191 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  074-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  the  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  an  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a 
currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

25-03-2010  Doctoral  Dissertation  Sep  2002  -  Mar  2010 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


Time  Dependent  Channel  Packet  Calculation  of  Two  Nucleon  Scattering  Matrix  Elements 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Davis,  Brian  S.,  Major,  USAF 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(S) 

Air  Force  Institute  of  Technology 

Graduate  School  of  Engineering  and  Management  (AFIT/EN) 

2950  Hobson  Way 
WPAFB  OH  45433-7765 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Dr.  Charles  Brennan 

Air  Force  Technical  Applications  Center 

1030  South  Highway  AlA 

Patrick  AFB,  FI  32925 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


51.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/DS/ENP/10-M03 


10.  SPONSOR/MONITOR’S 
ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

A  new  approach  to  calculating  nucleon-nucleon  scattering  matrix  elements  using  a  proven  atomic  time-dependent  wave  packet  technique  is 
investigated.  Wave  packets  containing  centripetal  barrier  information  are  prepared  in  close  proximity  to  nuclear  well.  This  is  accomplished  by  first  using  an 
analytic  equation  to  determine  the  wave  packets  in  a  suitable  intermediate  asymptotic  state  where  the  centripetal  barrier  is  negligible.  Then,  the  split  operator 
technique  is  used  to  propagate  the  wave  packets  back  to  their  original  positions  under  the  lull  Hamiltonian.  Here,  one  wave  packet  is  held  stationary  while  the 
other  is  allowed  to  evolve  and  explore  the  nuclear  well.  Scattering  matrix  elements  are  computed  from  the  correlation  function  between  the  stationary  wave 
packet  and  the  evolving  wave-packet  after  it  has  interacted  with  the  nuclear  potential.  Determination  of  nucleon-nucleon  phase  shifts  follows  directly  from 
computation  of  the  scattering  matrix  elements.  This  technique  is  ideally  suited  for  determining  nuclear  scattering  matrix  elements  and  phase  shifts  as  it 
provides  a  high  degree  of  energy  resolution  with  lower  computational  effort  than  traditional  time  independent  methods.  These  advantages  will  lead  to  a 
greater  understanding  of  reactions  involving  nucleons  with  other  elementary  particles. 


15.  SUBJECT  TERMS 

Nucleon-Nucleon,  Scattering  Matrix,  Time  Dependent ,  Phase  Shift,  Wave  Packet,  Split  Operator 


16.  SECURITY  CLASSIFICATION 
OF: 


17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER 
OF 

PAGES 

206 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Dr.  David  E  Weeks,  AFIT/ENP 


19b.  TELEPHONE  NUMBER  (Include  area  code) 
(937)  255-3636,  ext  4561 ;  e-mail:  david.weeks@afit.edu 


Standard  Form  298  (Rev:  8-98) 

PiEsciibadbyANSISldZ39-18 


